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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,986,3560AT 5.7% 401.2 / 9.6 123intergenic (‑228/‑569)araF/ftnBarabinose ABC transporter periplasmic binding protein/putative ferritin‑like protein
Reads supporting (aligned to +/- strand):  ref base A (60/56);  new base T (4/3);  total (64/59)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.68e-01

Only 100 of 125 total aligned reads displayed.

TTACCTTTTAACTAAAAGATAAGTGACTGTGTTGACATAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGCGTTGAGTTTTTGAAATGGTGTAATAACCGC  >  NC_000913/1986207‑1986505
                                                                                                                                                     |                                                                                                                                                     
 tACCTTTTAACTAAAAGATAAGTGACTGTGTTGACATAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAg                                                                                                                                                      <  1:958247/150‑1 (MQ=255)
 tACCTTTTAACTAAAAGATAAGTGACTGTGTTGACATAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAg                                                                                                                                                      >  1:380611/1‑150 (MQ=255)
  aCCTTTTAACTAAAAGATAAGTGACTGTGTTGACATAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGa                                                                                                                                                     >  2:1096723/1‑150 (MQ=255)
    ctattAACTAAAAGATAAGTGACTGTGTTGACATAGTTTTAGCGAGAAATTAATTCTACATAGGAGATCAATATCACATCGCAGAATTACAGTGAGAAAGTGCATAAATTTAGCGGGAAAAGAAATAAGGGAAAGCCAATTTGTCAGACa                                                                                                                                                   >  1:498201/4‑150 (MQ=255)
       ttAACTAAAAGATAAGTGACTGTGTTGACATAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAAt                                                                                                                                                <  1:811584/150‑1 (MQ=255)
       ttAACTAAAAGATAAGTGACTGTGTTGACATAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAAt                                                                                                                                                <  2:1149265/150‑1 (MQ=255)
          aCTAAAAGATAAGTGACTGTGTTGACATAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGt                                                                                                                                             <  1:315926/150‑1 (MQ=255)
           cTAAAAGATAAGTGACTGTGTTGACATAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCTGACAAATTggc                                                                                                                                            <  2:917745/150‑3 (MQ=255)
           cTAAAAGATAAGTGACTGTGTTGACATAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTc                                                                                                                                            >  2:30228/1‑150 (MQ=255)
            tAAAAGATAAGTGACTGTGTTGACATAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCg                                                                                                                                           >  1:520747/1‑150 (MQ=255)
             aaaaGATAAGTGACTGTGTTGACATAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCTGACAAATTggctt                                                                                                                                          <  1:255752/150‑5 (MQ=255)
              aaaGATAAGTGACTGTGTTGACATAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCTGACAAATTggcttt                                                                                                                                         <  1:216349/150‑6 (MQ=255)
                aGATAAGTGACTGTGTTGACATAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATg                                                                                                                                       <  2:259631/150‑1 (MQ=255)
                   tAAGTGACTGTGTTGACATAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGcac                                                                                                                                    <  1:1084834/150‑1 (MQ=255)
                   tAAGTGACTGTGTTGACATAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGcac                                                                                                                                    <  1:188930/150‑1 (MQ=255)
                   tAAGTGACTGTGTTGACATAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGGGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGcac                                                                                                                                    <  1:783325/150‑1 (MQ=255)
                     aGTGACTGTGTTGACATAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCAcag                                                                                                                                  >  1:996126/1‑150 (MQ=255)
                        gACTGTGTTGACATAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCAcagcag                                                                                                                               >  2:969370/1‑150 (MQ=255)
                        gACTGTGTTGACATAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCAcagcag                                                                                                                               <  1:262241/150‑1 (MQ=255)
                           tgtgTTGACATAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGAtt                                                                                                                            <  1:1214404/150‑1 (MQ=255)
                              gtTGACATAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAAt                                                                                                                         >  2:863955/1‑150 (MQ=255)
                               ttGACATAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATc                                                                                                                        <  2:801379/150‑1 (MQ=255)
                                tGACATAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATcc                                                                                                                       <  1:1151214/150‑1 (MQ=255)
                                   cATAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATa                                                                                                                    >  1:46083/1‑150 (MQ=255)
                                     tAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAg                                                                                                                  <  1:716476/150‑1 (MQ=255)
                                     tAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAg                                                                                                                  <  1:1120219/150‑1 (MQ=255)
                                       gTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACaaa                                                                                                                                                 >  1:502134/1‑117 (MQ=255)
                                       gTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACaaa                                                                                                                                                 <  2:502134/117‑1 (MQ=255)
                                       gTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGAt                                                                                                                >  1:502997/1‑150 (MQ=255)
                                         tttAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTa                                                                                                              <  1:5033/150‑1 (MQ=255)
                                              cGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTg                                                                                                         >  1:1443121/1‑150 (MQ=255)
                                               gagaAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTgg                                                                                                        <  1:52940/150‑1 (MQ=255)
                                                 gaAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATa                                                                                                                    >  1:1082791/1‑136 (MQ=255)
                                                   aattaattCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAAt                                                                                                    >  1:633426/1‑150 (MQ=255)
                                                      taattCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCt                                                                                                 >  2:112534/1‑150 (MQ=255)
                                                         ttCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCttgt                                                                                              <  2:394406/150‑1 (MQ=255)
                                                           ctcCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCttgttg                                                                                            <  2:1043784/150‑1 (MQ=255)
                                                            tcCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTtgttgt                                                                                           >  1:315940/1‑150 (MQ=255)
                                                              cATAGGAGAGCAATATCACATCGCAGAATTAAAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAAATTGTCAGACAAATTTTCGAATGAACAGCAGATTAATCCATAAGATTACAATGGAAATCCAtgatgtaa                                                                                         >  1:987184/1‑148 (MQ=255)
                                                                tAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCttt                                                                                       >  2:1537301/1‑150 (MQ=255)
                                                                 aGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTg                                                                                      >  2:44705/1‑150 (MQ=255)
                                                                 aGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTg                                                                                      <  2:553289/150‑1 (MQ=255)
                                                                  ggAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTgg                                                                                     <  1:259997/150‑1 (MQ=255)
                                                                  ggAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTgg                                                                                     >  2:1434420/1‑150 (MQ=255)
                                                                    agagCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGt                                                                                    >  1:931038/1‑149 (MQ=255)
                                                                    agagCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTa                                                                                   <  1:366502/150‑1 (MQ=255)
                                                                         aaTATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCAt                                                                              >  2:169902/1‑150 (MQ=255)
                                                                         aaTATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCAt                                                                              >  1:652421/1‑150 (MQ=255)
                                                                             tCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCgg                                                                          >  1:147881/1‑150 (MQ=255)
                                                                              cacaTCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCggg                                                                         >  2:897020/1‑150 (MQ=255)
                                                                              cacaTCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCggg                                                                         >  1:376556/1‑150 (MQ=255)
                                                                               acaTCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGa                                                                        <  1:1340855/150‑1 (MQ=255)
                                                                                  tCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGt                                                                     >  1:1013001/1‑150 (MQ=255)
                                                                                    gCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGAGGGATGtct                                                                   >  1:688922/1‑150 (MQ=255)
                                                                                       gAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGtcttct                                                                <  2:537089/150‑1 (MQ=255)
                                                                                       caaTTACACTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGTAAGCCAATTTGTCAGACAAATTGTCGAATGCAGAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGtcttct                                                                <  1:855985/149‑1 (MQ=255)
                                                                                        aaTTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTcttctt                                                               >  1:1514580/1‑150 (MQ=255)
                                                                                         aTTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAAACCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCttt                                                              >  2:771513/1‑150 (MQ=255)
                                                                                            aCAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTa                                                           <  2:582052/150‑1 (MQ=255)
                                                                                               gTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAAcc                                                        >  2:1416715/1‑150 (MQ=255)
                                                                                               gTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAAcc                                                        >  1:76855/1‑150 (MQ=255)
                                                                                               gTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGCCTTTGGTACCCATGCGGGATGTCTTCTTTTTAAcc                                                        <  2:380611/150‑1 (MQ=255)
                                                                                                   gaACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTc                                                    <  2:633426/150‑1 (MQ=255)
                                                                                                   gaACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTc                                                    <  1:439717/150‑1 (MQ=255)
                                                                                                   gaACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTc                                                    <  1:1149437/150‑1 (MQ=255)
                                                                                                      cGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAAt                                                 <  1:1442688/150‑1 (MQ=255)
                                                                                                       gTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATa                                                <  1:1434420/150‑1 (MQ=255)
                                                                                                         gCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAgg                                              >  2:1351993/1‑150 (MQ=255)
                                                                                                          cATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGc                                             >  1:917739/1‑150 (MQ=255)
                                                                                                             aaaTTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGc                                          >  2:2770/1‑150 (MQ=255)
                                                                                                               aTTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCAt                                        <  2:986929/150‑1 (MQ=255)
                                                                                                                tttAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCAtt                                       >  2:330695/1‑150 (MQ=255)
                                                                                                                         aaaaGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGCGt                              <  2:502997/150‑1 (MQ=255)
                                                                                                                           aaGACTTAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGTAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGCGTTg                            <  2:1409881/150‑1 (MQ=255)
                                                                                                                           aaGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGCGTTg                            <  2:1113030/150‑1 (MQ=255)
                                                                                                                             gACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGCGTTGAg                          >  1:364289/1‑150 (MQ=255)
                                                                                                                             caCATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGCGTTGAg                          <  2:1074649/149‑1 (MQ=255)
                                                                                                                                 tAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGCGTTGAGtttt                      >  1:469422/1‑150 (MQ=255)
                                                                                                                                 tAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGCGTTGAGtttt                      <  2:992009/150‑1 (MQ=255)
                                                                                                                                   aGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGCGTTGAGTTTTTg                    >  2:498913/1‑150 (MQ=255)
                                                                                                                                   aGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCACAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGCGTTGAGTTTTTg                    >  1:1258270/1‑150 (MQ=255)
                                                                                                                                    gggAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGCGTTGAGTTTTTGa                   >  2:263756/1‑150 (MQ=255)
                                                                                                                                    gggAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGCGTTGAGTTTTTGa                   >  2:631089/1‑150 (MQ=255)
                                                                                                                                     cattcgacAATTTGTCTGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGCGTTGAGTTTTTGaa                  >  1:121266/8‑150 (MQ=255)
                                                                                                                                      gAAGGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGCGTTGAGTTTTTGaaa                 >  1:1272760/1‑150 (MQ=255)
                                                                                                                                          gacAATTTGTCTGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGCGTTGAGTTTTTGAAATGgt             >  1:814158/3‑150 (MQ=255)
                                                                                                                                          gacAATTTGTCTGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGCGTTGAGTTTTTGAAATGgt             >  1:384927/3‑150 (MQ=255)
                                                                                                                                          gCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGCGTTGAGTTTTTGAAATGgt             >  1:448094/1‑150 (MQ=255)
                                                                                                                                          gCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGCGTTGAGTTTTTGAAATGgt             >  1:200007/1‑150 (MQ=255)
                                                                                                                                           ccAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGCGTTGAGTTTTTGAAATGgtg            <  2:484035/150‑1 (MQ=255)
                                                                                                                                            cAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGCGTTGAGTTTTTGAAATGgtgt           <  1:969370/150‑1 (MQ=255)
                                                                                                                                            cAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCAGGAAATCCTAGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTAAACCAGTCAATAGGCCGCATTACCGCGCGTTGAGTATTTGAAATGgtgt           >  2:615287/1‑150 (MQ=255)
                                                                                                                                             aaTTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGCGTTGAGTTTTTGAAATGGTGta          >  1:695675/1‑150 (MQ=255)
                                                                                                                                             aaTTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGCGTTGAGTTTTTGAAATGGTGta          >  2:924971/1‑150 (MQ=255)
                                                                                                                                                ttGTCAGACAAATTGTCGAATACACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGCGTTGAGTTTTTGAAATGGTGtaata       >  1:1047677/1‑150 (MQ=255)
                                                                                                                                                  gTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGCGTTGAGTTTTTGAAATGGTGTAATAAc     >  1:44830/1‑150 (MQ=255)
                                                                                                                                                    cAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAGTGCTTGATGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGCGTTGAGTTTTTGAAATGGTGTAATAACCg   >  2:1509106/1‑150 (MQ=255)
                                                                                                                                                    cAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGCGTTGAGTTTTTGAAATGGTGTAATAACCg   >  1:1104068/1‑150 (MQ=255)
                                                                                                                                                    cAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATACTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGCGTTGAGTTTTTGAAATGGTGTAATAACCg   >  2:1503910/1‑150 (MQ=255)
                                                                                                                                                     aGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGCGTTGAGTTTTTGAAATGGTGTAATAACCGc  >  2:96293/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTACCTTTTAACTAAAAGATAAGTGACTGTGTTGACATAGTTTTAGCGAGAAATTAATTCTCCATAGGAGAGCAATATCACATCGCAGAATTACAGTGAGAACGTGCATAAATTTAGCGGGAAAAGACATAAGGGAAAGCCAATTTGTCAGACAAATTGTCGAATGCACAGCAGATTAATCCATAAGATTAGCCTGGAAATCCTTGTTGTCTTTGGTACCCATGCGGGATGTCTTCTTTTTAACCAGTCAATAGGCCGCATTACCTGGCGTTGAGTTTTTGAAATGGTGTAATAACCGC  >  NC_000913/1986207‑1986505

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