Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,986,356 A→T 5.8% 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.8% 361.8 / 9.0 103intergenic (‑228/‑569)araF/ftnBL‑arabinose ABC transporter periplasmic binding protein/ferritin B, putative ferrous iron reservoir
Reads supporting (aligned to +/- strand):  ref base A (48/49);  new base T (2/4);  total (50/53)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

Only 100 of 104 total aligned reads displayed.

TACCTTTTAACTAAAAGATAAGTGACTGTGTTGACATAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGCGTTGAGTTTTTGAAATGGTGTA  >  NC_000913/1986208‑1986497
                                                                                                                                                    |                                                                                                                                             
    ttttAACTAAAAGATAAGTGACTGTGTTGACATAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACa                                                                                                                                           >  3:43925/1‑149 (MQ=255)
    ttttAACTAAAAGATAAGTGACTGTGTTGACATAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACa                                                                                                                                           >  1:796012/1‑149 (MQ=255)
    ttttAACTAAAAGATAAGTGACTGTGTTGACATAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACa                                                                                                                                           >  1:619728/1‑149 (MQ=255)
     tttAACTAAAAGATAAGTGACTGTGTTGACATAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACaa                                                                                                                                          >  1:864393/1‑149 (MQ=255)
         aCTAAAAGATAAGTGACTGTGTTGACATAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCTGACAAATTg                                                                                                                                      <  3:136335/149‑1 (MQ=255)
           tAAAAGATAAGTGACTGTGTTGACATAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATACATCTAGCCGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTc                                                                                                                                    <  4:206353/149‑1 (MQ=255)
           tAAAAGATAAGTGACTGTGTTGACATAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTc                                                                                                                                    <  4:664282/149‑1 (MQ=255)
           tAAAAGATAAGTGACTGTGTTGACATAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTc                                                                                                                                    >  1:198546/1‑149 (MQ=255)
           tAAAAGATAAGTGACTGTGTTGACATAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTc                                                                                                                                    >  3:191885/1‑149 (MQ=255)
            aaaaGATAAGTGACTGTGTTGACATAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCTGACAAATTggct                                                                                                                                   <  2:754310/149‑4 (MQ=255)
             aaaGATAAGTGACTGTGTTGACATAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGa                                                                                                                                  >  4:216547/1‑149 (MQ=255)
              aaGATAAGTGACTGTGTTGACATAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGaa                                                                                                                                 >  4:109587/1‑149 (MQ=255)
                  tAAGTGACTGTGTTGACATAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGca                                                                                                                             <  1:195535/149‑1 (MQ=255)
                     gTGACTGTGTTGACATAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCAcag                                                                                                                          <  2:619728/149‑1 (MQ=255)
                     gTGACTGTGTTGACATAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCAcag                                                                                                                          >  2:425732/1‑149 (MQ=255)
                     gTGACTGTGTTGACATAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCAcag                                                                                                                          >  3:330376/1‑149 (MQ=255)
                        aCTGTGTTGACATAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCAcagcag                                                                                                                       <  1:412706/149‑1 (MQ=255)
                              ttGACATAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAAt                                                                                                                 <  3:217337/149‑1 (MQ=255)
                               tGACATAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATc                                                                                                                <  2:574045/149‑1 (MQ=255)
                                gACATAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATcc                                                                                                               <  3:420395/149‑1 (MQ=255)
                                 aCATAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCa                                                                                                              >  2:867166/1‑149 (MQ=255)
                                   aTAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATa                                                                                                            >  1:680902/1‑149 (MQ=255)
                                    tAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATaa                                                                                                           <  2:159858/149‑1 (MQ=255)
                                    tAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCAGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATaa                                                                                                           <  3:642208/149‑1 (MQ=255)
                                       ttttAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAAt                                                                                                                 <  1:48788/140‑1 (MQ=255)
                                       ttttAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAAt                                                                                                                 >  2:48788/1‑140 (MQ=255)
                                             cGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCt                                                                                                  >  3:65765/1‑149 (MQ=255)
                                              gagaAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTg                                                                                                 >  2:517143/1‑149 (MQ=255)
                                                 aaaTTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGaa                                                                                              <  1:246051/149‑1 (MQ=255)
                                                  aattaattCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGaaa                                                                                             <  1:517143/149‑1 (MQ=255)
                                                    ttaattCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATc                                                                                           <  1:325844/149‑1 (MQ=255)
                                                    ttaattCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATc                                                                                           >  2:477988/1‑149 (MQ=255)
                                                     taattCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATcc                                                                                          <  3:561271/149‑1 (MQ=255)
                                                     taattCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATcc                                                                                          <  2:864393/149‑1 (MQ=255)
                                                        ttCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCttg                                                                                       <  3:383620/149‑1 (MQ=255)
                                                         tctcCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCttgt                                                                                      >  4:449293/1‑149 (MQ=255)
                                                         tctcCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCttgt                                                                                      <  3:216547/149‑1 (MQ=255)
                                                         tctcCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCttgt                                                                                      <  1:867166/149‑1 (MQ=255)
                                                             cATAGGAGAGCAATATCACATCGCAGAATTACAGTCAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTc                                                                                  >  1:705222/1‑149 (MQ=255)
                                                               tAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCtt                                                                                >  3:100563/1‑149 (MQ=255)
                                                                aGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCttt                                                                               >  1:212237/1‑149 (MQ=255)
                                                                              acaTCGCAGACTTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCggg                                                                 >  2:590795/1‑149 (MQ=255)
                                                                               caTCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGa                                                                <  3:279119/149‑1 (MQ=255)
                                                                               caTCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGa                                                                <  1:425732/149‑1 (MQ=255)
                                                                                aTCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGAt                                                               >  3:532370/1‑149 (MQ=255)
                                                                                      gAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGtcttc                                                         <  1:195986/149‑1 (MQ=255)
                                                                                      gAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGtcttc                                                         <  2:390974/149‑1 (MQ=255)
                                                                                       aaTTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATa                                                                                                            <  4:140697/97‑1 (MQ=255)
                                                                                       aaTTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATa                                                                                                            >  3:140697/1‑97 (MQ=255)
                                                                                       aaTTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGtcttct                                                        <  4:752545/149‑1 (MQ=255)
                                                                                          tACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCtttt                                                     <  1:873812/149‑1 (MQ=255)
                                                                                            cAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTa                                                   >  4:63845/1‑149 (MQ=255)
                                                                                             aGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTaa                                                  >  1:415964/1‑149 (MQ=255)
                                                                                              gTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAAc                                                 >  1:150224/1‑149 (MQ=255)
                                                                                               tGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAAcc                                                >  2:248957/1‑149 (MQ=255)
                                                                                                  gaACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGt                                             >  2:492637/1‑149 (MQ=255)
                                                                                                       tGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATa                                        >  3:571791/1‑149 (MQ=255)
                                                                                                        gCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAg                                       <  2:514157/149‑1 (MQ=255)
                                                                                                           tAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGcc                                    >  4:107816/1‑149 (MQ=255)
                                                                                                              aTTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCa                                 >  4:671926/1‑149 (MQ=255)
                                                                                                                     gggAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCttttt                                                    >  4:756723/1‑123 (MQ=255)
                                                                                                                     gggAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCttttt                                                    <  3:756723/123‑1 (MQ=255)
                                                                                                                       gAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGc                        <  4:205478/149‑1 (MQ=255)
                                                                                                                         aaaGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGCGt                      >  4:377090/1‑149 (MQ=255)
                                                                                                                         aaaGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGCGt                      >  4:331745/1‑149 (MQ=255)
                                                                                                                          aaGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGCGtt                     <  2:413781/149‑1 (MQ=255)
                                                                                                                            gACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGCGTTGa                   >  4:566437/1‑149 (MQ=255)
                                                                                                                            gACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGCGTTGa                   <  4:300875/149‑1 (MQ=255)
                                                                                                                               aTAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGCGTTGAGtt                >  1:486353/1‑149 (MQ=255)
                                                                                                                                tAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGCGTTGAGttt               >  4:393822/1‑149 (MQ=255)
                                                                                                                                tAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGCGTTGAGttt               >  4:808076/1‑149 (MQ=255)
                                                                                                                                  aGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACTCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGCGTTGAGttttt             <  3:19520/149‑1 (MQ=255)
                                                                                                                                  aGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGCGTTGAGttttt             <  3:808076/149‑1 (MQ=255)
                                                                                                                                  aGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGCGTTGAGttttt             <  1:266349/149‑1 (MQ=255)
                                                                                                                                  aGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGCGTTGAGttttt             <  3:331745/149‑1 (MQ=255)
                                                                                                                                  aGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGCGTTGAGttttt             <  3:377090/149‑1 (MQ=255)
                                                                                                                                    ggAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGCGTTGAGTTTTTGa           >  3:704493/1‑149 (MQ=255)
                                                                                                                                    ggAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGCGTTGAGTTTTTGa           >  1:56608/1‑149 (MQ=255)
                                                                                                                                    ggAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGCGTTGAGTTTTTGa           <  2:486353/149‑1 (MQ=255)
                                                                                                                                    ggAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGCGTTGAGTTTTTGa           >  4:329205/1‑149 (MQ=255)
                                                                                                                                     gAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGCGTTGAGTTTTTGaa          >  4:726415/1‑149 (MQ=255)
                                                                                                                                      aaaGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCgg                                                                  >  1:233231/1‑92 (MQ=255)
                                                                                                                                      aaaGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCgg                                                                  <  2:233231/92‑1 (MQ=255)
                                                                                                                                      aaaGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGCGTTGAGTTTTTGaaa         <  1:272206/149‑1 (MQ=255)
                                                                                                                                        aGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGCGTTGAGTTTTTGAAATg       >  1:398728/1‑149 (MQ=255)
                                                                                                                                        aGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGCGTTGAGTTTTTGAAATg       <  3:329205/149‑1 (MQ=255)
                                                                                                                                        aGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGCGTTGAGTTTTTGAAATg       <  1:248957/149‑1 (MQ=255)
                                                                                                                                        aGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGCGTAGAGTTTTTGAAATg       <  3:63845/149‑1 (MQ=255)
                                                                                                                                         gacAATTTGTCTGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGAt                                                               >  4:566499/3‑92 (MQ=255)
                                                                                                                                         gacAATTTGTCTGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGAt                                                               <  3:566499/90‑1 (MQ=255)
                                                                                                                                         gCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAACATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGCGTTGAGTTTTTGAAATgg      <  3:726415/149‑1 (MQ=255)
                                                                                                                                          ccAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGCGTTGAGTTTTTGAAATGgt     <  3:55848/149‑1 (MQ=255)
                                                                                                                                          ccAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGCGTTGAGTTTTTGAAATGgt     >  2:421436/1‑149 (MQ=255)
                                                                                                                                          ccAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGCGTTGAGTTTTTGAAATGgt     <  4:151739/149‑1 (MQ=255)
                                                                                                                                          ccAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGCGTTGAGTTTTTGAAATGgt     <  1:590795/149‑1 (MQ=255)
                                                                                                                                          acAATTTGTCTGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGCGTTGAGTTTTTGAAATGgt     >  1:715792/2‑149 (MQ=255)
                                                                                                                                            aaTTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAg                                       >  3:529431/1‑113 (MQ=255)
                                                                                                                                            aaTTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAg                                       <  4:529431/113‑1 (MQ=255)
                                                                                                                                            aaTTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGCGTTGAGTTTTTGAAATGgtgt   <  4:43925/149‑1 (MQ=255)
                                                                                                                                             aTTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGCGTTGAGTTTTTGAAATGGTGta  <  3:566437/149‑1 (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: