Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,446,362 Δ1 bp 5.3% intergenic (+156/+16) ydbL → / ← feaR DUF1318 family protein/transcriptional activator for tynA and feaB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,446,3620A. 5.3% 315.2 / 9.4 94intergenic (+156/+16)ydbL/feaRDUF1318 family protein/transcriptional activator for tynA and feaB
Reads supporting (aligned to +/- strand):  ref base A (44/45);  new base . (3/2);  total (47/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.61e-01

GCAATATCGGGATAACGCTACAGGAGAAGGCGATTGATCTATTTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGGCATG  >  NC_000913/1446217‑1446509
                                                                                                                                                  |                                                                                                                                                   
gCAATATCGGGATAACGCTACAGGAGAAGGCGATTGATCTATTTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATg                                                                                                                                                   >  3:208051/1‑148 (MQ=255)
 cAATATCGGGATAACGCTACAGGAGAAGGCGATTGATCTATTTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGc                                                                                                                                                  >  2:311178/1‑148 (MQ=255)
     atCGGGATAACGCTACAGGAGAAGGCGATTGATCTATTTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGc                                                                                                                                             >  2:299463/1‑149 (MQ=255)
     atCGGGATAACGCTACAGGAGAAGGCGATTGATCTATTTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGc                                                                                                                                             >  3:175047/1‑149 (MQ=255)
      tCGGGATAACGCTACAGGAGAAGGCGATTGATCTATTTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTACTGCAACGc                                                                                                                                             <  4:61409/148‑1 (MQ=255)
         ggATAACGCTACAGGAGAAGGCGATTGATCTATTTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAAt                                                                                                                                         >  1:218408/1‑149 (MQ=255)
         ggATAACGCTACAGGAGAAGGCGATTGATCTATTTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAAt                                                                                                                                         >  6:60764/1‑149 (MQ=255)
            tAACGCTACAGGAGAAGGCGATTGATCTATTTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAAttag                                                                                                                                      <  2:305239/149‑1 (MQ=255)
              aCGCTACAGGAGAAGGCGATTGATCTATTTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAAttagt                                                                                                                                     <  2:302675/148‑1 (MQ=255)
                gCTACAGGAGAAGGCGATTGATCTATTTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAAttagttag                                                                                                                                  >  5:258450/1‑149 (MQ=255)
                 cTACAGGAGAAGGCGATTGATCTATTTTCCTGAAACAAGGTGAATATTCAAACACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGACAATTAGTTAGc                                                                                                                                 <  6:329788/149‑1 (MQ=255)
                  tACAGGAGAAGGCGATTGATCTATTTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAAttagtt                                                                                                                                    >  8:248208/1‑145 (MQ=255)
                  tACAGGAGAAGGCGATTGATCTATTTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAAttagtt                                                                                                                                    <  7:248208/145‑1 (MQ=255)
                   aCAGGAGAAGGCGATTGATCTATTTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCg                                                                                                                                >  2:256654/1‑148 (MQ=255)
                    cAGGAGAAGGCGATTGATCTATTTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCg                                                                                                                                <  2:48296/147‑1 (MQ=255)
                     aGGAGAAGGCGATTGATCTATTTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGaa                                                                                                                             <  6:305633/149‑1 (MQ=255)
                      ggAGCACGCGATTGATCTCTTTTCCTGAAACCAGCTGAATATTCAAAAACTCCTGTCAAATTGCCTTATGCCCTGAACAATGCCTAGTCAAATCTCTGGTTTCACTACGTGGCGTGGCCGTACATGCAACGCCAATTAGTTAGCGGAAt                                                                                                                            <  8:265992/149‑1 (MQ=255)
                      ggAGAAGGCGATTGATCTATTTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAAt                                                                                                                            <  3:32666/149‑1 (MQ=255)
                       gagaAGGCGATTGATCTATTTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAAtt                                                                                                                           <  7:108381/149‑1 (MQ=255)
                       gagaAGGCGATTGATCTATTTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAAtt                                                                                                                           >  8:194345/1‑149 (MQ=255)
                         gaAGGCGATTGATCTATTTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTa                                                                                                                         <  4:208051/149‑1 (MQ=255)
                          aaGGCGATGGATCTATTTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTATTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTAc                                                                                                                        <  5:278234/149‑1 (MQ=255)
                            ggCGATTGATCTATTTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCaa                                                                                                                                          >  1:155910/1‑129 (MQ=255)
                            ggCGATTGATCTATTTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCaa                                                                                                                                          <  2:155910/129‑1 (MQ=255)
                            ggCGATTGATCTATTTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGt                                                                                                                      <  1:189761/149‑1 (MQ=255)
                            ggCGATTGATCTATTTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGt                                                                                                                      <  4:175047/149‑1 (MQ=255)
                                 ttGATCTATTTTCCTGAAACAATGTGAATATTCAAAAACTCCTGTCAAATTTCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGa                                                                                                                   <  3:91496/147‑1 (MQ=255)
                                  tGATCTATTTTCCTGAAACAAGCTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGACAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCAAATTACTTAGCGGAATATACTTCGATa                                                                                                                 >  7:89108/1‑148 (MQ=255)
                                     tCTATTTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTACTGCAACGCCAAttagttaa                                                                                                                                  <  8:130968/128‑2 (MQ=255)
                                     tCTATTTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTACTGCAACGCCAAttagttaa                                                                                                                                  >  7:130968/1‑127 (MQ=255)
                                     tCTATTTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCg                                                                                                             >  2:7352/1‑149 (MQ=255)
                                         ttttCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTc                                                                                                              <  7:151825/144‑1 (MQ=255)
                                          tttCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTATTAACTAATAGGCGTTGCAGTACATGAAACGACAATTAGTTAGCGGATTTTACGTCGATACTCGACTgg                                                                                                        >  3:41134/1‑149 (MQ=255)
                                              cTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGt                                                                                                     >  1:143868/1‑148 (MQ=255)
                                               tGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCa                                                                                                   >  7:11464/1‑149 (MQ=255)
                                                gAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCAt                                                                                                  >  7:335407/1‑149 (MQ=255)
                                                 aaaCAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATc                                                                                                 >  7:207577/1‑149 (MQ=255)
                                                 aaaCAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATc                                                                                                 <  8:263691/149‑1 (MQ=255)
                                                 aaaCAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATc                                                                                                 <  8:212947/149‑1 (MQ=255)
                                                 aaaCAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATc                                                                                                 >  7:212947/1‑149 (MQ=255)
                                                  aaCAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCaa                                                                                                                                          <  4:207567/107‑1 (MQ=255)
                                                  aaCAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCaa                                                                                                                                          >  3:207567/1‑107 (MQ=255)
                                                  aaCAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATc                                                                                                 <  8:207577/148‑1 (MQ=255)
                                                   aCAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATccc                                                                                               <  1:256654/149‑1 (MQ=255)
                                                               ttCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCtt                                                                                   <  2:143868/149‑1 (MQ=255)
                                                                  aaaaaCTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTTACTAATTGGCGTTGCAGTACATGGAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGTCGTCATCCCAAAGCGTTGCTTaaa                                                                                <  4:124460/149‑1 (MQ=255)
                                                                            gTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGa                                                                       >  8:167414/1‑148 (MQ=255)
                                                                             tCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGTAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGaa                                                                      >  1:254505/1‑148 (MQ=255)
                                                                             tCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGa                                                                       <  5:60764/147‑1 (MQ=255)
                                                                                taTTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAAt                                                                  <  8:95880/148‑1 (MQ=255)
                                                                                aaTTGCCTTTTGCCCTGAAAAATGCATAGTCATATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAAt                                                                  >  2:161034/1‑149 (MQ=255)
                                                                                aaTTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAAt                                                                  >  8:180945/1‑149 (MQ=255)
                                                                                aaTTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAAt                                                                  <  7:167414/149‑1 (MQ=255)
                                                                                     ccTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGActc                                                             <  2:218408/149‑1 (MQ=255)
                                                                                       ttttGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTg                                                           >  1:1925/1‑149 (MQ=255)
                                                                                         ttGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTg                                                           <  8:335407/147‑1 (MQ=255)
                                                                                          tGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGt                                                         <  2:332337/148‑1 (MQ=255)
                                                                                            cccTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGt                                                         >  5:188132/1‑146 (MQ=255)
                                                                                               tGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACTCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGa                                                                <  8:303695/136‑1 (MQ=255)
                                                                                               tGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATAc                                                                                                                >  8:354416/1‑88 (MQ=255)
                                                                                               tGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGa                                                                >  7:303695/1‑136 (MQ=255)
                                                                                               tGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGCAATTTACGTCGATAc                                                                                                                <  7:354416/88‑1 (MQ=255)
                                                                                                 aaaaaTGCATAGTCAAATATCTGTTTTAACTTATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGaaa                                                   <  7:180945/147‑1 (MQ=255)
                                                                                                    aaTGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATccc                                              >  7:133301/1‑149 (MQ=255)
                                                                                                        cATAGTCAAATATCTGTTTTAACTAATTGGCGTTGAAGTACATGCAACGCCAATTAGTTAGAGGAATTTACGTCGATACTCGACTGGCGTCATCCCAAAGCGTTGCTTAAATACGGTAGAAAAATGACTCTGGTCAGAAAATCCACAAt                                          >  2:308935/1‑149 (MQ=255)
                                                                                                            gTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGaaa                                      <  1:7352/149‑1 (MQ=255)
                                                                                                             tCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATa                                                                                                                 <  5:173978/73‑1 (MQ=255)
                                                                                                             tCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATa                                                                                                                 >  6:173978/1‑73 (MQ=255)
                                                                                                             tCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAg                                     >  4:55280/1‑149 (MQ=255)
                                                                                                               aaaTATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGaa                                       <  3:159808/145‑1 (MQ=255)
                                                                                                               aaaTATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGaa                                       >  4:159808/1‑145 (MQ=255)
                                                                                                               aaaTATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGcc                                   >  6:227300/1‑149 (MQ=255)
                                                                                                                aaTATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAAttagttag                                                                                                                                  <  4:75964/53‑1 (MQ=255)
                                                                                                                aaTATCTGTTTTAACTAATTGGCGTAGCAGTACATGCAACGCCAAttagttag                                                                                                                                  >  3:75964/1‑53 (MQ=255)
                                                                                                                    tCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGAtg                               <  1:311178/148‑1 (MQ=255)
                                                                                                                      tccgctAACTAATTGGCGTTGCATGTACTGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGAtgcct                            >  2:228039/6‑149 (MQ=38)
                                                                                                                       gTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAAt                                                                  >  5:155843/1‑110 (MQ=255)
                                                                                                                       gTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAAt                                                                  <  6:155843/110‑1 (MQ=255)
                                                                                                                        ttttAACTAATTGGCGTTGCAGTACATGCAACGCGAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGAtgc                              <  1:299463/145‑1 (MQ=255)
                                                                                                                          ttAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGAtgcctgc                          >  7:249698/1‑147 (MQ=255)
                                                                                                                           tAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGtt                                                                         >  8:1450/1‑99 (MQ=255)
                                                                                                                           tAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGtt                                                                         <  7:1450/99‑1 (MQ=255)
                                                                                                                             aCTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGtt                     <  6:188132/149‑1 (MQ=255)
                                                                                                                                 aTTGGCGTTGCATGTACTGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGttttt                  >  8:355403/1‑148 (MQ=255)
                                                                                                                                 aTTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGttttt                  >  2:148932/1‑148 (MQ=255)
                                                                                                                                 aTTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGttttt                  <  7:194345/148‑1 (MQ=255)
                                                                                                                                     gCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGTTTTtcatca             >  4:179983/1‑149 (MQ=255)
                                                                                                                                     gCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGTTTTtcatca             >  3:192236/1‑149 (MQ=255)
                                                                                                                                       gTTTCAGTAAATGCAAAGCCAATTAGTTAGCGGAATTTACGTAGATACTCGCCTGGCGTAATCCCAATGCGTAGATTAAATACACTTGAAAAAATACTCGGGTAAGAAAATCCCCAATGAAAGCCGCTGCCTGCCAGTTATTCatcatc           >  4:99390/1‑149 (MQ=255)
                                                                                                                                               tACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGGCa    >  7:74435/1‑149 (MQ=255)
                                                                                                                                               tACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGGCa    >  4:264388/1‑149 (MQ=255)
                                                                                                                                                 cATGCAACGCCAATTAGTTATCGTAATTTACGTCGATACTGGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGGCGATGCCTGCCAGTTTTTCATCATCTGCGGCAt   <  4:41134/148‑1 (MQ=255)
                                                                                                                                                 cATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGGCATg  <  8:11464/149‑1 (MQ=255)
                                                                                                                                                 cATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGGCATg  <  2:1925/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GCAATATCGGGATAACGCTACAGGAGAAGGCGATTGATCTATTTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGGCATG  >  NC_000913/1446217‑1446509

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