Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,614,214 T→G N164H (AAT→CAT)  sad ← succinate semialdehyde dehydrogenase, NAD(P)+‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,614,2140TG98.4% 244.0 / ‑5.2 63N164H (AAT→CAT) sadsuccinate semialdehyde dehydrogenase, NAD(P)+‑dependent
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (34/28);  total (34/29)
Fisher's exact test for biased strand distribution p-value = 4.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

CCGGTCACCGTGACAGCAGCAATGCGCGAGTCTTTAATCATCTGACTGACACCGTCGTTGTCGGCATTCAGCCAGCCATATACGCCTTGTGGGATACCCGCATCTTTAAACACCTGGGCAATGAGCTGTGCACAGCCCATCACATTCGGCGCATGTTTAAGTAAGTAGCCGTTACCTGCAAGAATGATGGGAACAGCGCCACGCATCACCTGCCATAACGGAAAATTCCACGGCATAATCGCCAGAATCGTCCCCAACGGTCGATACTCAATAACCGCCTGCTGATTTTCC  >  NC_000913/1614069‑1614359
                                                                                                                                                 |                                                                                                                                                 
ccGGTCACCGTGACAGCGGCAATGCGCGAGTCTTTAATCATCTGACTGACACCGTCGTTGTCGGCATTCAGCCAGCCATATACGCCTTGTGGGATACCCGCATCTTTAAACACCTGGGCAATGAGCTGTGCACAGCCCATCACATGCGgc                                                                                                                                               <  1:207284/150‑1 (MQ=255)
 cGGTCACCGTGACAGCAGCAATGCGCGAGTCTTTAATCATCTGACTGACACCGTCGTTGTCGGCATTCAGCCAGCCATATACGCCTTGTGGGATACCCGCATCTTTAAACACCTGGGCAATGAGCTGTGCACAGCCCATCACATGCGgcgc                                                                                                                                             >  1:35968/1‑151 (MQ=255)
        cGTGACAGCAGCAATGCGCGAGTCTTTAATCATCTGACTGACACCGTCGTTGTCGGCATTCAGCCAGCCATATACGCCTTGTGGGATACCCGCATCTTTAAACACCTGGGCAATGAGCTGTGCACAGCCCATCACATGCGGCGCATGttt                                                                                                                                       >  1:1085941/1‑150 (MQ=255)
           gACAGCAGCAATGCGCGAGTCTTTAATCATCTGACTGACACCGTCGTTGTCGGCATTCAGCCAGCCATATACGCCTTGTGGGATACCCGCATCTTTAAACACCTGGGCAATGAGCTGTGCACAGCCCATCACATTCGGCGCATGTTtaag                                                                                                                                    <  2:637432/150‑1 (MQ=255)
            aCAGCAGCAATGCGCGAGTCTTTAATCATCTGACTGACACCGTCGTTGTCGGCATTCAGCCAGCCATATACGCCTTGTGGGATACCCGCATCTTTAAACACCTGGGCAATGAGCTGTGCACAGCCCATCACATGCGGCGCATGTTtaagt                                                                                                                                   >  2:139423/1‑150 (MQ=255)
            aCAGCAGCAATGCGCGAGTCTTTAATCATCTGACTGACACCGTCGTTGTCGGCATTCAGCCAGCCATATACGCCTTGTGGGATACCCGCATCTTTAAACACCTGGGCAATGAGCTGTGCACAGCCCATCACATGCGGCGCATGTTtaagt                                                                                                                                   >  1:904092/1‑150 (MQ=255)
             cagcagCAATGCGCGAGTCTTTAATCATCTGACTGACACCGTCGTTGTCGGCATTCAGCCAGCCATATACGCCTTGTGGGATACCCGCATCTTTAAACACCTGGGCAATGAGCTGTGCACAGCCCATCACATGCGGCGCATGTTtaagtaa                                                                                                                                 <  2:363587/151‑1 (MQ=255)
             cagcagCAATGCGCGAGTCTTTAATCATCTGACTGACACCGTCGTTGTCGGCATTCAGCCAGCCATATACGCCTTGTGGGATACCCGCATCTTTAAACACCTGGGCAATGAGCTGTGCACAGCCCATCACATGCGGCGCATGTTtaagta                                                                                                                                  >  2:121679/1‑150 (MQ=255)
                cagcaATGCGCGAGTCTTTAATCATCTGACTGACACCGTCGTTGTCGGCATTCAGCCAGCCATATACGCCTTGTGGGATACCCGCATCTTTAAACACCTGGGCAATGAGCTGTGCACAGCCCATCACATGCGGCGCATGTTTaagtaagt                                                                                                                               >  1:820785/1‑150 (MQ=255)
                 agcaATGCGCGAGTCTTTAATCATCTGACTGACACCGTCGTTGTCGGCATTCAGCCAGCCATATACGCCTTGTGGGATACCCGCATCTTTAAACACCTGGGCAATGAGCTGTGCACAGCCCATCACATGCGGCGCATGTTTAAGTAAGTAg                                                                                                                             <  2:966083/151‑1 (MQ=255)
                 agcaATGCGCGAGTCTTTAATCATCTGACTGACACCGTCGTTGTCGGCATTCAGCCAGCCATATACGCCTTGTGGGATACCCGCATCTTTAAACACCTGGGCAATGAGCTGTGCACAGCCCATCACATGCGGCGCATGTTTAAGTAAGTAg                                                                                                                             <  2:211217/151‑1 (MQ=255)
                    aaTGCGCGAGTCTTTAATCATCTGACTGACACCGTCGTTGTCGGCATTCAGCCAGCCATATACGCCTTGTGGGATACCCGCATCTTTAAACACCTGGGCAATGAGCTGTGCACAGCCCATCACATGCGGCGCATGTTTAAGTAAGTAGCCg                                                                                                                          >  2:415708/1‑151 (MQ=255)
                       gcgcGAGTCTTTAATCATCTGACTGACACCGTCGTTGTCGGCATTCAGCCAGCCATATACGCCTTGTGGGATACCCGCATCTTTAAACACCTGGGCAATGAGCTGTGCACAGCCCATCACATGCGGCGCATGTTTAAGTAAGTAGCCGtt                                                                                                                        <  1:37842/150‑1 (MQ=255)
                       gcgcGAGTCTTTAATCATCTGACTGACACCGTCGTTGTCGGCATTCAGCCAGCCATATACGCCTTGTGGGATACCCGCATCTTTAAACACCTGGGCAATGAGCTGTGCACAGCCCATCACATGCGGCGCATGTTTAAGTAAGTAGCCGTTa                                                                                                                       <  1:659189/151‑1 (MQ=255)
                         gcgAGTCTTTAATCATCTGACTGACACCGTCGTTGTCGGCATTCAGCCAGCCATATACGCCTTGTGGGATACCCGCATCTTTAAACACCTGGGCAATGAGCTGTGCACAGCCCATCACATGCGGCGCATGTTTAAGTAAGTAGCCGTTAcc                                                                                                                     >  2:865666/1‑151 (MQ=255)
                           gAGTCTTTAATCATCTGACTGACACCGTCGTTGTCGGCATTCAGCCAGCCATATACGCCTTGTGGGATACCCGCATCTTTAAACACCTGGGCAATGAGCTGTGCACAGCCCATCACATGCGGCGCATGTTTAAGTAAGTAGCCGTTACCTg                                                                                                                   >  2:1083038/1‑151 (MQ=255)
                            aGTCTTTAATCATCTGACTGACACCGTCGTTGTCGGCATTCAGCCAGCCATATACGCCTTGTGGGATACCCGCATCTTTAAACACCTGGGCAATGAGCTGTGCACAGCCCATCACATGCGGCGCATGTTTAAGTAAGTAGCCGTTACCTGc                                                                                                                  <  2:864624/151‑1 (MQ=255)
                             gTCTTTAATCATCTGACTGACACCGTCGTTGTCGGCATTCAGCCAGCCATATACGCCTTGTGGGATACCCGCATCTTTAAACACCTGGGCAATGAGCTGTGCACAGCCCATCACATGCGGCGCATGTTTAAGTAAGTAGCCGTTACCTGc                                                                                                                  >  2:565443/1‑150 (MQ=255)
                                tttAATCATCTGACTGACACCGTCGTTGTCGGCATTCAGCCAGCCATATACGCCTTGTGGGATACCCGCATCTTTAAACACCTGGGCAATGAGCTGTGCACAGCCCATCACATGCGGCGCATGTTTAAGTAAGTAGCCGTTACCTGCa                                                                                                                 >  2:1217/1‑148 (MQ=255)
                                  tAATCATCTGACTGACACCGTCGTTGTCGGCATTCAGCCAGCCATATACGCCTTGTGGGATACCCGCATCTTTAAACACCTGGGCAATGAGCTGTGCACAGCCCATCACATGCGGCGCATGTTTAAGTAAGTAGCCGTTACCTGCAAGAat                                                                                                            >  1:561533/1‑151 (MQ=255)
                                   aaTCATCTGACTGACACCGTCGTTGTCGGCATTCAGCCAGCCATATACGCCTTGTGGGATACCCGCATCTTTAAACACCTGGGCAATGAGCTGTGCACAGCCCATCACATGCGGCGCATGTTTAAGTAAGTAGCCGTTACCTGCAAGAatg                                                                                                           >  1:143804/1‑151 (MQ=255)
                                   aaTCATCTGACTGACACCGTCGTTGTCGGCATTCAGCCAGCCATATACGCCTTGTGGGATACCCGCATCTTTAAACACCTGGGCAATGAGCTGTGCACAGCCCATCACATGCGGCGCATGTTTAAGTAAGTAGCCGTTACCTGCAAGAatg                                                                                                           <  2:994468/151‑1 (MQ=255)
                                    atcatcTGACTGACACCGTCGTTGTCGGCATTCAGCCAGCCATATACGCCTTGTGGGATACCCGCATCTTTAAACACCTGGGCAATGAGCTGTGCACAGCCCATCACATGCGGCGCATGTTTAAGTAAGTAGCCGTTACCTGCAAGAatg                                                                                                           >  2:409794/1‑150 (MQ=255)
                                           gactgacACCGTCGTTGTCGGCATTCAGCCAGCCATATACGCCTTGTGGGATACCCGCATCTTTAAACACCTGGGCAATGAGCTGTGCACAGCCCATCACATGCGGCGCATGTTTAAGTAAGTAGCCGTTACCTGCAAGAATGATGGGaa                                                                                                    <  2:722341/150‑1 (MQ=255)
                                           gacgggcacCGTCGTTGTCGGCATTCAGCCAGCCCTATACGCCTTGTGGGATACCCGCATCTTTAAACTCCTGGGCAATGAGCTGTGCACAGCCCATCACATGCGGCGCATGTTTAAGTAAGTAGCCGTTACCTGCAAGAATGATggg                                                                                                      <  2:650280/142‑1 (MQ=255)
                                            actgacACCGTCGTTGTCGGCATTCAGCCAGCCATATACGCCTTGTGGGATACCCGCATCTTTAAACACCTGGGCAATGAGCTGTGCACAGCCCATCACATGCGGCGCATGTTTAAGTAAGTAGCCGTTACCTGCAAGAATGATGGGaa                                                                                                    >  1:3165/1‑149 (MQ=255)
                                               gccacCGTCGTTGTCGGCATTCAGCCAGCCATATCCGCCTTGTGGGATACCCGCATCTTTAAACACCTGGGCAATGAGCTGTGCACAGCCCATCACATGCGGCGCATGTTTAAGTAAGTAGCCGTTACCTGCAAGAATGATGGGAACAgc                                                                                                <  2:647686/148‑1 (MQ=255)
                                                        gtTGTCGGCATTCAGCCAGCCATATACGCCTTGTGGGATACCCGCATCTTTAAACACCTGGGCAATGAGCTGTGCACAGCCCATCACATGCGGCGCATGTTTAAGTAAGTAGCCGTTACCTGCAAGAATGATGGGAACAGCGCCACGCATc                                                                                      >  1:833829/1‑151 (MQ=255)
                                                          tGTCGGCATTCAGCCAGCCATATACGCCTTGTGGGATACCCGCATCTTTAAACACCTGGGCAATGAGCTGTGCACAGCCCATCACATGCGGCGCATGTTTAAGTAAGTAGCCGTTACCTGCAAGAATGATGGGAACAGCGCCACGCATCa                                                                                     <  2:21768/150‑1 (MQ=255)
                                                           gTCGGCATTCAGCCAGCCATATACGCCTTGTGGGATACCCGCATCTTTAAACACCTGGGCAATGAGCTGTGCACAGCCCATCACATGCGGCGCATGTTTAAGTAAGTAGCCGTTACCTGCAAGAATGATGGGAACAGCGCCACGCATCAc                                                                                    <  2:752754/150‑1 (MQ=255)
                                                               gCATTCAGCCAGCCATATACGCCTTGTGGGATACCCGCATCTTTAAACACCTGGGCAATGAGCTGTGCACAGCCCATCACATGCGGCGCATGTTTAAGTAAGTAGCCGTTACCTGCAAGAATGATGGGAACAGCGCCACGCATCACCTGcc                                                                               <  2:585886/151‑1 (MQ=255)
                                                                cATTCAGCCAGCCATATACGCCTTGTGGGATACCCGCATCTTTAAACACCTGGGCAATGAGCTGTGCACAGCCCATCACATGCGGCGCATGTTTAAGTAAGTAGCCGTTACCTGCAAGAATGATGGGAACAGCGCCACGCATCACCTGcc                                                                               >  2:9940/1‑150 (MQ=255)
                                                                         agccaTATACGCCTTGTGGGATACCCGCATCTTTAAACACCTGGGCAATGAGCTGTGCACAGCCCATCACATGCGGCGCATGTTTAAGTAAGTAGCCGTTACCTGCAAGAATGATGGGAACAGCGCCACGCATCACCTGCCATAACGGaa                                                                      >  2:583281/1‑150 (MQ=255)
                                                                          gccaTATACGCCTTGTGGGATACCCGCATCTTTAAACACCTGGGCAATGAGCTGTGCACAGCCCATCACATGCGGCGCATGTTTAAGTAAGTAGCCGTTACCTGCAAGAATGATGGGAACAGCGCCACGCATCACCTGCCATAACGGaaa                                                                     <  2:210516/150‑1 (MQ=255)
                                                                           ccaTATACGCCTTGTGGGATACCCGCATCTTTAAACACCTGGGCAATGAGCTGTGCACAGCCCATCACATGCGGCGCATGTTTAAGTAAGTAGCCGTTACCTGCAAGAATGATGGGAACAGCGCCACGCATCACCTGCCATAACGGAAAAt                                                                   <  2:791547/151‑1 (MQ=255)
                                                                             atatACGCCTTGTGGGATACCCGCATCTTTAAACACCTGGGCAATGAGCTGTGCACAGCCCATCACATGCGGCGCATGTTTAAGTAAGTAGCCGTTACCTGCAAGAATGATGGGAACAGCGCCACGCATCACCTGCCATAACGGAAAATTc                                                                 >  2:968040/1‑151 (MQ=255)
                                                                              tataCGCCTTGTGGGATACCCGCATCTTTAAACACCTGGGCAATGAGCTGTGCACAGCCCATCACATGCGGCGCATGTTTAAGTAAGTAGCCGTTACCTGCAAGAATGATGGGAACAGCGCCACGCATCACCTGCCATAACGGAAAATTcc                                                                >  2:893228/1‑151 (MQ=255)
                                                                                             aTACCCGCATCTTTAAACACCTGGGCAATGAGCTGTGCACAGCCCATCACATGCGGCGCATGTTTAAGTAAGTAGCCGTTACCTGCAAGAATGATGGGAACAGCGCCACGCATCACCTGCCATAACGGAAAATTCCACGGCATAATCGCCa                                                 <  1:387547/151‑1 (MQ=255)
                                                                                              tACCCGCATCTTTAAACACCTGGGCAATGAGCTGTGCACAGCCCATCACATGCGGCGCATGTTTAAGTAAGTAGCCGTTACCTGCAAGAATGATGGGAACAGCGCCACGCATCACCTGCCATAACGGAAAATTCCACGGCATAATCGcc                                                  >  1:183783/1‑149 (MQ=255)
                                                                                                 ccGCATCTTTAAACACCTGGGCAATGAGCTGTGCACAGCCCATCACATGCGGCGCATGTTTAAGTAAGTAGCCGTTACCTGCAAGAATGATGGGAACAGCGCCACGCATCACCTGCCATAACGGAAAATTCCACGGCATAATCGCCAGAAt                                             <  2:16744/151‑1 (MQ=255)
                                                                                                 ccGCATCTTTAAACACCTGGGCAATGAGCTGTGCACAGCCCATCACATGCGGCGCATGTTTAAGTAAGTAGCCGTTACCTGCAAGAATGATGGGAACAGCGCCACGCATCACCTGCCATAACGGAAAATTCCACGGCATAATCGCCAGAAt                                             <  2:594182/151‑1 (MQ=255)
                                                                                                        tttAAACACCTGGGCAATGAGCTGTGCACAGCCCATCACATGCGGCGCATGTTTAAGTAAGTAGCCGTTACCTGCAAGAATGATGGGAACAGCGCCACGCATCACCTGCCATAACGGAAAATTCCACGGCATAATCGCCAGAATCGTcccc                                      >  1:528442/1‑151 (MQ=255)
                                                                                                         ttAAACACCTGGGCAATGAGCTGTGCACAGCCCATCACATGCGGCGCATGTTTAAGTAAGTAGCCGTTACCTGCAAGAATGATGGGAACAGCGCCACGCATCACCTGCCATAACGGAAAATTCCACGGCATAATCGCCAGAATCGTcccc                                      >  2:183311/1‑150 (MQ=255)
                                                                                                         ttAAACACCTGGGCAATGAGCTGTGCACAGCCCATCACATGCGGCGCATGTTTAAGTAAGTAGCCGTTACCTGCAAGAATGATGGGAACAGCGCCACGCATCACCTGCCATAACGGAAAATTCCACGGCATAATCGCCAGAATCGTcc                                        <  1:1102659/148‑1 (MQ=255)
                                                                                                            aaCACCTGGGCAATGAGCTGTGCACAGCCCATCACATGCGGCGCATGTTTAAGTAAGTAGCCGTTACCTGCAAGAATGATGGGAACAGCGCCACGCATCACCTGCCATAACGGAAAATTCCACGGCATAATCGCCAGAATCGTCCCCAACg                                  >  2:600880/1‑151 (MQ=255)
                                                                                                             acacCTGGGCAATGAGCTGTGCACAGCCCATCACATGCGGCGCATGTTTAAGTAAGTAGCCGTTACCTGCAAGAATGATGGGAACAGCGCCACGCATCACCTGCCATAACGGAAAATTCCACGGCATAATCGCCAGAATCGTCCCCAACgg                                 >  2:1020819/1‑151 (MQ=255)
                                                                                                                   gggCAATGAGCTGTGCACAGCCCATCACATGCGGCGCATGTTTAAGTAAGTAGCCGTTACCTGCAAGAATGATGGGAACAGCGCCACGCATCACCTGCCATAACGGAAAATTCCACGGCATAATCGCCAGAATCGTCCCCAACGGTCGATa                           <  1:457882/151‑1 (MQ=255)
                                                                                                                    ggCAATGAGCTGTGCACAGCCCATCACATGCGGCGCATGTTTAAGTAAGTAGCCGTTACCTGCAAGAATGATGGGAACAGCGCCACGCATCACCTGCCATAACGGAAAATTCCACGGCATAATCGCCAGAATCGTCCCCAACGGTCGTTAc                          >  1:436217/1‑151 (MQ=255)
                                                                                                                    ggCAATGAGCTGTGCACAGCCCATCACATGCGGCGCATGTTTAAGTAAGTAGCCGTTACCTGCAAGAATGATGGGAACAGCGCCACGCATCACCTGCCATAACGGAAAATTCCACGGCATAATCGCCAGAATCGTCCCCAACGGTCGATAc                          <  1:415770/151‑1 (MQ=255)
                                                                                                                       aaTGAGCTGTGCACAGCCCATCACATGCGGCGCATGTTTAAGTAAGTAGCCGTTACCTGCAAGAATGATGGGAACAGCGCCACGCATCACCTGCCATAACGGAAAATTCCACGGCATAATCGCCAGAATCGTCCCCAACGGTCGATACTc                        >  2:808070/1‑150 (MQ=255)
                                                                                                                       aaTGAGCTGTGCACAGCCCATCACATGCGGCGCATGTTTAAGTAAGTAGCCGTTACCTGCAAGAATGATGGGAACAGCGCCACGCATCACCTGCCATAACGGAAAATTCCACGGCATAATCGCCAGAATCGTCCCCAACGGTCGATACTCa                       <  2:208408/151‑1 (MQ=255)
                                                                                                                        aTGAGCTGTGCACAGCCCATCACATGCGGCGCATGTTTAAGTAAGTAGCCGTTACCTGCAAGAATGATGGGAACAGCGCCACGCATCACCTGCCATAACGGAAAATTCCACGGCATAATCGCCAGAATCGTCCCCAACGGTCGATACTCa                       >  1:135845/1‑150 (MQ=255)
                                                                                                                          gAGCTGTGCACAGCCCATCCCATGCGGCGCATGTTTAAGTAAGCAGCCGTGACCTGCAAGAATGATGGGAACAGCGCCAAGCATCACCTGCCATAACGGAAAATTCCACGGCATAATCGCCAGAATCGTCCCCAACGGTCGATACTCAATa                    <  2:573486/151‑1 (MQ=255)
                                                                                                                          gAGCTGTGCACAGCCCATCACATGCGGCGCATGTTTAAGTAAGTAGCCGTTACCTGCAAGAATGATGGGAACAGCGCCACGCATCACCTGCCATAACGGAAAATTCCACGGCATAATCGCCAGAATCGTCCCCAACGGTCGATACTCAAt                     >  2:418297/1‑150 (MQ=255)
                                                                                                                            gCTGTGCACAGCCCATCACATGCGGCGCATGTTTAAGTAAGTAGCCGTTACCTGCAAGAATGATGGGAACAGCGCCACGCATCACCTGCCATAACGGAAAATTCCACGGCATAATCGCCAGAATCGTCCCCAACGGTCGATACTCAATAAc                  <  2:693174/151‑1 (MQ=255)
                                                                                                                               gtgCACAGCCCATCACATGCGGCGCATGTTTAAGTAAGTAGCCGTTACCTGCAAGAATGATGGGAACAGCGCCACGCATCCCCTGCCATAACGGAAAATTCCACGGCATAATCGCCAGAATCGTCCCCAACGGTCGATACTCAATAACCg                <  2:772010/150‑1 (MQ=255)
                                                                                                                                  cacaGCCCATCACATGCGGCGCATGTTTAAGTAAGTAGCCGTTACCTGCAAGAATGATGGGAACAGCGCCACGCATCACCTGCCATAACGGAAAATTCCACGGCATAATCGCCAGAATCGTCCCCAACGGTCGATACTCAATAACCGCctg            >  2:1081488/1‑151 (MQ=255)
                                                                                                                                   acaGCCCATCACATGCGGCGCATGTTTAAGTAAGTAGCCGTTACCTGCAAGAATGATGGGAACAGCGCCACGCATCACCTGCCATAACGGAAAATTCCACGGCATAATCGCCAGAATCGTCCCCAACGGTCGATACTCAATAACCGcc              <  2:666414/148‑1 (MQ=255)
                                                                                                                                       cccATCACATGCGGCGCATGTTTAAGTAAGTAGCCGTTACCTGCAAGAATGATGGGAACAGCGCCACGCATCACCTGCCATAACGGAAAATTCCACGGCATAATCGCCAGAATCGTCCCCAACGGTCGATACTCAATAACCGCCTGCTGa        <  2:515112/150‑1 (MQ=255)
                                                                                                                                       cccATCACATGCGGCGCATGTTTAAGTAAGTAGCCGTTACCTGCAAGAATGATGGGAACAGCGCCACGCATCACCTGCCATAACGGAAAATTCCACGGCATAATCGCCAGAATCGTCCCCAACGGTCGATACTCAATAACCGCCTGCTGAt       >  2:478790/1‑151 (MQ=255)
                                                                                                                                        ccATCACATGCGGCGCATGTTTAAGTAAGTAGCCGTTACCTGCAAGAATGATGGGAACAGCGCCACGCATCACCTGCCATAACGGAAAATTCCACGGCATAATCGCCAGAATCGTCCCCAACGGTCGATACTCAATAACCGCCTGCTGAtt      >  2:592205/1‑151 (MQ=255)
                                                                                                                                          aTCACATGCGGCGCATGTTTAAGTAAGTAGCCGTTACCTGCAAGAATGATGGGAACAGCGCCACGCATCACCTGCCATAACGGAAAATTCCACGGCATAATCGCCAGAATCGTCCCCAACGGTCGATACTCAATAACCGCCTGCTGAtttt    >  1:748925/1‑151 (MQ=255)
                                                                                                                                            cacaTGCGGCGCATGTTTAAGTAAGTAGCCGTTACCTGCAAGAATGATGGGAACAGCGCCACGCATCACCTGCCATAACGGAAAATTCCACGGCATAATCGCCAGAATCGTCCCCAACGGTCGATACTCAATAACCGCCTGCTGATTTTcc  >  2:609526/1‑151 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
CCGGTCACCGTGACAGCAGCAATGCGCGAGTCTTTAATCATCTGACTGACACCGTCGTTGTCGGCATTCAGCCAGCCATATACGCCTTGTGGGATACCCGCATCTTTAAACACCTGGGCAATGAGCTGTGCACAGCCCATCACATTCGGCGCATGTTTAAGTAAGTAGCCGTTACCTGCAAGAATGATGGGAACAGCGCCACGCATCACCTGCCATAACGGAAAATTCCACGGCATAATCGCCAGAATCGTCCCCAACGGTCGATACTCAATAACCGCCTGCTGATTTTCC  >  NC_000913/1614069‑1614359

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: