Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,194,790 A→G 5.7% noncoding (76/145 nt)
I9V (ATT→GTT) 		 sibD ←
ibsD →		 small RNA SibD
putative toxic peptide IbsD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,194,7900AG 5.7% 266.3 / 3.7 86noncoding (76/145 nt)
I9V (ATT→GTT) 		sibD
ibsD		small RNA SibD
putative toxic peptide IbsD
Reads supporting (aligned to +/- strand):  ref base A (38/43);  new base G (3/2);  total (41/45)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAAAACTGCATCAGAGAAAGCATTGCGTGGGGTCTGGCGCTCATCAACAATCGGGGAGTCAGCAGGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCTTATTAACAGCCAATCAGAGGGGGGAGAAATCCTCCCTCACCCTTGTTCCTTTACTTTAGGTTGAAAAAAAATCAGCGTCAAGAAGGTTTCATTAACACATATCCCTGGCCCAAAAAGCC  >  NC_000913/3194643‑3194938
                                                                                                                                                   |                                                                                                                                                    
gCAAAACTGCATCAGAGAAAGCATTGCGTGGGGTCTGGCGCTCATCAACAATCGGGGAGTCAGCAGGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGCAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGAtt                                                                                                                                                    >  1:528764/1‑150 (MQ=255)
 cAAAACTGCATCAGTGAAAGGATTGGGTGGGGTCTGGCGCTCATCAAGAATCGGGGAGTCATCAGGGGGCTGAAACGGGAAAGCGCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATtg                                                                                                                                                   <  2:306219/150‑1 (MQ=18)
 cAAAACTGCATCAGAGAAAGCATTGCGTGGGGTCTGGCGCTCATCAACAATCGGGGAGTCAGCAGGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATtg                                                                                                                                                   >  1:574448/1‑150 (MQ=255)
 cAAAACTGCATCAGAGAAAGCATTGCGTGGGGTCTGGCGCTCATCAACAATCGGGGAGGCAGCAGGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATtg                                                                                                                                                   <  1:1572576/150‑1 (MQ=255)
 cAAAACTGCATCACCCAAAGCCTTGCGTGGGGTCTGGCCCTCATCAACAATCGGGGAGTCAGCAGGGGGCTGAAACGGGAAACCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATtg                                                                                                                                                   <  1:293934/150‑1 (MQ=18)
 cAAAAATGCATCAGAGAAAGCATTGCGTGGGGTCTGGCGCTCATCAACAATCGGGGAGTCAGCAGGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATtg                                                                                                                                                   <  1:66401/150‑1 (MQ=255)
    aaCTGCATCAGAGAAAGCATTGCGTGGCGTCTGGCGCTCATCAACAATCGGGGAGTCAGCAGGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTgtgt                                                                                                                                                <  1:102974/150‑1 (MQ=255)
        gCATCAGAGAAAGCATTGCGTGGGGTCTGGCGCTCATCAACAATCGGGGAGTCAGCAGGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTgttgtt                                                                                                                                            >  2:1331848/1‑150 (MQ=255)
         cATCAGAGAAAGCATTGCGTGGGGTCTGGCGCTCATCAACAATCGGGGAGTCAGCAGGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTa                                                                                                                                           <  2:1327627/150‑1 (MQ=255)
         cATCAGAGAAAGCATTGCGTGGGGTCTGGCGCTCATCAACAATCGGGGAGTCAGCAGGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTa                                                                                                                                           >  1:1272257/1‑150 (MQ=255)
          aTCAGAGAAAGCATTGCGTGGGGTCTGGCGCTCATCAACAATCGGGGAGTCAGCAGGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTAc                                                                                                                                          >  2:1463538/1‑150 (MQ=255)
            cAGAGAAAGCATTGCGTGGGGTCTGGCGCTCATCAACAATCGGGGAGTCAGCAGGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTc                                                                                                                                        >  1:332598/1‑150 (MQ=255)
              gagaAAGCATTGCGTGGGGTCTGGCGCTCATCAACAATCGGGGAGTCAGCAGGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGt                                                                                                                                      >  1:932020/1‑150 (MQ=255)
                gaAAGCATTGCGTGGGGTCTGGCGCTCATCAACAATCGGGGAGTCAGCAGGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTaa                                                                                                                                    <  2:1520903/150‑1 (MQ=255)
                gaAAGCATTGCGTGGGGTCTGGCGCTCATCAACAATCGGGGAGTCAGCAGGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGAGTGTGTTGTTACTCGTat                                                                                                                                    >  2:1037996/1‑149 (MQ=255)
                    gCATTGCGTGGGGTCTGGCGCTCATCAACAATCGGGGAGTCAGCAGGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGttt                                                                                                                                >  1:1282724/1‑150 (MQ=255)
                       ttGCGTGGGGTCTGGCGCTCATCAACAATCGGGGAGTCAGCAGGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGc                                                                                                                             <  2:197985/150‑1 (MQ=255)
                       ttGCGTGGGGTCTGGAGCTCATCAACAATCGGGGAGTCAGCAGGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGc                                                                                                                             <  2:1050637/150‑1 (MQ=255)
                        tGCGTGGGGTCTGGCGCTCATCAACAATCGGGGAGTCAGCAGGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGTTTcccg                                                                                                                            >  1:1202292/1‑147 (MQ=255)
                        tGCGTGGGGTCTGGCGCTCATCAACAATCGGGGAGTCAGCAGGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCa                                                                                                                            <  2:987871/150‑1 (MQ=255)
                         gCGTGGGGTCTGGCGCTCATCAACAATCGGGGAGTCAGCAGGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAg                                                                                                                           >  1:1474940/1‑150 (MQ=255)
                          cGTGGGGTCTGGCGCTCATCAACAATCGGGGAGTCAGCAGGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGc                                                                                                                          >  2:1067016/1‑150 (MQ=255)
                           gTGGGGTCTGGCGCTCATCAACAATCGGGGAGTCAGCAGGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCt                                                                                                                         <  1:1525000/150‑1 (MQ=255)
                                gTCTGGCGCTCATCAACAATCGGGGAGTCAGCAGGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCttatta                                                                                                                    <  2:567353/150‑1 (MQ=255)
                                gTCTGGCGCTCATCAACAATCGGGGAGTCAGCAGGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCttatta                                                                                                                    <  1:1330073/150‑1 (MQ=255)
                                 tCTGGCGCTCATCAACAATCGGGGAGTCAGCAGGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGTGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCTTATTaa                                                                                                                   <  2:1354064/150‑1 (MQ=255)
                                    ggCGCTCATCAACAATCGGGGAGTCAGTAGGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCTTATTAACAg                                                                                                                >  2:1114230/1‑150 (MQ=255)
                                          catcaACAATCGGGGAGTCAGCAGGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCTTATTAACAGCCAATc                                                                                                          >  2:37871/1‑150 (MQ=255)
                                           atcaACAATCGGGGAGTCAGCAGGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCTTATTAACAGCCAATCa                                                                                                         <  1:212180/150‑1 (MQ=255)
                                             caacaaTCGGGGAGTCAGCAGGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCTTATTAACAGCCAATCaga                                                                                                       >  1:35153/1‑150 (MQ=255)
                                                  aTCGGGGAGTCAGCAGGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCTTATTAACAGCCAATCAGAggggg                                                                                                  <  1:1298434/150‑1 (MQ=255)
                                                     ggggAGTCAGCAGGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCTTATTAACAGCCAATCAGAGGGGGgag                                                                                               <  2:202111/150‑1 (MQ=255)
                                                     ggggAGTCAGCAGGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCTTATTAACAGCCAATCAGAGGGGGgag                                                                                               <  1:788985/150‑1 (MQ=255)
                                                            cagcagGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCTTATTAACAGCCAATCAGAGCGGGGAGAAAtcct                                                                                        >  2:435593/1‑150 (MQ=255)
                                                              gcagGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCTTATTAACAGCCAATCAGAGGGGGGAGAAAtcctcc                                                                                      >  2:1513646/1‑150 (MQ=255)
                                                               cagGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCTTATTAACAGCCAATCAGAGGGGGGAGAAATcctccc                                                                                     <  2:981225/150‑1 (MQ=255)
                                                               cagGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCTTATTAACAGCCAATCAGAGGGGGGAGAAATcctccc                                                                                     <  2:1386083/150‑1 (MQ=255)
                                                               cagGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCTTATTAACAGCCAATCAGAGGGGGGAGAAATcctccc                                                                                     <  2:1282724/150‑1 (MQ=255)
                                                                agGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCTTATTAACAGCCAATCAGAGGGGGGAGAAATcctccct                                                                                    <  2:25523/150‑1 (MQ=255)
                                                                     gCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCTTATTAACAGCCAATCAGAGGGGGGAGAAATCCTCCCTCAccc                                                                               <  1:1513646/150‑1 (MQ=255)
                                                                      cTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCTTATTAACAGCCAATCAGAGGGGGGAGAAATCCTCCCTCACCCt                                                                              <  1:245448/150‑1 (MQ=255)
                                                                      cTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCTTATTAACAGCCAATCAGAGGGGGGAGAAATCCTCCCTCACCCt                                                                              >  1:1249458/1‑150 (MQ=255)
                                                                       tGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCTTATTAACAGCCAATCAGAGGGGGGAGAAATCCTCCCTCACCCtt                                                                             >  2:1092594/1‑150 (MQ=255)
                                                                        gAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCTTATTAACAGCCAATCAGAGGGGGGAGAAATCCTCCCTCACCCTTg                                                                            <  2:1272257/150‑1 (MQ=255)
                                                                         aaaCGGGAAATCCACTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCTTATTAACAGCCAATCAGAGGGGGGAGAAATCCTCCCTCACCCTTGt                                                                           >  2:856611/1‑150 (MQ=255)
                                                                           aCGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCTTATTAACAGCCAATCAGAGGGGGGAGAAATCCTCCCTCACCCTTGTTc                                                                         >  2:835077/1‑150 (MQ=255)
                                                                            cGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCTTATTAACAGCCAATCAGAGGGGGGAGAAATCCTCCCTCACCCTTGTTcc                                                                        <  1:37871/150‑1 (MQ=255)
                                                                             gggAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCTTATTAACAGCCAATCAGAGGGGGGAGAAATCCTCCCTCACCCTTGTTc                                                                         >  2:635447/1‑148 (MQ=255)
                                                                             gggAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCTTATTAACAGCCAATCAGAGGGGGGAGAAATCCTCCCTCACCCTTGTTc                                                                         <  1:635447/148‑1 (MQ=255)
                                                                                aaaGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCTTATTAACAGCCAATCAGAGGGGGGAGAAATCCTCCCTCACCCTTGTTCcttta                                                                    <  1:866384/150‑1 (MQ=255)
                                                                                aaaGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCTTATTAACAGCCAATCAGAGGGGGGAGAAATCCTCCCTCACCCTTGTTCcttta                                                                    <  1:509544/150‑1 (MQ=255)
                                                                                 aaGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCTTATTAACAGCCAATCAGAGGGGGGAGAAATCCTCCCTCACCCTTGTTCctttac                                                                   <  2:1474940/150‑1 (MQ=255)
                                                                                    ccccTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCTTATTAACAGCCAATCAGAGGGGGGAGAAATCCTCCCTCACCCTTGTTCctttacttt                                                                <  1:1331848/150‑1 (MQ=255)
                                                                                      ccTCCCGAAGAAGGGGGCCTTACAGAAGCAAAGGGTTATGCTGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCTTATTAACAGCCAATCAGAGGGGGGAGAAATCCTCCCTCACCCTTGTTCCTTTACTTTAg                                                              <  2:950318/150‑1 (MQ=255)
                                                                                       cTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCTTATTAACAGCCAATCAGAGGGGGGAGAAATCCTCCCTCACCCTTGTTCCTTTACTTTAgg                                                             <  1:1463538/150‑1 (MQ=255)
                                                                                        gcccGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCTTATTAACAGCCAATCAGAGGGGGGAGAAATCCTCCCTCACCCTTGTTCCTTTACTTTAGGt                                                            <  1:188235/149‑1 (MQ=255)
                                                                                            gaagaaGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCTTATTAACAGCCAATCAGAGGGGGGAGAAATCCTCCCTCACCCTTGTTCCTTTACTTTAGGtt                                                           <  1:1276891/147‑1 (MQ=255)
                                                                                            gaagaaGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCTTATTAACAGCCAATCAGAGGGGGGAGAAATCCTCCCTCACCCTTGTTCCTTTACTTTAGGtt                                                           >  2:1276891/1‑147 (MQ=255)
                                                                                               gaagGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCTTATTAACAGCCAATCAGAGGGGGGAGAAATCCTCCCTCACCCTTGTTCCTTTACTTTAGGTTGaaaaa                                                     >  2:995357/1‑150 (MQ=255)
                                                                                                 agGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCTTATTAACAGCCAATCAGAGGGGGGAGAAATCCTCCCTCACCCTTGTTCCTTTACTTTAGGTTGaaaaaaa                                                   >  1:217124/1‑150 (MQ=255)
                                                                                                    gggCCTTACAGAAGGAAAGGGTTATGATGAAGCACGTCATCATACTGGTTGtgtt                                                                                                                                               >  2:1084892/1‑55 (MQ=14)
                                                                                                    gggCCTTACAGAAGGAAAGGGTTATGATGAAGCACGTCATCATACTGGTTGtgtt                                                                                                                                               >  2:1079310/1‑55 (MQ=14)
                                                                                                    gggCCTTACAGAAGGAAAGGGTTATGATGAAGCACGTCATCATACTGGTTGtgtt                                                                                                                                               <  1:1084892/55‑1 (MQ=14)
                                                                                                    gggCCTTACAGAAGGAAAGGGTTATGATGAAGCACGTCATCATACTGGTTGtgtt                                                                                                                                               <  1:1079310/55‑1 (MQ=14)
                                                                                                       ccTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCTTATTAACAGCCAATCAGAGGGGGGAGAAATCCTCCCTAACCCTTGTTCCTTTACATTAGGTTGAAAAAAAATCAGc                                             >  1:1532888/1‑150 (MQ=255)
                                                                                                        cTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCTTATTAACAGCCAATCAGAGGGGGGAGAAATCCTCCCTCACCCTTGTTCCTTTACTTTAGGTTGAAAAAAAATCAGCg                                            <  2:236138/150‑1 (MQ=255)
                                                                                                        cTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCTTATTAACAGCCAATCAGAGGGGGGAGAAATCCTCCCTCACCCTTGTTCCTTTACTTTAGGTTGAAAAAAAATCAGCg                                            >  2:337041/1‑150 (MQ=255)
                                                                                                          tACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCTTATTAACAGCCAATCAGAGGGGGGAGAAATCCTCCCTCACCCTTGTTCCTTTACTTTAGGTTGAAAAAAAATCACCGTc                                          >  1:706399/1‑150 (MQ=255)
                                                                                                           aCAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCTTATTAACAGCCAATCAGAGGGGGGAGAAATCCTCCCTCACCCTTGTTCCTTTACTTTAGGTTGAAAAAAAATCAGCGTCa                                         <  2:332598/150‑1 (MQ=255)
                                                                                                            cAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCTTATTAACAGCCAATCAGAGGGGGGAGAAATCCTCCCTCACCCTTGTTCCTTTACTTTAGGTTGAAAAAAAATCAGCGTCaa                                        >  1:890054/1‑150 (MQ=255)
                                                                                                            cAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCTTATTAACAGCCAATCAGAGGGGGGAGAAATCCTCCCTCACCCTTGTTCCTTTACTTTAGGTTGAAAAAAAATCAGCGTCaa                                        >  2:542444/1‑150 (MQ=255)
                                                                                                            cAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCTTATTAACAGCCAATCAGAGGGGGGAGAAATCCTCCCTCACCCTTGTTCCTTTACTTTAGGTTGAAAAAAAATCAGCGTCaa                                        >  2:342071/1‑150 (MQ=255)
                                                                                                              gAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCTTATTAACAGCCAATCAGAGGGGGGAGAAATCCTCCCTCACCCTTGTTCCTTTACTTTAGGTTGAAAAAAAATCAGCGTCaaga                                      >  1:894954/1‑150 (MQ=255)
                                                                                                                      gggTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCTTATTAACAGCCAATCAGAGGGGGGAGAAATCCTCCCTCACCCTTGTTCCTTTACTTTAGGTTGAAAAAAAATCAGCGTCAAGAAGGTTTCa                              >  1:928349/1‑150 (MQ=255)
                                                                                                                         ttATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCTTATTAACAGCCAATCAGAGGGGGGAGAAATCCTCCCTCACCCTTGTTCCTTTACTTTAGGTTGAAAAAAAATCAGCGTCAAGAAGGTTTCATTa                           >  1:1302488/1‑150 (MQ=255)
                                                                                                                          tATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCTTATTAACAGCCAATCAGAGGGGGGAGAAATCCTCCCTCACCCTTGTTCCTTTACTTTAGGTTGAAAAAAAATCAGCGTCAAGAAGGTTTCATTaa                          >  2:122111/1‑150 (MQ=255)
                                                                                                                           atgatgAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCTTATTAACAGCCAATCAGAGGGGGGAGAAATCCTCCCTCACCCTTGTTCCTTTACTTTAGGTTGAAAAAAAATCAGCGTCAAGAAGGTTTCATTAac                         <  1:93667/150‑1 (MQ=255)
                                                                                                                           atgatgAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCTTATTAACAGCCAATCAGAGGGGGGAGAAATCCTCCCTCACCCTTGTTCCTTTACTTTAGGTTGAAAAAAAATCAGCGTCAAGAAGGTTTCATTAac                         <  1:835077/150‑1 (MQ=255)
                                                                                                                             gatgaAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCTTATTAACAGCCAATCAGAGGGGGGAGAAATCCTCCCTCCCCCTTGTTCCTTTACTTTAGGTTGAAAAAAAATCAGCGTCAAGAAGGTTTCATTAacac                       >  2:714587/1‑150 (MQ=255)
                                                                                                                                      cGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCTTATTAACAGCCAATCAGAGGGGGGAGAAATCCTCCCTCACCCTTGTTCCTTTACTTTAGGTTGAAAAAAAATCAGCGTCAAGAAGGTTTCATTAACACATATCCCTg              >  1:814346/1‑150 (MQ=255)
                                                                                                                                      cGTCATCATACTGATTGTGGTGTTACTCGTAATTTTCGCAGCTTATTAACAGCCAATCAGAGGGGGGAGAAATCCTCCCTCTCGCTTGTTCCTTTACTTTAGGTTGAAAAAAAATCAGCGTCAAGAAGGTTTCATTAACACATTTccgtt              <  2:979225/150‑4 (MQ=25)
                                                                                                                                           tcatACTGATTGTTTTGTTACTCGTAAGTTTCGCAGCTTATTAACAGCCAATCAGAGGGTGGAGAAATCCTCCCTCACCCTTGTTCCTTTACTTTAGGTTGAAAAAAAATCAGCGTCAAGAAGGTTTCATTAACACATATCCCTGGCCCa         <  1:1038473/150‑1 (MQ=255)
                                                                                                                                           tcatACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCTTATTAACAGCCAATCAGAGGGGGGAGAAATCCTCCCTCACCCTTGTTCCTTTACTTTAGGTTGAAAAAAAATCAGCGTCAAGAAGGTTTCATTAACACATATCCCTGGCCCa         <  1:856611/150‑1 (MQ=255)
                                                                                                                                            catACTCATTGTGTTGTTACTCGTAAGTTTCGCAGCTTATTAACAGCCAATCAGAGGGGGGAGAAATCCTCCCTCACCCTTGTTCCTTTACTTTAGGTTTAAAAAAAATCAGCGTCAATAAGGATTCATTAACACATATACCTGGCACaa        >  2:1038473/1‑150 (MQ=25)
                                                                                                                                                 tGATTGTGTTTTTAGTCGTAATTTTCGCAGCTTTTTAACAGCCAATCTGTGGGGGGAGAAATCCTCCCTCACCCTTGTTCCTTTACTTTAGTTTGAAAAAAAATCAGCTTCAAGAAGGTTTAATTAACACATATCGCTGGCCCAAAAAGc   <  2:528764/150‑1 (MQ=255)
                                                                                                                                                  gATTGTGTTGTTACTCGTAAGCTTCGCAGCTTATTAACAGCCAATCAGAGGGGGGAGAAATCCTCCCTCACCCTTGTTCCTTTACTTTAGGTTGAAAAAAAATCAGCGTCAAGAAGGTTTCATTAACACATATCCCTGGCCCAAAAAGcc  >  1:1255772/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GCAAAACTGCATCAGAGAAAGCATTGCGTGGGGTCTGGCGCTCATCAACAATCGGGGAGTCAGCAGGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGATTGTGTTGTTACTCGTAAGTTTCGCAGCTTATTAACAGCCAATCAGAGGGGGGAGAAATCCTCCCTCACCCTTGTTCCTTTACTTTAGGTTGAAAAAAAATCAGCGTCAAGAAGGTTTCATTAACACATATCCCTGGCCCAAAAAGCC  >  NC_000913/3194643‑3194938

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