Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,986,356 A→T 5.6% intergenic (‑228/‑569) araF ← / → ftnB L‑arabinose ABC transporter periplasmic binding protein/ferritin B, putative ferrous iron reservoir

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,986,3560AT 5.6% 427.3 / 10.7 121intergenic (‑228/‑569)araF/ftnBL‑arabinose ABC transporter periplasmic binding protein/ferritin B, putative ferrous iron reservoir
Reads supporting (aligned to +/- strand):  ref base A (49/64);  new base T (3/4);  total (52/69)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 125 total aligned reads displayed.

TACCTTTTAACTAAAAGATAAGTGACTGTGTTGACATAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGCGTTGAGTTTTTGAAATGGTGTA  >  NC_000913/1986208‑1986497
                                                                                                                                                    |                                                                                                                                             
tACCTTTTAACTAAAAGATAAGTGACTGTGTTGACATAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCa                                                                                                                                               <  2:952563/149‑1 (MQ=255)
  ccTTTTAACTAAAAGATAAGTGACTGTGTTGACATAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGa                                                                                                                                             >  2:963213/1‑149 (MQ=255)
     tttAACTAAAAGATAAGTGACTGTGTTGACATAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACa                                                                                                                                           >  4:92984/1‑148 (MQ=255)
         aCTAAAAGATAAGTGACTGTGTTGACATAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCTGACAAATTg                                                                                                                                      <  1:97399/149‑1 (MQ=255)
          cTAAAAGATAAGTGACTGTGTTGACATAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCTGACAAATTgg                                                                                                                                     <  1:102317/149‑2 (MQ=255)
          cTAAAAGATAAGTGACTGTGTTGACATAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCTGACAAATTgg                                                                                                                                     <  2:660550/149‑2 (MQ=255)
           tAAAAGATAAGTGACTGTGTTGACATAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGCACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGt                                                                                                                                     <  3:612064/148‑1 (MQ=255)
           tAAAAGATAAGTGACTGTGTTGACATAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGt                                                                                                                                     <  4:488538/148‑1 (MQ=255)
           tAAAAGATAAGTGACTGTGTTGACATAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTc                                                                                                                                    >  1:843233/1‑149 (MQ=255)
           tAAAAGATAAGTGACTGTGTTGACATAGTTTTAGCGAGAAATTAATTCTCCACAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCCGACAAACTGt                                                                                                                                     <  4:946265/148‑1 (MQ=255)
            aaaaGATAAGTGACTGTGTTGACATAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTc                                                                                                                                    >  4:523466/1‑148 (MQ=255)
            aaaaGATAAGTGACTGTGTTGACATAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTc                                                                                                                                    >  4:467874/1‑148 (MQ=255)
            aaaaGATAAGTGACTGTGTTGACATAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCg                                                                                                                                   >  2:566695/1‑149 (MQ=255)
             aaaGATAAGTGACTGTGTTGACATAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGa                                                                                                                                  >  1:227241/1‑149 (MQ=255)
             aaaGATAAGTGACTGTGTTGACATAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGa                                                                                                                                  >  1:181747/1‑149 (MQ=255)
             aaaGATAAGTGACTGTGTTGACATAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGa                                                                                                                                  <  1:467713/149‑1 (MQ=255)
              aaGATAAGTGACTGTGTTGACATAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGaa                                                                                                                                 <  1:933357/149‑1 (MQ=255)
              aaGATAAGTGACTGTGTTGACATAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTATCGaa                                                                                                                                 >  1:862305/1‑149 (MQ=255)
               aGATAAGTGACTGTGTTGACATAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAAt                                                                                                                                >  1:685138/1‑149 (MQ=255)
                    aGTGACTGTGTTGACATAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGcac                                                                                                                            <  4:26336/148‑1 (MQ=255)
                     gTGACTGTGTTGACATAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCAcag                                                                                                                          >  2:369961/1‑149 (MQ=255)
                       gACTGTGTTGACATAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCAcagca                                                                                                                        <  1:45677/149‑1 (MQ=255)
                           gtgtTGACATAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGAt                                                                                                                     >  3:527102/1‑148 (MQ=255)
                             gtTGACATAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTa                                                                                                                   <  3:396480/148‑1 (MQ=255)
                              ttGACATAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAAt                                                                                                                 <  1:510728/149‑1 (MQ=255)
                               tGACATAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATc                                                                                                                <  1:418155/149‑1 (MQ=255)
                                   aTAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATa                                                                                                            >  1:322751/1‑149 (MQ=255)
                                    tAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCa                                                                                                                                               >  1:716782/1‑113 (MQ=255)
                                    tAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCa                                                                                                                                               <  2:716782/113‑1 (MQ=255)
                                    tAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATaa                                                                                                           <  2:322751/149‑1 (MQ=255)
                                    tAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATaa                                                                                                           <  2:897275/149‑1 (MQ=255)
                                    tAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATa                                                                                                            >  3:827531/1‑148 (MQ=255)
                                    tAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATa                                                                                                            <  3:818081/148‑1 (MQ=255)
                                    tAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATa                                                                                                            >  3:127461/1‑148 (MQ=255)
                                          tAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTa                                                                                                      <  4:827531/148‑1 (MQ=255)
                                             cGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCt                                                                                                  <  2:436916/149‑1 (MQ=255)
                                              gagaAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTg                                                                                                 >  1:580762/1‑149 (MQ=255)
                                                gaAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGa                                                                                               >  2:543180/1‑149 (MQ=255)
                                                  aattaattCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCTGACAAAtt                                                                                                                                       <  3:313762/107‑1 (MQ=255)
                                                  aattaattCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCTGACAAAtt                                                                                                                                       >  4:313762/1‑107 (MQ=255)
                                                   attaattCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAAt                                                                                            >  1:582716/1‑149 (MQ=255)
                                                      aattCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCt                                                                                         >  2:169807/1‑149 (MQ=255)
                                                      aattCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCt                                                                                         <  2:843233/149‑1 (MQ=255)
                                                        ttCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCttg                                                                                       <  2:823183/149‑1 (MQ=255)
                                                        ttCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCttg                                                                                       >  1:869549/1‑149 (MQ=255)
                                                        ttCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAACCCttg                                                                                       <  1:369961/149‑1 (MQ=255)
                                                          ctcCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCttgt                                                                                      <  3:468856/148‑1 (MQ=255)
                                                               tAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCtt                                                                                >  1:6943/1‑149 (MQ=255)
                                                                aGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCttt                                                                               <  2:1074741/149‑1 (MQ=255)
                                                                 ggAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCTGACAAATTCTCGAATGCTCAGCAGATTAATCCATAAGATTAGCCTGGTAATCCTTGTTGTCttt                                                                               >  4:631212/1‑148 (MQ=255)
                                                                 ggAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCttt                                                                               <  3:331071/148‑1 (MQ=255)
                                                                  gagaGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGa                                                                                               <  2:21259/131‑1 (MQ=255)
                                                                  gagaGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGa                                                                                               >  1:21259/1‑131 (MQ=255)
                                                                  gagaGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTgg                                                                             <  2:869549/149‑1 (MQ=255)
                                                                  gagaGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTg                                                                              <  4:127461/148‑1 (MQ=255)
                                                                     agCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGca                                                                                                                             <  1:721230/98‑1 (MQ=255)
                                                                     agCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCAcag                                                                                                                          >  2:414591/1‑101 (MQ=255)
                                                                        aaTATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCa                                                                       <  1:417505/149‑1 (MQ=255)
                                                                          tatCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCAt                                                                      >  4:531119/1‑148 (MQ=255)
                                                                              acaTCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCgg                                                                  <  3:438236/148‑1 (MQ=255)
                                                                               caTCGCAGAATTACAGTGAGAACGTGCAAAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAACTGTCGAATGCATAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCggg                                                                 >  3:720337/1‑148 (MQ=255)
                                                                                aTCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAAACCTTGTTGTCTTTGGTACCCATGCGGGa                                                                <  3:631212/148‑1 (MQ=255)
                                                                                   gCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGtc                                                            <  2:546015/149‑1 (MQ=255)
                                                                                      gAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGtctt                                                          <  3:252148/148‑1 (MQ=255)
                                                                                        aTTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGtcttct                                                        <  3:241154/148‑1 (MQ=255)
                                                                                        aTTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTcttctt                                                       <  2:1044691/149‑1 (MQ=255)
                                                                                        aTTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTcttctt                                                       <  1:722372/149‑1 (MQ=255)
                                                                                           aCAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCtttt                                                     <  4:212321/148‑1 (MQ=255)
                                                                                           aCAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCtttt                                                     <  3:626714/148‑1 (MQ=255)
                                                                                               tGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAAcc                                                >  2:824844/1‑149 (MQ=255)
                                                                                                  gaACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATcc                                                                                          >  1:366646/1‑104 (MQ=255)
                                                                                                  gaACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGt                                             <  2:685138/149‑1 (MQ=255)
                                                                                                       tGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAAt                                         <  3:1031294/148‑1 (MQ=255)
                                                                                                        gCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAg                                       >  1:722446/1‑149 (MQ=255)
                                                                                                         cATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAgg                                      <  2:582716/149‑1 (MQ=255)
                                                                                                         cATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAg                                       <  3:392978/148‑1 (MQ=255)
                                                                                                           tAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCt                                                                                 <  2:11002/104‑1 (MQ=255)
                                                                                                           tAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCt                                                                                 >  1:11002/1‑104 (MQ=255)
                                                                                                             aaTTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCg                                   <  3:115035/148‑1 (MQ=255)
                                                                                                                ttAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCAt                                >  4:75169/1‑148 (MQ=255)
                                                                                                                   gCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTAcc                            <  1:814185/149‑1 (MQ=255)
                                                                                                                     gggAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTg                          >  1:1045601/1‑149 (MQ=255)
                                                                                                                       gAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGc                        >  2:421438/1‑149 (MQ=255)
                                                                                                                          aaGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGCGtt                     <  1:563246/149‑1 (MQ=255)
                                                                                                                                tAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGCGTTGAGtt                >  3:453462/1‑148 (MQ=255)
                                                                                                                                tAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGCGTTGAGtt                <  4:203393/148‑1 (MQ=255)
                                                                                                                                 aaGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGCGTTGAGttt               >  4:882627/1‑148 (MQ=255)
                                                                                                                                  aGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGCGTTGAGttttt             <  2:181747/149‑1 (MQ=255)
                                                                                                                                     gAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGCGTTGAGTTTTTGa           >  4:645024/1‑148 (MQ=255)
                                                                                                                                     gAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCGTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGCGTTGAGTTTTTGaa          >  1:1021284/1‑149 (MQ=255)
                                                                                                                                      aaaGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGCGTTGAGTTTTTGaa          >  3:109336/1‑148 (MQ=255)
                                                                                                                                       aaGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGCGTTGAGTTTTTGaaa         <  3:488203/148‑1 (MQ=255)
                                                                                                                                        aGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATGAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTAGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGCGTTGAGTTTTTGAAATg       <  2:722036/149‑1 (MQ=255)
                                                                                                                                         gacAATTTGTCTGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGCGTTGAGTTTTTGAAATgg      >  1:961817/3‑149 (MQ=255)
                                                                                                                                         gCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGCGTTAAGTTTTTGAAATg       <  3:75169/148‑1 (MQ=255)
                                                                                                                                          ccAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGCGTTGAGTTTTTGAAATGgt     <  1:170018/149‑1 (MQ=255)
                                                                                                                                           cAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGCGTTGAGTTTTTGAAATGgt     <  3:531119/148‑1 (MQ=255)
                                                                                                                                            aaTTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGCGTTGAGTTTTTGAAATGgtgt   >  1:348419/1‑149 (MQ=255)
                                                                                                                                             aTTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGCGTTGAGTTTTTGAAATGGTGta  >  2:552543/1‑149 (MQ=255)
                                                                                                                                              tttGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGCGTTGAGTTTTTGAAATGGTGta  >  4:757709/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                             
TACCTTTTAACTAAAAGATAAGTGACTGTGTTGACATAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGCGTTGAGTTTTTGAAATGGTGTA  >  NC_000913/1986208‑1986497

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