Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,638,9110AG45.9% 16.9 / 82.0 61M424V (ATG→GTG) pitAphosphate transporter, low‑affinity; tellurite importer
Reads supporting (aligned to +/- strand):  ref base A (15/18);  new base G (10/18);  total (25/36)
Fisher's exact test for biased strand distribution p-value = 6.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.74e-01

GACATGCTTAGCACCATCGAGTATGCACCGGTGTGGATCATCATGGCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACGCTCAGGGGATGTCTGCCCAGATGACGGCGGCAGTGTCTATCGGCCTGGCGAGTTATACCGGGATGCCGGTTTCCACTACTCACGTACTCTCCTCTTCTGTCGCGGGGACGATGGTGGTAGATGGTGGCGGCTTACAGCGTAAAACCGTGACCAG  >  NC_000913/3638761‑3639056
                                                                                                                                                      |                                                                                                                                                 
gACATGCTTAGCACCATCGAGTATGCACCGGTGTGGATCATCATGGCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACGCTCAGGGGa                                                                                                                                                   <  2:4080/151‑1 (MQ=255)
       ttAGCACCATCGAGTATGCACCGGTGTGGATCATCATGGCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACGCTCAGGGGGTGTCTGc                                                                                                                                            >  1:446051/1‑151 (MQ=255)
       ttAGCACCATCGAGTATGCACCGGTGTGGATCATCATGGCGGTCGCGCCGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGGGTGGCAACGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACGCTCAGGGGGTGTCTGc                                                                                                                                            >  2:193346/1‑151 (MQ=255)
        tAGCACCATCGAGTATGCACCGGTGTGGATCATCATGGCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACGCTCAGGGGGTGTCTg                                                                                                                                             <  2:235961/149‑1 (MQ=255)
        tAGCACCATCGAGTATGCACCGGTGTGGATCATCATGGCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACGCTCAGGGGGTGTCTg                                                                                                                                             <  2:255854/149‑1 (MQ=255)
          gCACCATCGAGTATGCACCGGTGTGGATCATCATGGCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACGCTCAGGGGATGTCTGccc                                                                                                                                          >  2:336118/1‑150 (MQ=255)
                 cGAGTATGCACCGGTGTGGATCATCATGGCGGCCGCGCTGGCGTTAGGGATCGGTACGATGGTTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACGCTCAGGGGGTGTCTGCCCAGATGAc                                                                                                                                   >  1:386949/1‑150 (MQ=255)
                  gAGTATGCACCGGTGTGGATCATCATGGCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACGCTCAGGGGATGTCTGCCCAGATGAcgg                                                                                                                                 >  1:9313/1‑151 (MQ=255)
                  gAGTATGCACCGGTGTGGATCATCATGGCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACGCTCAGGGGATGTCTGCCCAGATGAcgg                                                                                                                                 <  2:451748/151‑1 (MQ=255)
                     tATGCACCGGTGTGGATCATCATGGCGGTCGCGTCCGCGTTAAGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATTGGTGAGAAAATCGGTAAGAAAGGCATGACCCACGCTCAGGGGATGTCTGCCCAGATGAcggcgg                                                                                                                              >  2:553493/1‑151 (MQ=38)
                        gCACCGGTGTGGATCATCATGGCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACGCTCAGGGGGTGTCTGCCCAGATGACGGCGGCAg                                                                                                                           <  1:67609/151‑1 (MQ=255)
                                gtgGATCATCATGGCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGGGCCGTGTGGCAACGACTATCGGGGAGAAAATCGGTAAGAAAGGCATGACCTACGCTCAGGGGGGGTCTGCCCAGATGACGGCGGCAGTGTCTATc                                                                                                                   >  2:204196/1‑151 (MQ=255)
                                          atGGCGGTCGCGCTGGCGTTAGGTATCGGTACTATGATTGGCTGGCGCCGTGTGGCAACGACCATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACGCTCAGGGGATGTCTGCCCAGATGACGGCGGCAGTGTCTATCGGCCTGGCg                                                                                                          >  2:360236/1‑150 (MQ=255)
                                            ggCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACGCTCAGGGGATGTCTGCCCAGATGACGGCGGCAGTGTCTATCGGCCTGGCGAg                                                                                                        <  1:301795/150‑1 (MQ=255)
                                               ggTCGCGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACGCTCAGGGGGTGTCTGCCCAGATGACGGCGGCAGTGTCTATCGGCCTGGCGAGTtat                                                                                                    >  2:285904/1‑151 (MQ=255)
                                                  cgcgCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGGCCTACGCTCAGGGGGTGTCTGCCCAGATGACGGCGGCAGTGTCTATCGGCCTGGCGAGTTATAcc                                                                                                 <  2:366149/151‑1 (MQ=255)
                                                  cgcgCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACGCTCAGGGGGTGTCTGCCCAGATGACGGCGGCAGTGTCTATCGGCCTGGCGAGTTATAcc                                                                                                 <  1:309287/151‑1 (MQ=255)
                                                   gcgcTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACGCTCAGGGGATGTCTGCCCAGATGACGGCGGCAGTGTCTATCGGCCTGGCGAGTTATAcc                                                                                                 >  2:570778/1‑150 (MQ=255)
                                                   gcgcTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACGCTCAGGGGATGTCTGCCCAGATGACGGCGGCAGTGTCTATCGGCCTGGCGAGTTATACCg                                                                                                <  1:336125/151‑1 (MQ=255)
                                                    cgcTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACGCTCAGGGGATGTCTGCCCAGATGACGGCGGCAGTGTCTATCGGCCTGGCGAGTTATACCgg                                                                                               <  1:256657/151‑1 (MQ=255)
                                                        ggCGGTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACGCTCAGGGGATGTCTGCCCAGATGACGGCGGCAGTGTCTATCGGCCTGGCGAGTTATACCGGGATg                                                                                           >  1:228007/1‑151 (MQ=255)
                                                           gTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACGCTCAGGGGGTGTCTGCCCAGATGACGGCGGCAGTGTCTATCGGCCTGGCGAGTTATACCGGGATGCCg                                                                                        <  1:377163/151‑1 (MQ=255)
                                                           gTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACGCTCAGGGGGTGTCTGCCCAGATGACGGCGGCAGTGTCTATCGGCCTGGCGAGTTATACCGGGATGCCg                                                                                        <  2:592989/151‑1 (MQ=255)
                                                              aGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACGCTCAGGGGATGTCTGCCCAGATGACGGCGGCAGTGTCTATCGGCCTGGCGAGTTATACCGGGATGCCGGt                                                                                      <  1:445355/150‑1 (MQ=255)
                                                              aGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACGCTCAGGGGATGTCTGCCCAGATGACGGCGGCAGTGTCTATCGGCCTGGCGAGTTATACCGGGATGCCGGt                                                                                      <  1:248741/150‑1 (MQ=255)
                                                                    cGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACGCTCAGGGGATGTCTGCCCAGATGACGGCGGCAGTGTCTATCGGCCTGGCGAGTTATACCGGGATGCCGGTTTCCact                                                                               >  1:28403/1‑151 (MQ=255)
                                                                        aCGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACGCTCAGGGGATGTCTGCCCAGATGACGGCGGCAGTGTCTATCGGCCTGGCGAGTTATACCGGGATGCCGGTTTCCACTACTc                                                                           <  2:507687/151‑1 (MQ=255)
                                                                              atTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACGCTCAGGGGATGTCTGCCCAGATGACGGCGGCAGTGTCTATCGGCCTGGCGAGTTATACCGGGATGCCGGTTTCCACTACTCACGTAc                                                                     <  2:435464/151‑1 (MQ=255)
                                                                                       cgcCGTGTGGCAACGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACGCTCAGGGGGTGTCTGCCCAGATGACGGCGGCAGTGTCTATCGGCCTGGCGAGTTATACCGGGATGCCGGTTTCCACTACTCACGTACTCTCCtctt                                                            >  1:34418/1‑151 (MQ=255)
                                                                                       cgcCGTGTGGCAACGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACGCTCAGGGGATGTCTGCCCAGATGACGGCGGCAGTGTCTATCGGCCTGGCGAGTTATACCGGGATGCCGGTTTCCACTACTCACGTACTctcctc                                                              <  2:372383/149‑1 (MQ=255)
                                                                                       cgcCGTGTGGCAACGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACGCTCAGGGGATGTCTGCCCAGATGACGGCGGCAGTGTCTATCGGCCTGGCGAGTTATACCGGGATGCCGGTTTCCACTACTCACGTACTCTCctct                                                             >  2:144745/1‑150 (MQ=255)
                                                                                       cgcCGTGTGGCAACGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACGCTCAGGGGATGTCTGCCCAGATGACGGCGGCAGTGTCTATCGGCCTGGCGAGTTATACCGGGATGCCGGTTTCCACTACTCACGTACTCTCCtctt                                                            <  1:297759/151‑1 (MQ=255)
                                                                                         ccGTGTGGCAACGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACGCTCAGGGGATGTCTGCCCAGATGACGGCGGCAGTGTCTATCGGCCTGGCGAGTTATACCGGGATGCCGGTTTCCACTACTCACGTACTCTCCtcttct                                                          <  2:184241/151‑1 (MQ=255)
                                                                                         ccGTGTGGCAACGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACGCTCAGGGGATGTCTGCCCAGATGACGGCGGCAGTGTCTATCGGCCTGGCGAGTTATACCGGGATGCCGGTTTCCACTACTCACGTACTCTCCtcttct                                                          >  2:527001/1‑151 (MQ=255)
                                                                                           gtgtGGCAACGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACGCTCAGGGGATGTCTGCCCAGATGACGGCGGCAGTGTCTATCGGCCTGGCGAGTTATACCGGGATGCCGGTTTCCACTACTCACGTACTCTCCTCTTCTGt                                                        >  2:163492/1‑151 (MQ=255)
                                                                                              tgGCAACGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACGCTCAGGGGATGTCTGCCCAGATGACGGCGGCAGTGTCTATCGGCCTGGCGAGTTATACCGGGATGCCGGTTTCCACTACTCACGTACTCTCCTCTTCTGTcgc                                                     >  2:171103/1‑151 (MQ=255)
                                                                                                 cAACGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACGCTCAGGGGATGTCTGCCCAGATGACGGCGGCAGTGTCTATCGGCCTGGCGAGTTATACCGGGATGCCGGTTTCCACTACTCACGTACTCTCCTCTTCTGTCGCggg                                                  >  1:561120/1‑151 (MQ=255)
                                                                                                   aCGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACGCTCAGGGGGTGTCTGCCCAGATGACGGCGGCAGTGTCTATCGGCCTGGCGAGTTATACCGGGATGCCGGTTTCCACTACTCACGTACTCTCCTCTTCTGTCGCGGGGa                                                <  1:156234/151‑1 (MQ=255)
                                                                                                    cGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACGCTCAGGGGGTGTCTGCCCAGATGACGGCGGCAGTGTCTATCGGCCTGGCGAGTTATACCGGGATGCCGGTTTCCACTACTCACGTACTCTCCTCTTCTGTCGCGGGGAc                                               <  2:285506/151‑1 (MQ=255)
                                                                                                    cGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACGCTCAGGGGATGTCTGCCCAGATGACGGCGGCAGTGTCTATCGGCCTGGCGAGTTATACCGGGATGCCGGTTTCCACTACTCACGTACTCTCCTCTTCTGTCGCGGGGAc                                               >  2:617398/1‑151 (MQ=255)
                                                                                                    cGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACGCTCAGGGGATGTCTGCCCAGATGACGGCGGCAGTGTCTATCGGCCTGGCGAGTTATACCGGGATGCCGGTTTCCACTACTCACGTACTCTCCTCTTCTGTCGCGGGGAc                                               <  2:429407/151‑1 (MQ=255)
                                                                                                     gACTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACGCTCAGGGGATGTCTGCCCAGATGACGGCGGCAGTGTCTATCGGCCTGGCGAGTTATACCGGGATGCCGGTTTCCACTACTCACGTACTCTCCTCTTCTGTCGCGGGGACg                                              <  1:307239/151‑1 (MQ=255)
                                                                                                      aCTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACGCTCAGGGGATGTCTGCCCAGATGACGGCGGCAGTGTCTATCGGCCTGGCGAGTTATACCGGGATGCCGGTTTCCACTACTCACGTACTCTCCTCTTCTGTCGCGGGGACGa                                             <  2:573257/151‑1 (MQ=255)
                                                                                                            ggTGAGAAAATCGGTAAGAAAGGCATGACCTACGCTCAGGGGATGTCTGCCCAGATGACGGCGGCAGTGTCTATCGGCCTGGCGAGTTATACCGGGATGCCGGTTTCCACTACTCACGTACTCTCCTCTTCTGTCGCGGGGACGAtggtgg                                       >  1:124391/1‑151 (MQ=255)
                                                                                                               gagaAAATCGGTAAGAAAGGCATGACCTACGCTCAGGGGATGTCTGCCCAGATGACGGCGGCAGTGTCTATCGGCCTGGCGAGTTATACCGGGATGCCGGTTTCCACTACTCACGTACTCTCCTCTTCTGTCGCGGGGACGATGGTGGTa                                     <  1:17608/150‑1 (MQ=255)
                                                                                                                     aTCGGTAAGAAAGGCATGACCTACGCTCAGGGGGTGTCTGCCCAGATGACGGCGGCAGTGTCTATCGGCCTGGCGAGTTATACCGGGATGCCGGTTTCCACTACTCACGTACTCTCCTCTTCTGTCGCGGGGACGATGGTGGTAGAtggtg                              <  1:322384/151‑1 (MQ=255)
                                                                                                                       cGGTAAGAAAGGCATGACCTACGCTCAGGGGGTGTCTGCCCAGATGACGGCGGCAGTGTCTATCGGCCTGGCGAGTTATACCGGGATGCCGGTTTCCACTACTCACGTACTCTCCTCTTCTGTCGCGGGGACGATGGTGGTAGAtggtgg                             <  2:118038/150‑1 (MQ=255)
                                                                                                                           aaGAAAGGCATGACCTACGCTCAGGGGGTGTCTGCCCAGATGACGGCGGCAGTGTCTATCGGCCTGGCGAGTTATACCGGGATGCCGGTTTCCACTACTCACGTACTCTCCTCTTCTGTCGCGGGGACGATGGTGGTAGATGGTggcggc                         <  2:582932/150‑1 (MQ=255)
                                                                                                                            aGAAAGGCATGACCTACGCTCAGGGGGTGTCTGCCCAGATGACGGCGGCAGTGTCTATCGGCCTGGCGAGTTATACCGGGATGCCGGTTTCCACTACTCACGTACTCTCCTCTTCTGTCGCGGGGACGATGGTGGTAGATGGTGGCGGCtt                       <  1:338678/151‑1 (MQ=255)
                                                                                                                            aGAAAGGCAAGACCTACGCTCAGGGGGTGTCTGCCCAGATGACGGCGGCAGTGTCTATCGGCCTGGCGAGTTATACCGGGATGCCGGTTTCCACTACTCACGTACTCTCCTCTTCTGTCGCGGGGACGATGGTGGTAGATGGTGGCGGCtt                       <  2:114268/151‑1 (MQ=255)
                                                                                                                             gAAAGGCATGACCTACGCTCAGGGGGTGTCTGCCCAGATGACGGCGGCAGTGTCTATCGGCCTGGCGAGTTATACCGGGATGCCGGTTTCCACTACTCACGTACTCTCCTCTTCTGTCGCGGGGACGATGGTGGTAGATGGTGGCGGCtt                       >  1:449667/1‑150 (MQ=255)
                                                                                                                                aGGCATGACCTACGCTCAGGGGGTGTCTGCCCAGATGACGGCGGCAGTGTCTATCGGCCTGGCGAGTTATACCGGGATGCCGGTTTCCACTACTCACGTACTCTCCTCTTCTGTCGCGGGGACGATGGTGGTAGATGGTGGCGGCTTAc                     >  2:64326/1‑149 (MQ=255)
                                                                                                                                aGGCATGACCTACGCTCAGGGGGTGTCTGCCCAGATGACGGCGGCAGTGTCTATCGGCCTGGCGAGTTATACCGGGATGCCGGTTTCCACTACTCACGTACTCTCCTCTTCTGTCGCGGGGACGATGGTGGTAGATGGTGGCGGCTTAc                     >  2:313759/1‑149 (MQ=255)
                                                                                                                                 ggCATGACCTACGCTCAGGGGATGTCTGCCCAGATGACGGCGGCAGTGTCTATCGGCCTGGCGAGTTATACCGGGATGCCGGTTTCCACTACTCACGTACTCTCCTCTTCTGTCGCGGGGACGATGGTGGTAGATGGTGGCGGCTTACAGc                  <  1:376521/151‑1 (MQ=255)
                                                                                                                                      gACCTACGCTCAGGGGGTGTCTGCCCAGATGACGGCGGCAGTGTCTATCGGCCTGGCGAGTTATACCGGGATGCCGGTTTCCGCTACTCACGTACTCTCCTCTTCTGTCGCGGGGACGATGGTGGTAGATGGTGGCGGCTTACAGCGTaaa             <  1:177272/151‑1 (MQ=255)
                                                                                                                                      gACCTACGCTCAGGGGATGTCTGCCCAGATGACGGCGGCAGTGTCTATCGGCCTGGCGAGTTATACCGGGATGCCGGTTTCCACTACTCACGTACTCTCCTCTTCTGTCGCGGGGACGATGGTGGTAGATGGTGGCGGCTTACAGCGTaa              <  2:373820/150‑1 (MQ=255)
                                                                                                                                       aCCTACGCTCAGGGGGTGTCTGCCCAGATGACGGCGGCAGTGTCTATCGGCCTGGCGAGTTATACCGGGATGCCGGTTTCCACTACTCACGTACTCTCCTCTTCTGTCGCGGGGACGATGGTGGTAGATGGTGGCGGCTTACAGCGTaaaa            <  2:356920/151‑1 (MQ=255)
                                                                                                                                              cTCAGGGGATGTCTGCCCAGATGACGGCGGCAGTGTCTATCGGCCTGGCGAGTTATACCGGGATGCCGGTTTCCACTACTCACGTACTCTCCTCTTCTGTCGCGGGGACGATGGTGGTAGATGGTGGCGGCTTACAGCGTAAAACCGTGAc     >  1:169278/1‑151 (MQ=255)
                                                                                                                                                 aGGGGGTGTCTGCCCAGATGACGGCGGCAGTGTCTATCGGCCTGGCGAGTTATACCGGGATGCCGGTTTCCACTACTCACGTACTCTCCTCTTCTGTCGCGGGGACGATGGTGGTAGATGGTGGCGGCTTACAGCGTAAAACCGTGAcc    >  2:81158/1‑149 (MQ=255)
                                                                                                                                                 aGGGGGTGTCTGCCCAGATGACGGCGGCAGTGTCTATCGGCCTGGCGAGTTATACCGGGATGCCGGTTTCCACTACTCACGTACTCTCCTCTTCTGTCGCGGGGACGATGGTGGTAGATGGTGGCGGCTTACAGCGTAAAACCGTGACCAg  <  1:326901/151‑1 (MQ=255)
                                                                                                                                                 aGGGGGTGTCTGCCCAGATGACGGCGGCAGTGTCTATCGGCCTGGCGAGTTATACCGGGATGCCGGTTTCCACTACTCACGTACTCTCCTCTTCTGTCGCGGGGACGATGGTGGTAGATGGTGGCGGCTTACAGCGTAAAACCGTGACCAg  <  1:530327/151‑1 (MQ=255)
                                                                                                                                                 aGGGGATGTCTGCCCAGATGACGGCGGCAGTGTCTATCGGCCTGGCGAGTTATACCGGGATGCCGGTTTCCACTACTCACGTACTCTCCTCTTCTGTCGCGGGGACGATGGTGGTAGATGGTGGCGGCTTACAGCGTAAAACCGTGACCAg  <  1:433498/151‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                 
GACATGCTTAGCACCATCGAGTATGCACCGGTGTGGATCATCATGGCGGTCGCGCTGGCGTTAGGTATCGGTACGATGATTGGCTGGCGCCGTGTGGCAACGACTATCGGTGAGAAAATCGGTAAGAAAGGCATGACCTACGCTCAGGGGATGTCTGCCCAGATGACGGCGGCAGTGTCTATCGGCCTGGCGAGTTATACCGGGATGCCGGTTTCCACTACTCACGTACTCTCCTCTTCTGTCGCGGGGACGATGGTGGTAGATGGTGGCGGCTTACAGCGTAAAACCGTGACCAG  >  NC_000913/3638761‑3639056

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