Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,478,3240TG44.9% 122.2 / 12.6 70G33G (GGT→GGGynbBputative CDP‑diglyceride synthase
Reads supporting (aligned to +/- strand):  ref base T (8/30);  new base G (31/0);  total (40/30)
Fisher's exact test for biased strand distribution p-value = 1.65e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTTATATCTGTACATTTT  >  NC_000913/1478174‑1478474
                                                                                                                                                         |                                                                                                                                                       
ggACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGt                                                                                                                                                         <  2:813112/151‑1 (MQ=255)
     gCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAgg                                                                                                                                                     <  2:1020651/150‑1 (MQ=255)
           aaTCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACAGGGGGGAAAAAGGGGGGAGGGTACAtt                                                                                                                                              >  1:1116209/1‑151 (MQ=255)
               gTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGGGGGAAAAGGGGGGGGGGGTACa                                                                                                                                                >  1:198097/1‑145 (MQ=255)
               gTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGt                                                                                                                                          <  2:829045/151‑1 (MQ=255)
                    cGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCg                                                                                                                                     <  2:648637/151‑1 (MQ=255)
                     ggAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGGGAGAAAAGGGGGGGAGGGTACATTaaacaattgt                                                                                                                                    >  1:1063334/1‑142 (MQ=255)
                      gAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGGGGGGGTACATTACGTAtt                                                                                                                                       >  2:465643/1‑147 (MQ=255)
                       aGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTAACGGAGCGGAAAGGGGGGGAGGGTAAAtttccgatt                                                                                                                                       >  2:904550/1‑139 (MQ=255)
                           tgCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCtg                                                                                                                              <  1:515192/150‑1 (MQ=255)
                               tcttATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCtggtg                                                                                                                           <  1:479227/149‑1 (MQ=255)
                                 ttATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGGGAGAAAAGGGGGGGGGGGTACATTACGTATTTGTACCTGGTgggg                                                                                                                        >  2:739729/1‑150 (MQ=255)
                                 ttATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTAACGGGGAGGAAAAGGGGGGGGGGTACATTACGGATTTGTACCTTggggggg                                                                                                                       >  2:858003/1‑150 (MQ=255)
                                 ttATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTAACGGGGAGAAAAGGGGGGGGGGGTACATTACGTATTTGTACCTGGTggggg                                                                                                                       >  2:311359/1‑150 (MQ=255)
                                     gAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCg                                                                                                                   <  2:1031013/151‑1 (MQ=255)
                                          cTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGATGATAAACCGTTTTCTGCTATGGCGGTTTACTGAGAGAAAAGGGGGGGAGGTTTAATTAAGTATTCGTGCCTGGTgggggtaagttatt                                                                                                              >  1:1150972/1‑141 (MQ=255)
                                            gAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCAtt                                                                                                              <  2:954503/149‑1 (MQ=255)
                                            gAAATATAATGATGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACACGAGGGGAAAAGGGGGGGGGGTTCATTTCGTATTTCTTCCCTGTGGGGGATCGTCAttt                                                                                                             >  1:1145543/1‑150 (MQ=255)
                                             aaaTATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTAACGGGGGGAAAAGGGGGGGAGGGTATATTACGTTTTTGTACCTGGTGGGGGAGCGTCAttt                                                                                                             >  1:140788/1‑149 (MQ=255)
                                                tataATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTGACGGGGGGGAAAGGGGGGGAGGGTTCATTTATTTTTCGTATCCTGTGGGGGAtgctgccttttttt                                                                                                        >  1:82140/1‑138 (MQ=255)
                                                tataATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTAACCGAGAGAAAAGGGGGGGAGGGTACATTTCGTATTCGGACCCGGTGGGGGATCGTCATTTGtttt                                                                                                        >  1:461627/1‑151 (MQ=255)
                                                  taATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTAACGGGGGGGAAAAGGGGGGGGGTTTCATTTCGTTTTTGTTCCTGGGGGGGGCTCGTTATTTGtttttt                                                                                                      >  1:204264/1‑150 (MQ=255)
                                                      gCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTATCCGCGGGGAAAAGGGGGGAGGATAGATTTAGTTTTCGTTCCCGGGGGGGGATTGTCCTTTGttttt                                                                                                       >  1:159777/1‑146 (MQ=255)
                                                      gCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATg                                                                                                  <  1:854449/151‑1 (MQ=255)
                                                       cTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTTACGGAGAGAAAAGGGGGGGGGGGTACATTACTGATTTGTTCCTGGTGGGGGATCGTCATTTGtttttttat                                                                                                   >  1:921624/1‑145 (MQ=255)
                                                        tGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTGTCCGGGGGAAAAAGGGGGGGGGGTTGATTTTCGATTTTTACCTGGGGGGGCGTCTTTCTTTGttttttcattgt                                                                                                >  2:997158/1‑144 (MQ=255)
                                                        tGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTTCCGGAGAGAAAAGGGGGGGGGGGTACATTACGTATTCGTACCTGGTGGGGGATCGTCATTTTTTTTTTAATgg                                                                                                 >  2:256438/1‑150 (MQ=255)
                                                          gAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGGGGGAAAAGGGGGGGGGGGTTCATTACGTATTCGTCCCTGGGGGGGCATCGTCATTTGttttt                                                                                                       >  2:723642/1‑142 (MQ=255)
                                                          gAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGa                                                                                              <  2:381891/151‑1 (MQ=255)
                                                          gAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGa                                                                                              <  1:917062/151‑1 (MQ=255)
                                                              aaTCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACAGGAGGGAAAAAGGGGGGAGGGTTTATTACGACTTGTTCCCCGGTGGGGGATGCGCCTTTTTTTTTTTATGGttttttt                                                                                          >  1:897678/1‑145 (MQ=255)
                                                               aTCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGGGAAAAGGGGGGGGGGTTACATTACGTATTCGTACCTGGGGGGGGATCGTCATTTGTTTTTTTATGGTGAttt                                                                                           >  2:905228/1‑149 (MQ=255)
                                                               aTCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTAACGGGGAGAAAAGGGGGGGGGGGTACATTTCGTATTTGTACCTGGGGGGGGATCGGCATttttttttttattgggatt                                                                                            >  2:1072172/1‑137 (MQ=255)
                                                                tctcTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGGGAGAAAAAGGGGGGGGGGTACATTTCGTATTTGTACCTGGGGGGGGATCGTCATTTGTTTTTTTATGGTGAttt                                                                                           >  1:716607/1‑148 (MQ=255)
                                                                  tctGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTTACGGGGAGAAAAGGGGGGGAGGGTACATTAAGTATTTGTAACTGGTGGGGGATCGTCATTTGTTTTTCTATTGTGGTTTTAgg                                                                                       >  2:516306/1‑150 (MQ=255)
                                                                  tctGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTATCGGGGAGAAAAGGGGGGGAGGGTTCATTAACGATTCGTACCTGGTGGGGGGTCGTCATTTTTTTTTTTATGGGGAttttaag                                                                                       >  2:175849/1‑146 (MQ=255)
                                                                   ctGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCa                                                                                         <  1:718002/147‑1 (MQ=255)
                                                                   ctGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAg                                                                                        <  1:1020519/148‑1 (MQ=255)
                                                                    tGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTAACCGGGGGAAAAAGGGGGGAGGTTACATTACGTTTTTGTTCCCGGCGGGGGGTGGGCATTTTTTTTTTTGTGGTGGTTTTAGGGcc                                                                                    >  2:211494/1‑151 (MQ=255)
                                                                      gCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTAACGGGGAGAAAAGGGGGGGGGGGTACATTTCGTATTCGTACCTGGTGGGGGATCGTCATTTGTTTTTTAATGGTGATTTTAGGGCCAc                                                                                  >  2:572827/1‑151 (MQ=255)
                                                                         acacTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCt                                                                               <  2:1007956/151‑1 (MQ=255)
                                                                           acTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTAAAGGAGAGAAAAGGGGGTGAGGTGTTATTACCGATTCGTTACTTGGGGGGGGACGGTATTTGTTTTTTTATGTGGATTTTAAGTCCACGCTgg                                                                             >  1:589244/1‑151 (MQ=255)
                                                                                 cgcgTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGAtgac                                                                        <  1:32152/150‑1 (MQ=255)
                                                                                          aaTTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGGGGGGAAAGGGGGGGAGGGTACATTATGTATTCGTACCTGGGGGGGGGTCGTCATTTGTTTTTTAATGGGGGTTTTAGGTCCCCGCTGGATGACGTTGACGttt                                                              >  2:994549/1‑151 (MQ=255)
                                                                                             tttAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGGGGGAAAAGGGGGGGGGGGTACATTACGTATTCGTACCTGGTGGGGGATCGTCCTTTTTTTTTTTATGGTGGTTTTAGGTCCACGCTGGGTGACGTTGACGtttttt                                                           >  2:376950/1‑151 (MQ=255)
                                                                                                 gCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTgcgc                                                       <  2:1076665/151‑1 (MQ=255)
                                                                                                     cGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGAt                                                   <  2:29515/151‑1 (MQ=255)
                                                                                                      gCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGGGAGAAAAGGGGGGGGGGGTACATTACGTATTCGTACCTGGGGGGGGGTCGTCATTTTTTTTTTTATGGTGGTTTTCGGGCCCCGCTGGGTGACGTTGACGTTTTTTTTGCTg                                                     >  2:341330/1‑148 (MQ=255)
                                                                                                         gATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAg                                                <  1:256259/150‑1 (MQ=255)
                                                                                                         gATAAACCGTTTTCTGCTATGGCGGTTAACGGCGAGAAAAAGGGGTGAGGGTACATTACCTATTCGTACCCTGGGGGGGAGCGTCATTTGTTTTTTAATGGTGATTTTAGGGCCACGCTGGGTGACGTTGACGGTTTTTTTGCTGATa                                                  >  1:494000/1‑148 (MQ=255)
                                                                                                           tAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGGGGGGGTACATTACGTATTTGTACCTGGGGGGGGGTCGTCATTTTTTTTTTTATGGGGATTTTAGGTCCACGCTGGGTGACGTTTACGGTTTTTTCGGTGGTAAGttt                                             >  2:108149/1‑151 (MQ=255)
                                                                                                            aaaCCGTTTTCTGCTATGGCGGTTAACGGGGAGAAAAGGGGGGGAGGGTACATTACGTATTTGTACCTGGGGGGGGGTCGTCCTTTGTTTTTTTATGGGGATTTTAGGTCCACGCTGGGTGACGTTTGCGTTTTTTTCGCTGATAAGtttt                                            >  2:442164/1‑151 (MQ=255)
                                                                                                              aCCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCt                                          <  1:457316/151‑1 (MQ=255)
                                                                                                                 gTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGc                                       <  1:1106003/151‑1 (MQ=255)
                                                                                                                    ttCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCAt                                     <  1:1118364/150‑1 (MQ=255)
                                                                                                                        gCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGa                                  <  1:572527/149‑1 (MQ=255)
                                                                                                                                ggTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTg                        >  2:39436/1‑151 (MQ=255)
                                                                                                                                   tACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTAc                     >  2:670634/1‑151 (MQ=255)
                                                                                                                                     aCCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTAc                     <  2:459756/150‑1 (MQ=255)
                                                                                                                                     aCCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACg                    <  2:45512/151‑1 (MQ=255)
                                                                                                                                         cGGAGAGAAACGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACg                    >  1:151067/1‑149 (MQ=255)
                                                                                                                                         cGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCt                  <  2:788587/151‑1 (MQ=255)
                                                                                                                                          ggAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCtt                 >  1:1119672/1‑151 (MQ=255)
                                                                                                                                            agagaAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTtat               <  2:198272/151‑1 (MQ=255)
                                                                                                                                             gagaAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTtata              <  2:250301/151‑1 (MQ=255)
                                                                                                                                                aaaaGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTTATATCt           <  2:486508/151‑1 (MQ=255)
                                                                                                                                                aaaaGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTTATATCt           >  1:889237/1‑151 (MQ=255)
                                                                                                                                                 aaaGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTTATATCTg          >  1:376834/1‑151 (MQ=255)
                                                                                                                                                     ggggTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTTATATCTGTACa      <  2:758694/151‑1 (MQ=255)
                                                                                                                                                      gggTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTTATATCTGTACAt     <  1:1017672/151‑1 (MQ=255)
                                                                                                                                                       ggTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTTATATCTGTACAtt    <  2:1086377/151‑1 (MQ=255)
                                                                                                                                                         tGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTTATATCTGTACAtttt  >  2:406247/1‑151 (MQ=255)
                                                                                                                                                         |                                                                                                                                                       
GGACTGCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTTATATCTGTACATTTT  >  NC_000913/1478174‑1478474

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