Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,478,3240TG45.5% 137.8 / 18.3 77G33G (GGT→GGGynbBputative CDP‑diglyceride synthase
Reads supporting (aligned to +/- strand):  ref base T (4/38);  new base G (35/0);  total (39/38)
Fisher's exact test for biased strand distribution p-value = 9.07e-18
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTTATATCTGTAC  >  NC_000913/1478179‑1478469
                                                                                                                                                  |                                                                                                                                                 
gCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGt                                                                                                                                              <  1:36222/151‑1 (MQ=255)
  gATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTATCGGAGAGGAAAAGGGGGGAg                                                                                                                                                >  1:740692/1‑147 (MQ=255)
      aaTCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACa                                                                                                                                          <  1:954245/149‑1 (MQ=255)
      aaTCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACAtt                                                                                                                                        <  2:1081300/151‑1 (MQ=255)
        tCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGGTTCCGGGGCGAAAAGGGGGGGATGTTACATTAc                                                                                                                                      >  2:801675/1‑151 (MQ=255)
             gTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGGGAGAAAAGGGGGGGGGGGTACATTAAGTAtt                                                                                                                                 >  2:465730/1‑151 (MQ=255)
                 gAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTa                                                                                                                             <  1:1074219/151‑1 (MQ=255)
                     gtgCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCtg                                                                                                                         <  1:564035/151‑1 (MQ=255)
                     gtgCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCt                                                                                                                          <  2:825825/150‑1 (MQ=255)
                      tgCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTTACTGCAGGACAAGGGGGTGAGGTGAGATTACATATTctt                                                                                                                              >  1:415961/1‑143 (MQ=255)
                       gCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTTCCGGGGAGAAAAAGGGGGGGGGGTTTATTAACTATTTTTACCtgg                                                                                                                        >  1:577048/1‑150 (MQ=255)
                             tATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTggg                                                                                                                    <  2:461977/148‑1 (MQ=255)
                             tATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGc                                                                                                                  <  2:50666/150‑1 (MQ=255)
                             tATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGGGAGGGTACATTACGTATTCGTACCTGGTgggg                                                                                                                   >  2:308141/1‑149 (MQ=255)
                              aTGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCAt                                                                                                                <  2:795647/151‑1 (MQ=255)
                              aTGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTAACGGCGGGGGAAGGGGGGGAGGGTGTATTTCGTTTTCTTagacgggtgg                                                                                                                    >  2:268292/1‑138 (MQ=255)
                              aTGAGCGCTGAAATATAATGATGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACAGTTTTCTGCTATGGCGGGTACCGGAGAGGAAAGGGGGGGGGGGTAAATTAAGGATTTGTTCCTGGCgggg                                                                                                                   >  2:105649/1‑148 (MQ=255)
                                gAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCg                                                                                                              <  2:130483/151‑1 (MQ=255)
                                 aGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGGGAGAAAAGGGGGGGGGGTTACATTACGTATTCGTACCTGGTGGGGCATCGt                                                                                                             >  2:521295/1‑151 (MQ=255)
                                    gcTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGGGAGGAAAGGGGGGGGGGGTTTATTTTGTATTTGTAACTGGTggg                                                                                                                    >  1:1051092/1‑141 (MQ=255)
                                    gcTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTAACCGAGAGAAAAGGGGGGGGGGGTACATTTCGTATTTGTACCTGGGGGGGGATCg                                                                                                              >  2:299245/1‑147 (MQ=255)
                                       gAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCAtt                                                                                                         <  2:477587/149‑1 (MQ=255)
                                       gAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGATACCGGGGAGAAAAGGGGGGGGGGGTACATTACGTATTCGTAACTTGGGGGGGATGGTCAttt                                                                                                        >  1:443043/1‑150 (MQ=255)
                                         aaTATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGGGAGAAAAGGGGGGGGGGGTACATTACGTATTTGTAACTGGTGGGGGATCGTCCTTTGtt                                                                                                     >  2:1078281/1‑151 (MQ=255)
                                              aaTGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTc                                                                                                 <  2:105236/150‑1 (MQ=255)
                                              aaTGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGGGGGGTTTCATTACGTATTCGTACCTGGTGGGGGATCGTCCTtttttttt                                                                                                  >  2:193586/1‑149 (MQ=255)
                                                   tGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGGGAGGGTTCATTACGTATTTGTACCTTGTGGGGGATCGTCATTTGttttt                                                                                                  >  1:94317/1‑144 (MQ=255)
                                                       aaaaTCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTTACGGAGGGAAAAGGGGGGGAGGGTACATTTCGTATTCGTACCTGGGGGGGGATCGTCCTTTTTTTTTTTATGGTggt                                                                                        >  2:977164/1‑148 (MQ=255)
                                                       aaaaTCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTGCCGGAGAGGAAAAGGGGGGAGGGTACATTACGTCTTTGTTGCTGGTGGGGGATCGTCCTtttttttttacttgggagt                                                                                       >  2:347919/1‑140 (MQ=255)
                                                       aaaaTCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGGGGGAAAAGGGGGGGAGGGTACATTTCGTATTCGTACCTGGGGGGGCATCGGCCTTTTTTTTTTAATGGGGAtt                                                                                       >  2:104295/1‑151 (MQ=255)
                                                       aaaaTCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTAACGGAGGGAAAAGGGGGGGAGGGTACATTACGGATTCGTACCTGGTGGGGGATCGTCATTTGTTTTTTAATTGTGAtt                                                                                       >  2:945418/1‑151 (MQ=255)
                                                        aaaTCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGGGGGAAAAGGGGGGGGGGGTACATTACGTATTCGTACCTGGTGGGGGGTCGTCATttttttttttattggggttt                                                                                      >  2:235863/1‑139 (MQ=255)
                                                          aTCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGGGAGAAAAAGGGGGGAGGTTAAATTACGTATTTGTACCTGGGGGGGGATCGTCATTTTTTTTTTAATGGTGAttt                                                                                      >  2:447601/1‑149 (MQ=255)
                                                              ctGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCa                                                                                    <  2:474791/147‑1 (MQ=255)
                                                              ctGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAgg                                                                                  <  1:550910/149‑1 (MQ=255)
                                                              ctGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTAACGGGGAGAAAAGGGGGGGGGGGTACATTACGTATTCGTACCTCGTGGGGGATCGTCATTTTTTTTTTTATGGTGATTTCAgg                                                                                  >  1:1092407/1‑149 (MQ=255)
                                                                 gCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGGGAGGAAAAGGGGGGAGGGTACATTACGTATTCGTACCTGGTGGGGGATCGTCATTTGTTTTTTTATGGTGAttttcgggcc                                                                               >  2:257094/1‑142 (MQ=255)
                                                                 gCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCAc                                                                             <  2:1040799/151‑1 (MQ=255)
                                                                  cAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTATCCGAGAGAAAAGGGGGGGGGGTTACATTACGTATTTTTACCTGGGGGGGGGACGTCATTTGTTTTTTAATGTGGAtttttag                                                                                  >  1:1024865/1‑141 (MQ=255)
                                                                     cacTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTTCCGGGGAGAAAAAGGGGGGAGGTTTCATTTAGTATTTCTACCTGGGGGGGGATCGTCATTTTTTTTTTTATTGGGAtttttggggcccgctg                                                                         >  1:365362/1‑138 (MQ=255)
                                                                              cgTTGTTAAGTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGGGGTTACCGGTGAGAAAAGGGGGTGATGTTACATTGCGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGAtggcgtt                                                                <  1:805175/151‑1 (MQ=255)
                                                                                  gttAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGAtgacgttgacg                                                            <  2:42726/151‑1 (MQ=255)
                                                                                  gttAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGAtgacgttgacg                                                            <  1:452729/151‑1 (MQ=255)
                                                                                  gttAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGAtgacgttgac                                                             <  2:137129/150‑1 (MQ=255)
                                                                                    tAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTAACGGAGAGGAAAGGGGGGGAGGGTACATTACGTATTTGTACCTGGTGGGGGCTCGTCATTTGGTTTTTAATGGGGATTTCAGGTCCACGCTGGATggcggtgacgtt                                                          >  2:716639/1‑151 (MQ=255)
                                                                                     aaTTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTTCCGGGGAGAAAAGGGGGGGAGGGTACATTACGGATTTGTACCCGGTGGGGGATCGTCATTTTTTTTTTCATTGGGATTTTTGGTCCCCGCTGGGTGGCGTTGTCGttt                                                         >  2:1092250/1‑151 (MQ=255)
                                                                                      aTTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTTACGGAGAGAAAAAGGGGGGGGGGTAAATTTCGGATTTGTACCCGGTGGGGGATCGACAATTTTTTTTTAATGGGGGTTTTTGGGTCCCGCTGGGtgacggttgcggttt                                                        >  1:704971/1‑141 (MQ=255)
                                                                                      aTTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGGGGGGAAAGGGGGGGGGGGTAAATTACGTATTCGTACCTGGTGGGGGATCGTCATTTTTTTTTTAATGGGGGTTTCAGGGCCACGCTGGATGACGGTGACGgttt                                                        >  2:1013739/1‑151 (MQ=255)
                                                                                      aTTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGtttt                                                        <  2:82747/151‑1 (MQ=255)
                                                                                         ttAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGtttttt                                                      <  2:164246/150‑1 (MQ=255)
                                                                                         ttAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGtttttt                                                      <  2:382805/150‑1 (MQ=255)
                                                                                          tAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGGGGGAAAAGGGGGGGAGGGTAAATTTCGTATTCGTACCCGGTGGGGGATCGGCATTTGTTTTTTAATTGTGATTTTAGGGTCACGGTGGGTGAAGGTTGCGtttttt                                                      >  1:961693/1‑149 (MQ=255)
                                                                                          tAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTg                                                     <  2:1115390/150‑1 (MQ=255)
                                                                                           aGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCCGGGAGGAAAGGGGGGGGGGGTACAATACCTATTCCTACCCGGGGGGGGCTCGGGCTTTGTTTTTTCATGGTGATTTTCGGGCCACGCTGGAGGAACTTGAAGgtttttt                                                     >  2:1077650/1‑148 (MQ=255)
                                                                                             cAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCt                                                 <  1:457775/151‑1 (MQ=255)
                                                                                              aaCGCTGATAAACCGTTTTCTGCTATGGCGGTTTACCGGGAGGAAAAGGGGGGGGGTTACAATTCGCATTTGTACCTGGGGGGGCATCGTCCTTTTTTTTTTAAAGGGGGTTTGATGTACCCACCGGGTGACGGTGACGtttttt                                                      >  2:261501/1‑145 (MQ=255)
                                                                                                cGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGGGAGAAAAAGGGGGGAGGGTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTTTTTTTTAATGGGGATTTCAGGTCCACCCTGGGTGACGTTGACGtttttt                                                      >  2:539994/1‑143 (MQ=255)
                                                                                                    gATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAg                                           <  2:824983/150‑1 (MQ=255)
                                                                                                    gATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAg                                           <  2:749103/150‑1 (MQ=255)
                                                                                                     aTAAACCGTTTTCTGCTATGGCGGTTACCGGGGAGAAAAGGGGGGGAGGGTACATTTCGTATTCGTACCTGGTGGGGGATCGTCATTTGTTTTTTAATGGTGATTTTAGGTCCACGCTGGGTGACGTTGACGTTTTTTGTGCTGATAAGtt                                         >  1:195241/1‑151 (MQ=255)
                                                                                                     aTAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGtt                                         <  1:814379/151‑1 (MQ=255)
                                                                                                     aTAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGtt                                         <  1:558887/151‑1 (MQ=255)
                                                                                                       aaaCCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGtttt                                       <  2:320681/151‑1 (MQ=255)
                                                                                                        aaCCGTTTTCTGCTATGGCGGTTACCCGGGAGAAAAAGGGGGGAGGGTACATTACGTATTCGTACCTGGTGGGGGATCGTCATTTGTTTTTCAATGGTGATTTCAGGGCCACGCTGGATGACGTTGACGGTTTTTTCGCTGATAAGtttt                                       >  2:641491/1‑150 (MQ=255)
                                                                                                            gTTTTCTGCTATGGCGGGTACCGGGGAGGAAAGGGGGGGGGGGTACATTACGTATTCGTACCTGGGGGGGGATCGTCCTTTGTTTTTTAATGGTGGTTTTAGGTCCACGCTGGGTGACGTTGACGTTTTTTTCGCTGATAAGTTTTCTg                                    >  1:1037333/1‑149 (MQ=255)
                                                                                                              tttCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCa                                 <  2:818881/150‑1 (MQ=255)
                                                                                                                   gCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTg                              <  2:757552/148‑1 (MQ=255)
                                                                                                                   gCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGa                             <  1:142071/149‑1 (MQ=255)
                                                                                                                       tGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAat                       >  1:660776/1‑151 (MQ=255)
                                                                                                                              ttACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGt                  >  2:9208/1‑150 (MQ=255)
                                                                                                                              ttACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGt                  <  2:938803/150‑1 (MQ=255)
                                                                                                                                    gagagaAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTta           <  1:244613/151‑1 (MQ=255)
                                                                                                                                    gagagaAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTta           <  2:820830/151‑1 (MQ=255)
                                                                                                                                    gagagaAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTta           <  1:801565/151‑1 (MQ=255)
                                                                                                                                        gaAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTTATATc       >  1:834322/1‑151 (MQ=255)
                                                                                                                                            aGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTTATATCTGTa   <  2:923593/151‑1 (MQ=255)
                                                                                                                                             gggggTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTTATATCTGTAc  <  1:188752/151‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
GCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTTATATCTGTAC  >  NC_000913/1478179‑1478469

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