Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,549,6100GT22.7% 98.4 / 13.3 66W737L (TGG→TTG) fdnGformate dehydrogenase‑N, alpha subunit, nitrate‑inducible
Reads supporting (aligned to +/- strand):  ref base G (32/19);  new base T (1/14);  total (33/33)
Fisher's exact test for biased strand distribution p-value = 2.09e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CGCAATCTGACGAAGTGGCTAAAGAGAACAACGGCTACGCGCTGGAAGATCTCTATGACGCTAATGGCGTGCTTATTGCGAAGAAAGGTCAGTTGCTGAGTAGCTTTGCGCATCTGCGTGATGACGGTACAACCGCATCTTCTTGCTGGATCTACACCGGTAGCTGGACAGAGCAGGGCAACCAGATGGCTAACCGCGATAACTCCGACCCGTCCGGTCTGGGGAATACGCTGGGATGGGCCTGGGCGTGGCCGCTCAACCGTCGCGTGCTGTACAACCGTGCTTCGGCGGATAT  >  NC_000913/1549463‑1549757
                                                                                                                                                   |                                                                                                                                                   
cGCAATCTGACGAAGTGGCTAAAGAGAACAACGGCTACGCGCTGGAAGATCTCTATGACGCTAATGGCGTGCTTATTGCGAAGAAAGGTCAGTTGCTGAGTAGCTTTGCGCATCGGCGTGATGCCGGTACAACCGCATCTTCTTGCTgg                                                                                                                                                    >  1:523120/1‑149 (MQ=255)
 gCAATCTCACGAAGTGTCTAAAGAGAACAACGTCTACGGGCTGGAAGATCTCTATGTAGCTAATGTCGTGCTTATTGCGAATAAAGTTCAGTTTATGAGTATCTTTGCTCATCTGCGTGATGACGGTACAAACTCATCTTATTGCTGGa                                                                                                                                                   <  2:154633/149‑1 (MQ=255)
      cTGACGAAGTGGCTAAAGAGAACAACGGCTACGCGCTGGAAGATCTCTATGACGCTAATGGCGTGCTTATTGCGAAGAAAGGTCAGTTGCTGAGTAGCTTTGCGCATCTGCGTGCTGACGGTACAACCGCATCTTCTTGCTGGATCTac                                                                                                                                              >  2:972732/1‑149 (MQ=255)
        gACGAAGTGGCTAAAGAGAACAACGGCTACGCGCTGGAAGATCTCTATGACGCTAATGGCGTGCTTATTGCGAAGAAAGGTCAGTTGCTGAGTAGCTTTGCGCATCTGCGTGATGACGGTACAACCGCATCTTCTTGCTGGATCTacac                                                                                                                                            >  2:280482/1‑149 (MQ=255)
          cGAAGTGGCTAAAGAGAACAACGGCTACGCGCTGGAAGATCTCTATGACGCTAATGGCGTGCTTATTCCGAAGAAACGTCAGTTGATCAGTAGCTTTGCGCATCTGCGTGATGACGGTACAACCGCATCTTAATGCTGGATCTACACCg                                                                                                                                          >  2:424603/1‑149 (MQ=255)
            aaGTGGCTAAAGAGAACAACGGCTACGCGCTGGAAGATCTCTATGACGCTAATGGCGTGCTTATTGCGAAGAAAGGTCAGTTGCTGAGTAGCTTTGCGCATCTGCGTGATGACGGTACAACCGCATCTTCTTGCTGGATCTACACCGGt                                                                                                                                        <  1:373638/149‑1 (MQ=255)
                  ctcacgcgaACACCGGCTCCGCGCTGGAAGAGCTCTATGACGCTAAGGCCGTGCTTCTTGCGAAGAACGGCCAGTTGCCGAGTAGCTTTGCGCATCTGCGTTATTTTTTTACATATTCTGCGACTAGCTTGATCTACACCGGTAGCTgg                                                                                                                                  <  2:75402/142‑1 (MQ=255)
                   tAAAGAGAACAACGGCTACGCGCTGGAAGATCTCTATGACGCTAATGGCGTGCTTATTGCGAAGAAAGGTCAGTTGCTGAGTAGCTTTGCGCATCTGCGTGATGACGGTACAACCGCATCTTCTTGCTGGTTCTACACCGGTAGCTGGa                                                                                                                                 >  2:850658/1‑149 (MQ=255)
                    aaaGAGAACAACGGCTACGCGCTGGAAGATCTCTATGACGCTAATGGCGTGCTTATTGCGAAGAAAGGTCAGTTGCTGAGTAGCTTTGCGCATCTGCGTTATGTTTTTAATTATTCTTCTTCTTGCTTGATCTACACCGGTAGCTGGAc                                                                                                                                <  2:654023/149‑1 (MQ=255)
                    aaaGAGAACAACGGCTACGCGCTGGAAGATCTCTATGACGCTAATGGCGTGCTTATTGCGAAGAAAGGTCAGTTGCTGAGTAGCTTTGCGCATCTGCGTGATGACGGTACAACCGCATCTTCTTGCTGGATCTACACCGGTAGCTGGAc                                                                                                                                >  2:321167/1‑149 (MQ=255)
                       gagaACAACGGCTAGGCGCTGGAAGATCTCTATGTCGCTAATGGCGTGCTTATTGCGAAGAATGTTCAGTTGCTGCGTAGCTTTGCGCATCTGCGTGATGACGGTAGAACCGCATCTTCTTGCTGGATCTACACCGGTAGCTGGACaga                                                                                                                             <  1:126984/149‑1 (MQ=255)
                       gagaACAACGGCTACGCGCTGGAAGATCTCTATGACGCTAATGGCGTGCTTATTGCGAAGAAAGGTCAGTTGCTGAGTAGCTTAGCGCATCTGCGTGATGACGGTACAAACGCATCTTATTGCTGGATCTACACCGGTAGATGGACaga                                                                                                                             >  1:750902/1‑149 (MQ=255)
                        agaACAACGGCTACGCGCTGGAAGATCTCTATGACGCTAATGGCGTGCTTATTGCGAAGAAAGGTCAGTTGCTGAGTAGCTTTGCGCATCTGCGTGATGACGGTACAACCGCATCTTCTTGCTGGATCTACACCGGTAGCTGGACagag                                                                                                                            >  1:1281142/1‑149 (MQ=255)
                          aacaacGGCTACGCGCTGGAAGATCTCTATGACGCTAATGGCGTGCTTATTGCGAAGAAAGGTCAGTTGCTGAGTAGCTTTGCGCATCTGCGTGATGACGGTACAACCGCATCTTCTTGCTGGATCTACACCGGTAGCTGGACAGAGcc                                                                                                                          >  1:678699/1‑148 (MQ=255)
                           acaacGGCTACGCGCAGGAAGATCTCTATGACGCTAATGGCGTGCTTATTGCGAAGAAAGGTCAGTTGCTGAGTAGCTTTGCGCATCTGAGTTATTTTTTTAATTATTCTTCTTCTTGCTTGATCTACACCGGTAGCTGGACAGAGCAg                                                                                                                         <  2:367736/149‑1 (MQ=255)
                            ccacGGCTACGCGCTGGAAGATCGCTTTGTCGCTAATGGCGGGCTGACTGAGAAGAAAGGTCAGGTGCTGAGTAGCTTTGCGCATCTGCGTTTTTTTTTTAATTATTCTTCTTCTTGCGTGATCTACAGCGGTAGCTGGACAGAGCAgg                                                                                                                        <  2:998796/147‑1 (MQ=255)
                            caacGGCTACGCGCTGGAAGATCTCTATGACGCTAATGGCGTGATTATTGCGAAGAAAGGTCAGTTGCTGAGTAGCTTTGCGCATCTGCGTGATGACGGTACAACCGCATCTACTTGCTGGATCTACACCCGTAGCTGGACATAGCAgg                                                                                                                        >  2:309229/1‑149 (MQ=255)
                             aacGGCTCCGCGCTGGAAGATCTCTATGACGCAAATGGCGCGCTTAATGCGAAGAAAGGTCAGTGGCAGAGTAGCTTTGCGCATCTGCGTTTTTTTTTTAATTATTCTTCTTCTTGCGTGATCTACACCGGTAGCTGGACAGAGCAggg                                                                                                                       <  2:864429/149‑1 (MQ=255)
                                  cTACGCGCTGGAAGCTCTCTATGCCGCTAATGGCGTGCTTATTGCGAAGAAAGGTCAGTTGCTGAGTAGCTTTGCGCATCTGCGTGATGACGGTACAACCGCATCTTCTTGCTGGATCTACACCGGTAGCTGGACAGAGCAGGGCAAcc                                                                                                                  <  2:1045119/149‑1 (MQ=255)
                                               gATCTCTATGACGCTAATGGCGTGCTTATTGCGAAGAAAGGTCAGTTGCTGAGTAGCTTTGCGCATCTGCGTGATGACGGTACAACCGCATCTTCTTGCTGGATCTACACCGGTAGCTGGACAGCGCAGGGCAACCAGATGGCTAACcg                                                                                                     >  1:1230903/1‑149 (MQ=255)
                                                 tctcTATGACGCTAATGGCGTGCTTATTGCGAAGAAAGGTCAGTTGCTGAGTAGCTTTGCGCATCTGCGTGATGACGGTACAACCGCATCTTCTTGCTGGATCTACACCGGTAGCTGGACAGAGCAGGGCAACCAGATGGCTAACcgcg                                                                                                   >  2:814271/1‑149 (MQ=255)
                                                  gtctATGACTCTAATGGCGGGCTTTTTGCGAAGAATGTTCAGTTGCTGAGTAGCTTTGCGCATCTGCGTGATTTCTGTACAACTGCATCTTCTTGCTGGATCTACACCGGTAGCTGGACAGAGCAGGGCAACCAGATGGCTAACCGCGa                                                                                                  <  1:716718/148‑1 (MQ=255)
                                                  ctctATGACGCTAATGGTGTGCTTATTGCGAAGAAAGGTCAGTTGCTGAGTAGCTTTGCGCATCTGCGTGATGACGGTACAACCGCATCTTCTTGCTGGATCTACACCGGTAGCTGGACAGAGCAGGGCAACCAGATGGCTAACCGCGa                                                                                                  <  1:690083/149‑1 (MQ=255)
                                                  ctctATGACGCTAATGGCGTGCTTAATGCGAAGAAAGGTCAGTTGCTGAGTAGCTTTGCGCATCTGCGTGATGTAGTTAATTAAGCATCTTCTTGCTTGATCTACACCGGTAGCTGGACAGAGCAGGGCAACCAGATGGCTAACCGCGa                                                                                                  <  2:382013/149‑1 (MQ=255)
                                                      aTGACCCTAATGGCGTGCTCATTGCGAAGAAAGGTCAGTTGCTGAGTAGCTTTGCGGATCTGATTTTTTTTTTTAATTATTCTTCTTCTTGCTTGATCTACACCGGTAGCTGGACAGAGCAGGGCAACCAGATGGCTAACCGCGATAAc                                                                                              <  2:290174/149‑1 (MQ=255)
                                                       tGACGCTAATGGCGTGCTTATTGAGAAGAAAGGTCAGTTGCTGAGTAGCTTTGCGCATCTGCGTGATGACGGTACAACCGCATCTTCTTGCTGGATCTACAACGGTAGCTGGACAGAGCAGGGCAACCAGATGGCTAACCGCGATAACt                                                                                             >  2:658832/1‑149 (MQ=255)
                                                        gACGCTAATGGCGTGCTTATTGCGAAGAAAGGTCAGTTGCTGAGTAGCTTTGCGCATCTGCGTGATGACGGTACAACCGCATCTTCTTGCTGGATCTACACCGGTAGCTg                                                                                                                                   >  1:1031851/1‑110 (MQ=255)
                                                        gACGCTAATGGCGTGCTTATTGCGAAGAAAGGTCAGTTGCTGAGTAGCTTTGCGCATCTGCGTGATGACGGTACAACCGCATCTTCTTGCTGGATCTACACCGGTAGCTg                                                                                                                                   <  2:1031851/110‑1 (MQ=255)
                                                         ccGCTAATGGCGTGCTTATTGCGAAGAAAGGGCAGTTGCTGAGTAGCTTTGCGCATCTGCGTGATGATTTTAATTATTCTTCTTCTTGCTTGATCTACACCGGTAGCTGGACAGAGCAGGGCAACCAGATGGCTAACCGCGATAACTcc                                                                                           <  2:354141/148‑1 (MQ=255)
                                                           gCTAATGGCGTGCTTATTGCGAAGAAAGGTCAGTTGCTGAGTAGCTTTGCGCATCTGCGTTTTTACTGTAAAAAATCTTCTTCTTGCTTGATCTACACCGGTAGCTGGACAGAGCAGGGCAACCAGATGGCTAACCGCGATAACTCCGa                                                                                         <  2:1152810/149‑1 (MQ=255)
                                                           gCTAATGGCGTGCTTATTGCGAAGAAAGGTCAGTTGCTGAGTAGCTTTGCGCATCTGCGTGATGACGGTACAACCGCATCTTCTTGCTTGATCTACACCGGTAGCTGGACAGAGCAGGGCAACCAGATGGCTAACCGCGATAACTCCGa                                                                                         >  1:1068289/1‑149 (MQ=255)
                                                           gCTAATGGCGTGCTTATTGCGAAGAAAGGTCAGTTGCTGAGTAGCTTTGCGCATCTGCGTGATGACGGTACAACCGCATCTTCTTGCTGGATCTACACCGGTAGCTGGACAGAGCAGGGCAACCAGATGGCTAACCGCGATAACTCCGa                                                                                         <  1:926757/149‑1 (MQ=255)
                                                            cTAATGGCGTGCTTATTGCGAAGAAAGGTCAGTTGCTGAGTAGCTTTGCGCATCTGCGTGATGACGGTACAACCGCATCTTCTTGCTGGATCTACACCGGTAGCTGGACAGAGCAGGGCAACCAGATGGCTAACCGCGATAACTCCGAc                                                                                        >  2:1162984/1‑149 (MQ=255)
                                                              aaTGCCGTGCTTATAGCGAAGAAAGGTCAGTTGCTGAGGAGCTTTGCGCATCTGCTTTATTTATTTAATTATTCTGCTTCTTGCTTGATCTACACCGGTAGCTGGACAGAGCAGGGCAACCAGATGGCTAACCGCGATAACTCCGAccc                                                                                      <  2:824965/149‑1 (MQ=255)
                                                                 ggCGTGCTTATTGCGAAGAAAGGTCAGTTGCTGAGTAGCTTTGCGCATCTGCGTGATGACGGTACAACCGCATCTTCTTGCTGGATCTACACCGGTAGCTGGACAGAGCAGGGCAACCAGATGGCTAACCGCGATAACTCCGACCCGTc                                                                                   <  1:500722/149‑1 (MQ=255)
                                                                         ttttGCGAAGAAAGGTCAGTTGCGGAGTATCTTTGCGGATCTGCGTGTTGACGGTACAACCGCATCTGCTTGGTGGATCTACACCGGTAGCTGGACAGAGCAGGGCAACCAGATGGCTAACCGCGATAACTCCGACCCGTCCGGTCTgg                                                                           <  1:1048253/147‑1 (MQ=255)
                                                                          tttGCGAAGAAAGGTCAGTTGCTGAGTAGGTTTGCGCATCTGCGTGAAGACGGTACAACCGCATCTTCTTGCTGGATCTACACCGGTAGCTGGACAGAGCAGGGCAACCACATGGCTAACCGCGATAACTCCGACCCGTCCGGTCTggg                                                                          <  1:544083/148‑1 (MQ=255)
                                                                                atgaaAGGTCAGTTGCTGACTAGCTTTGCGCATCTGCGTGATGCCGGTACAACCGCAGCTTCGTGCTGCATTTCCCCCGGTAGCTGGACAGATCAGGGCAACCACATGGCTAACCGCGAGAACTCCGACAACTCCCGTCTCGGGAATAc                                                                    >  2:405569/3‑149 (MQ=255)
                                                                                  gaaAGGTCAGTTGCTGAGTAGCTTTGCGCATCTGCGTGATGACGGTACAACCGCATCTTCTTGCTGGATCTACACCGGTAGCTGGACAGAGCAGGGCAACCAGATGGCTAACCGCGATAACTCCGACCCGTCCGGTCTGGGGAATACGc                                                                  <  1:280482/149‑1 (MQ=255)
                                                                                     aGGTCAGTTGCTGAGTAGCTTTGCGCATCTGCGTGATGACGGTACAACCGCATCTTCTTGCTGGATCTACACCGGTAGCTGGACAGAGCAGGGCAACCAGATGGCTACCCGCGATAACTCCGACCCGTCCGGTCTGGGGAATTCGCTgg                                                               >  1:1161168/1‑149 (MQ=255)
                                                                                     aGGTCAGTTGCTGAGTAGCTTTGCGCATCTGCGTGATGACGGTACAACCGCATCTTCTTGCTGGATCTACACCGGTAGCTGGACAGAGCAGGGCAACCAGATGGCTAACCGCGATAACTCCGACCCGTCCGGTCTGGGGAATACGCTgg                                                               >  2:213202/1‑149 (MQ=255)
                                                                                       gTCAGTTGCTGAGTAGCTTTGCGCATCTGCGTGATGACGGTACAAACGCATCTTCTTGCTGGATCTACACCGGTAGCTTGACAGAGCAGGGCAACCAGATGGCTAACCGCGATAACTCCGACCCGTCCGGTCTGGGGAATACTCTgaga                                                             >  1:301775/1‑149 (MQ=255)
                                                                                        tCAGTTGCTGAGTAGCTTTGCGCGTCTGCGTGATGACGGTACAACCGCATCTTCTTGCTGGATCTACACCGGTAGCTGGACAGAGCAGGGCAACCAGATGGCTAACCGCGATAACTCCGACCCGTCCGGTCTGGGGAATACGCTGGGAt                                                            >  2:203802/1‑149 (MQ=255)
                                                                                        tCAGTTGCTGAGTAGCTTTGCGCATCTGCGTGATGACGGTACTACCTCATCTTCTTGCTGGATCTACACCGGTAGCTGGACAGAGCAGGGCAACCAGATGGCTAACCGCGATAACTCCGACCCGTCCGGTCTGGGGAATACGCTGGGAt                                                            <  2:750902/149‑1 (MQ=255)
                                                                                        tCAGTTGCTGAGTAGCTTTGCGCATCTGCGTGATGACGGTACAACCGCATCTTCTTGCTGGATCTACACCGGTAGCTGGACAGAGCAGGGCAACCAGATGGCTAACCGCGATAACTCCGACCCGTCCGGTCTGGGGAATACGCTGGGAt                                                            >  1:939964/1‑149 (MQ=255)
                                                                                          aGTTGCTGAGTAGCGGTGCGCATCTGCGTGTTTATGTTAATTATTCTTCTTCTTGCGTGATCTACACCGGTAGCTGGACAGAGCAGGGCAACCAGATGGCTAACCGCGATAACTCCGACCCGTCCGGTCTGGGGAATACGCTGGGATgg                                                          <  2:203064/149‑1 (MQ=255)
                                                                                            ttGCTGAGTAGCTTTGCGCATCTGCGTGATGAGGGTACAACCGCATCTTCTTGCTGGATCTACACCGGTAGCTGGACAGAGCAGGGCAACCAGATGGCTAACCGCGATAACTCCGACCCGTCCGGTCTGGGGAATACGCTGGGATGGGc                                                        <  1:119089/149‑1 (MQ=255)
                                                                                                      gcagtGCGCAGCTGCGGGATGATGTTACATCAGCATCTTCTTGCTTGATCTACACCGGTAGCTGGACAGAGCAGGGCAACCAGATGGCTAACCGCGATAACTCCGACCCGTCCGGTCTGTGGAATACGCTGGGATTGGCCTGGGCGTgg                                              <  2:707596/145‑1 (MQ=255)
                                                                                                       cTTTGCGCATCTGCGTGATGACGGTACAACCGCATCTTCGTGCTGGATCTACAGCGGTAGCTGGACAGAGCAGGGGAACCAGATGGCTAACCGCGATAAATCCGACCCGTCCGGGGGGGGGTTTACGGTGGGCTGGGCCTGGGCGTgga                                             >  2:731501/1‑148 (MQ=255)
                                                                                                            cgcATCTGCGTGATGATGTTACTTCCGCATCTTCTTGCTGGATCTACACCGGTAGCTGGACAGAGCAGGGCAACCAGATGGCTAACCGCGATAACTCCGACCCGTCCGGTCTGGGGAATACGCTGGGATGGGCCTGGGCGTGGCCGCTc                                        <  2:523120/149‑1 (MQ=255)
                                                                                                               aTCTGCGTGATGTATGTATAAATGCTTCTTCTTGCTTGATCTACACCGGTAGCTGGACAGAGCAGGGCAACCAGATGGCTAACCGCGATAACTCCGACCCGTCCGGTCTGGGGAATACGCTGGGATGGGCCTGGGCGTGGCCGCTCAAc                                     <  2:933690/149‑1 (MQ=255)
                                                                                                               aTCTGCGTGATGACGGTACAACCGCATCTTCTTGCTGGATCTACACCGGTAGCTGGACAGAGCAGGGCAACCAGATGGCTTACCGCGATAACTCCGCCCCGTCAGGTCTGTGGAATAAGCTGGGATGGGCCTGGGAGTGGCAGcaaaaa                                     >  2:225189/1‑144 (MQ=255)
                                                                                                                   gCGTGATGACGGGACAACCGCATCTTCTTGCTGGATCAACCCCGGTAGCGGGACAGAGCACGGCAACCATATGGCTAACCGCGATAACTCCGACCCCGCCGGGCTGGGGAATACGCTGGGATGGGCCTGGGCGTGGCCGCTCAACCGTc                                 <  1:1018042/149‑1 (MQ=255)
                                                                                                                       gatgaCGGTACAACCGCATCTTCTTGCTGGATCTACACCGGTAGCTGGACAGAGCAGGGCAACCAGATGGCTAACCGCGATAACTACGACCCGTCCCGTCTGGGGAATACGCTGGGATGGGCCTCGGCGTGGCCGCTCAACCGTCGCGt                             >  1:628646/1‑149 (MQ=255)
                                                                                                                        atgaCGGTACAACCGCATCTTCTTGCTGGATCTACACCGGTAGCTGGACAGAGCAGGGCAACCAGATGGCTAACCGCGATAACTCCGACCCGTCCGGTCTGGGGAATACGCTGGGATGGGCCTGGGCGTGGCCGCTCAACCGTCGCGTg                            >  2:757149/1‑149 (MQ=255)
                                                                                                                        atgaCGGTACAACCGCATCTTCTTGCTGGATCTACACCGGTAGCTGGACAGAGCAGGGCAACCAGATGGCTAACCGCGATAACTCCGACCCGTCCGGTCTGGGGAATACGCTGGGATGGGCCTGGGCGTGGCCGCTCAACCGTCGCGTg                            >  1:1064004/1‑149 (MQ=255)
                                                                                                                          gaCGGTACAACCGCATCTTCTTGCTGGATCTACACCGGTAGCTGGACAGAGCAGGGCAACCAGATGGCTAACCGCGATAACTCCGACCCGTCCGGTCTGGGGAATACGCTGGGATGGGCCTGGGCGTGGCCGCTCAACCGTCGCGTGCt                          >  1:936509/1‑149 (MQ=255)
                                                                                                                           aCGGTACAACCGCATCTTCTTGCTGGATCTACACCGGTAGCTGGACAGAGCAGGGCAACCAGATGGCTAACCGCGATAACTCCGACCCGTCCGGTCTGGGGAATACGCTGGGATGGGCCTGGGCGTGGCCGCTCAACCGTCGCGTGCTg                         <  1:321167/149‑1 (MQ=255)
                                                                                                                                accaCCGCATCTTCCTGCTGGACCTCCCCCGGTAGCTGGCCAGAGCAGGGCAACCAGAGGGCTAACCGCGCTAACTCCGACCCGTCCGGTCTGGGGAATACGCTGGGATGGGCCTGGGCGTGGCCGCTCAACCGTCGCGTGCTGTACaa                    <  1:154157/146‑1 (MQ=255)
                                                                                                                                  aaCCGCATCTTCTTGCTGGATCTACACCGGTAGCTGGACAGAGCAGGGCAACCAGATGGCTAACCGCGATAACTCCGACACGTCCGGTCTGGGGAATACGCTGGGAAGGGCCTGGGAGTGGCCGCTCAACCGTCGCGTGCTGTACAAcc                  >  1:959056/1‑149 (MQ=255)
                                                                                                                                        aTCTTCTTGCTTGATCTACACCGGGAGCTGGACAGAGCAGGGCAACCAGATGGCTAACCGCGATAACTCCGACCCGTCCGGGCTGGGGAATACGCTGGGATGGGCCTGGGCGTGGCCGCTCAACCGTCGCGTGCTGTACAACCGTGCtt            <  1:424603/149‑1 (MQ=255)
                                                                                                                                              ttGCTGGATCTACACCGGTAGCTGGACAGAGCAGGGCAACCAGATGGCTAACCGCGATAACTCCGACCCGGCCGGTCTGGGGAATACGCTGGGATGGGCCTGGGCGTGGCCGCGCAACCGTCGCGTGCTGTACAACCGTGCTTcggcgg      <  1:331103/149‑1 (MQ=255)
                                                                                                                                                gCTGGATCTACACCGGTAGCTGGACAGAGCAGGGCAACCAGATGGCTAACCGCGATAACTCCGACCCGTCCGGTCTGGGGAATACGCTGGGATGGGCCTGGGCGTGGCCGCTCAACCGTCGCGTGCTGTACAACCGTGCTTCGGCGGat    >  2:188267/1‑149 (MQ=255)
                                                                                                                                                 cTGGATCTACACCGGTAGCTGGACAGAGCAGGGAAACCAGATGGCTAACCGCGATAACTCCGACCCGTCCGGTCTGGGGAATACGCTGGGATGGGCCTGGGCGTGGCCGCGCAACCGTCGCGTGCTGTACACCCGTGCTGCGGCGGata   >  2:325600/1‑149 (MQ=255)
                                                                                                                                                  tGGATCTACACCGGTAGCTGGACAGAGCAGGGCAACCAGATGGCTAACCGCGATAACTCCGACCCGTCCGGTCTGGGGAATACGCTGGGATGGGCCTGGGCGTGGCCGCTCAACCGTCGCGTGCTGTACAACCGTGCTTCGGCGGatat  >  1:391958/1‑149 (MQ=255)
                                                                                                                                                  tGGATCTACACCGGTAGCTGGACAGAGCAGGGCAACCAGATGGCTAACCGCGATAACTCCGACACGTCCGGTCTGGGGAATACGCTGGGATGGGCCTGGGCGTGGCCGCTCAACCGTCGCGAGCTGTACAACCGTGCTTCGGCGGatat  >  2:47278/1‑149 (MQ=255)
                                                                                                                                                  tGGATCTACACCGGTAGCTGGACAGAGCAGGGCAACCAGATGGCTAACAGCGATAACTCCGACCCGTCCGGTCTGGGGAATACGCTGGGATTGGCCTGGGCGTGGCCGCTCAACCGTCGCGTGCTGTACAACCGTGCTTCGGCGGatat  >  2:654723/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CGCAATCTGACGAAGTGGCTAAAGAGAACAACGGCTACGCGCTGGAAGATCTCTATGACGCTAATGGCGTGCTTATTGCGAAGAAAGGTCAGTTGCTGAGTAGCTTTGCGCATCTGCGTGATGACGGTACAACCGCATCTTCTTGCTGGATCTACACCGGTAGCTGGACAGAGCAGGGCAACCAGATGGCTAACCGCGATAACTCCGACCCGTCCGGTCTGGGGAATACGCTGGGATGGGCCTGGGCGTGGCCGCTCAACCGTCGCGTGCTGTACAACCGTGCTTCGGCGGATAT  >  NC_000913/1549463‑1549757

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑