Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,286,311 C→G 18.5% intergenic (+24/+79) proL → / ← yejO tRNA‑Pro/pseudogene, autotransporter outer membrane homology;putative transport; Not classified; putative ATP‑binding component of a transport system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,286,3110CG18.5% 333.8 / 65.3 135intergenic (+24/+79)proL/yejOtRNA‑Pro/pseudogene, autotransporter outer membrane homology;putative transport; Not classified; putative ATP‑binding component of a transport system
Reads supporting (aligned to +/- strand):  ref base C (61/49);  new base G (12/13);  total (73/62)
Fisher's exact test for biased strand distribution p-value = 5.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.48e-01

Only 100 of 137 total aligned reads displayed.

GGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACC  >  NC_000913/2286212‑2286410
                                                                                                   |                                                                                                   
ggCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAAccc                                                                                                     >  1:263086/1‑100 (MQ=255)
   aCGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTa                                                                                                  >  2:1646700/1‑100 (MQ=255)
        gcgcAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACggtt                                                                                             <  2:1865975/100‑1 (MQ=255)
         cgcAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACggttg                                                                                            >  2:2002168/1‑100 (MQ=255)
         cgcAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGATTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACggtt                                                                                             <  2:3782242/99‑1 (MQ=255)
          gcAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACggttgg                                                                                           >  1:367277/1‑100 (MQ=255)
           cAGCCTGGTATAGCACCGTCATGGGGTGTCGGGGGTCGGAGATTCAAATCCTCTCGTGCCGACCAAAAAAACAAAAAAAAAACCAACCCTTACggttg                                                                                            <  2:3152380/98‑1 (MQ=255)
             gCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCGTAAGGGTTGGttt                                                                                        <  1:3351706/100‑1 (MQ=39)
              ccTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCGTAAGGGTTGGtttt                                                                                       <  1:2946/100‑1 (MQ=255)
              ccTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGtttt                                                                                       <  1:2439892/100‑1 (MQ=255)
               cTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGttttt                                                                                      >  1:120312/1‑100 (MQ=255)
                tGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGACAAAACCAACCCTTACGGTTGGtttttt                                                                                     >  2:1630811/1‑100 (MQ=255)
                tGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGttttt                                                                                      <  1:2681363/99‑1 (MQ=255)
                 ggTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGttttttt                                                                                    <  1:3491525/100‑1 (MQ=255)
                  gTATCGCACCGTCATGGAGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCTACCAAAAATCCCAAGACAACACCAACCGTAAGGGTTGGtttttcct                                                                                   <  2:2794619/100‑4 (MQ=255)
                  gTAGCGCACCGTTATGGGGTGACGGGGGTCGGAGGTTGAAATCCTCTCGTGACGACCAAAAATCCCAAGAAAAAACCAACCGTAAGGGTTGGttttttct                                                                                   <  2:619861/100‑3 (MQ=37)
                  gTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCGTAAGGGTTGGttttttct                                                                                   <  2:212203/100‑3 (MQ=38)
                  gTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTta                                                                                   <  2:3648557/100‑1 (MQ=255)
                    aGCGCACCGTCATGGGGTGTCGGGGGTCGGCGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTtata                                                                                 >  2:2750472/1‑100 (MQ=255)
                    aGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAACTCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTtata                                                                                 >  2:1435837/1‑100 (MQ=255)
                    aGCGCACAGTCAGGGAGTATCGGAGCGCGCAGGTTCAAATCCTCTCGTGTCGACCAACTATCCCAAGAAAAAACCAAGCGCAATGGTTGGgtttttct                                                                                   <  2:659483/98‑3 (MQ=255)
                     gcgcACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCGTAAGGGTTGGttttttcttgg                                                                                <  1:1053302/100‑6 (MQ=37)
                      cgcACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATc                                                                               >  1:1841973/1‑100 (MQ=255)
                       gcACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTGTTATATc                                                                               >  2:366912/1‑99 (MQ=255)
                        cACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTg                                                                             >  2:3296991/1‑100 (MQ=255)
                        cACCCTCATGGGGTATCGGCGGTCCGAAGATCAAATCCGCACATGCCGACTAAAAATCCCAAGAAAAAACCAACCGTAAGGGTTGGttttttcttgg                                                                                <  2:2321511/97‑6 (MQ=12)
                            gTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCaat                                                                         <  1:1617334/100‑1 (MQ=255)
                                tGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCaattaat                                                                     >  2:2667287/1‑100 (MQ=255)
                                 ggggTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCaattaatt                                                                    <  1:727149/100‑1 (MQ=255)
                                  gggTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTc                                                                   >  1:2624993/1‑100 (MQ=255)
                                  gggTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTc                                                                   <  2:1605915/100‑1 (MQ=255)
                                     tgtCGGGCGGCGAAGGTACAAATTCTCTCGTGCCGACCAAAAATCCCAAGAATAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCGATTAATTCGa                                                                 <  2:1231843/99‑1 (MQ=255)
                                      gtCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATa                                                               >  1:2699858/1‑100 (MQ=255)
                                      gtCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTACAATTAATTCGATa                                                               >  2:879912/1‑100 (MQ=255)
                                        cGGGGGTCGGAGGTTCAAATCATCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATaaa                                                             >  2:1277009/1‑100 (MQ=255)
                                         gggggTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAAc                                                            <  2:1841963/100‑1 (MQ=255)
                                          ggggTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACa                                                           <  1:1787469/100‑1 (MQ=255)
                                            ggTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACa                                                           <  2:468754/98‑1 (MQ=255)
                                                ggAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGt                                                     >  2:2996086/1‑100 (MQ=255)
                                                 gAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTg                                                    >  1:572910/1‑100 (MQ=255)
                                                   ggTTCAAATCCTCTCGTGCCGACCAAAAATCGCAAGAAAAAACCAACACTTACGGTTGGTTTTTTTATATCTGCCATTAATTCGATAAACCGACCGTGac                                                  >  2:1690703/1‑100 (MQ=255)
                                                   ggTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGac                                                  >  1:3173519/1‑100 (MQ=255)
                                                   ggTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGac                                                  >  1:3115654/1‑100 (MQ=255)
                                                   ggTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGac                                                  <  1:602667/100‑1 (MQ=255)
                                                    tccgatctCCTCTCGTGCCGACCAAAAATCCCGAGAAAAAACCGACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGaca                                                 <  2:2447257/93‑1 (MQ=255)
                                                    gTTCAAATCCTCTCGTGCAGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGaca                                                 <  2:1011525/100‑1 (MQ=255)
                                                    gTTCAAATCCTCCCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGaca                                                 >  2:73207/1‑100 (MQ=255)
                                                     ttCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGacac                                                >  2:2229309/1‑100 (MQ=255)
                                                      tCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGAcaca                                               <  2:2123489/100‑1 (MQ=255)
                                                       cAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACAt                                              >  1:2152005/1‑100 (MQ=255)
                                                        aaaTCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAGTTAATTCGATAAACAGACCGTGACACATc                                             >  1:3541418/1‑100 (MQ=255)
                                                            ccTCTCGTGCCGACCAAAAATCCCGAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGacaagatcgga                                         >  1:2447273/1‑92 (MQ=255)
                                                            ccTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGacaagatcgga                                         >  1:2866452/1‑92 (MQ=255)
                                                             ctctCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGc                                        >  1:1226901/1‑100 (MQ=255)
                                                             ctctCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGc                                        <  1:2389358/100‑1 (MQ=255)
                                                               ctcGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCt                                      <  1:423356/100‑1 (MQ=255)
                                                                 cGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCt                                      >  2:3574546/1‑98 (MQ=255)
                                                                  gTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGtt                                   <  1:1542506/100‑1 (MQ=255)
                                                                  gTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGtt                                   >  1:1989215/1‑100 (MQ=255)
                                                                    gCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTa                                 >  1:966768/1‑100 (MQ=255)
                                                                    gCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTa                                 >  1:1003920/1‑100 (MQ=255)
                                                                    gCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTGTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACCGCGTGTTTa                                 >  2:2478340/1‑100 (MQ=255)
                                                                      cGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTAtt                               >  2:2493386/1‑100 (MQ=255)
                                                                       gACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTAttt                              >  2:1748812/1‑100 (MQ=255)
                                                                       gACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTAttt                              >  1:2395389/1‑100 (MQ=255)
                                                                         ccAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTAtttt                             >  1:2408091/1‑99 (MQ=255)
                                                                         ccAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTc                            >  2:903977/1‑100 (MQ=255)
                                                                         ccAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTGTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTc                            <  2:2685187/100‑1 (MQ=255)
                                                                          cAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCt                           >  2:1518245/1‑100 (MQ=255)
                                                                          cAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCt                           >  1:3706116/1‑100 (MQ=255)
                                                                          cAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGAATTTTTATATCTGCAATTAATTCGATAAACAGACAGTGACACATCACAGCCTGTTTAtttt                             >  2:3379552/1‑98 (MQ=255)
                                                                             aaaTCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGtt                        <  1:916887/100‑1 (MQ=255)
                                                                              aaTCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTa                       >  1:184009/1‑100 (MQ=255)
                                                                              aaTCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTa                       >  1:2908810/1‑100 (MQ=255)
                                                                               aTCCCAAGAAAAAACCAACCCTTACGGGTTGGTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTa                       <  2:3541391/99‑1 (MQ=255)
                                                                                tctcAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATc                     <  1:2611560/97‑1 (MQ=255)
                                                                                 gatataaaaaaaCCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCa                    >  1:1755322/7‑100 (MQ=255)
                                                                                  ccAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAg                   >  1:1036651/1‑100 (MQ=255)
                                                                                  ccAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAg                   <  1:1323159/100‑1 (MQ=255)
                                                                                  ccAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAg                   >  1:1905897/1‑100 (MQ=255)
                                                                                  ccAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAg                   >  1:2546267/1‑100 (MQ=255)
                                                                                   tataaaaaaaCCAACCGTAAGGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGa                  >  1:585807/5‑100 (MQ=38)
                                                                                   cAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGa                  >  2:2768900/1‑100 (MQ=255)
                                                                                   cAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGa                  >  1:811201/1‑100 (MQ=255)
                                                                                   cAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGa                  <  1:3440947/100‑1 (MQ=255)
                                                                                   cAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGa                  >  2:2611541/1‑100 (MQ=255)
                                                                                    aaGAGAAAACCAACCCTTACGGGTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGa                  <  1:3572490/99‑1 (MQ=255)
                                                                                     taaaaaaaCCAACCGTAAGGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAAc                >  1:2643920/3‑100 (MQ=38)
                                                                                     aGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAAc                >  1:441118/1‑100 (MQ=255)
                                                                                     aGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAAc                >  1:3610828/1‑100 (MQ=255)
                                                                                      aaaaaaaCCAACCGTAAGGGTTGGTTTTTTTATATCTGCAATTACATCGATAAACAGACCGTGACACATCACAGCCAGTTTATTTCCTGATATcccc                  >  2:412397/2‑94 (MQ=25)
                                                                                      aaaaaaaCCAACCGTAAGGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACg               >  1:1568585/2‑100 (MQ=38)
                                                                                      aaaaaaaCCAACCGTAAGGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACg               >  2:569074/2‑100 (MQ=38)
                                                                                       aaaaaaCCAACCGTAAGGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGt              >  1:3190332/1‑100 (MQ=38)
                                                                                        aaaaaCCAACCGTAAGGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTc             >  2:1049900/1‑100 (MQ=39)
                                                                                              caaacCTTACGGTTGGTTTTTTTATATCTGCAAGTAGTTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAGCGTCCAg          <  2:2303206/97‑1 (MQ=255)
                                                                                               aaccCTTAAGGTTGGTTATTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTTTTATCATAACGTCCAGAcc       <  2:107625/99‑1 (MQ=255)
                                                                                                accCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACcac     <  1:2765813/100‑1 (MQ=255)
                                                                                                accCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACcac     <  1:2134187/100‑1 (MQ=255)
                                                                                                   cTTACGGTTGGTTTTTTTTTTTTTTTTATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGTCCACAcc  >  2:860838/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
GGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACC  >  NC_000913/2286212‑2286410

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: