Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,883,098 (C)6→7 coding (1041/1509 nt) casA ← type I‑E CRISPR system Cascade subunit CasA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,883,0921.C97.6% 306.7 / ‑0.3 84coding (1047/1509 nt)casAtype I‑E CRISPR system Cascade subunit CasA
Reads supporting (aligned to +/- strand):  ref base . (2/0);  new base C (43/39);  total (45/39)
Fisher's exact test for biased strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.07e-01

GCCTTGCGTAAGGCTGTTTTATATCCCAAACCAACAGTCACTATTTCGTTTATCACATTGCCGTATTGTTGCCACCCCTGATTAAACATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCATAATCAATTCAAGAGGACTTTGCGGCGCAATATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGAATAATCTTATCTACCACAACTCGGCTGATTTGTGTCCATGATGGTGCG  >  NC_000913/2882955‑2883230
                                                                                                                                          |                                                                                                                                           
gCCTTGCGTAAGGCTGTTTTATATCCCAAACCAACAGTCACTATTTCGTTTATCACATTGCCGTATTGTTGCCACCCCTGATTAAACATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATcc                                                                                                                                           <  1:209564/140‑1 (MQ=255)
gCCTTGCGTAAGGCTGTTTTATATCCCAAACCAACAGTCACTATTTCGTTTATCACATTGCCGTATTGTTGCCACCCCTGATTAAACATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATcc                                                                                                                                           <  1:892974/140‑1 (MQ=255)
gCCTTGCGTAAGGCTGTTTTATATCCCAAACCAACAGTCACTATTTCGTTTATCACATTGCCGTATTGTTGCCACCCCTGATTAAACATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATcc                                                                                                                                           <  2:163164/140‑1 (MQ=255)
gCCTTGCGTAAGGCTGTTTTATATCCCAAACCAACAGTCACTATTTCGTTTATCACATTGCCGTATTGTTGCCACCCCTGATTAAACATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATcc                                                                                                                                           <  2:1077012/140‑1 (MQ=255)
 ccTTGCGTAAGGCTGTTTTATATCCCAAACCAACAGTCACTATTTCGTTTATCACATTGCCGTATTGTTGCCACCCCTGATTAAACATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATccc                                                                                                                                          >  1:727732/1‑140 (MQ=255)
 ccTTGCGTAAGGCTGTTTTATATCCCAAACCAACAGTCACTATTTCGTTTATCACATTGCCGTATTGTTGCCACCCCTGATTAAACATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATccc                                                                                                                                          >  1:637598/1‑140 (MQ=255)
  cTTGCGTAAGGCTGTTTTATATCCCAAACCAACAGTCACTATTTCGTTTATCACATTGCCGTATTGTTGCCACCCCTGATTAAACATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATcccc                                                                                                                                         >  2:955992/1‑140 (MQ=255)
  cTTGCGTAAGGCTGTTTTATATCCCAAACCAACAGTCACTATTTCGTTTATCACATTGCCGTATTGTTGCCACCCCTGATTAAACATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATcccc                                                                                                                                         >  1:1187685/1‑140 (MQ=255)
  cTTGCGTAAGGCTGTTTTATATCCCAAACCAACAGTCACTATTTCGTTTATCACATTGCCGTATTGTTGCCACCCCTGATTAAACATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATcccc                                                                                                                                         >  2:326072/1‑140 (MQ=255)
    tGCGTAAGGCTGTTTTATATCCCAAACCAACAGTCACTATTTCGTTTATCACATTGCCGTATTGTTGCCACCCCTGATTAAACATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATcccccc                                                                                                                                       >  2:903006/1‑140 (MQ=255)
    tGCGTAAGGCTGTTTTATATCCCAAACCAACAGTCACTATTTCGTTTATCACATTGCCGTATTGTTGCCACCCCTGATTAAACATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATcccccc                                                                                                                                       >  2:454662/1‑140 (MQ=255)
      cGTAAGGCTGTTTTATATCCCAAACCAACAGTCACTATTTCGTTTATCACATTGCCGTATTGTTGCCACCCCTGATTAAACATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATccccccca                                                                                                                                     >  2:746213/1‑138 (MQ=255)
          aGGCTGTTTTATATCCCAAACCAACAGTCACTATTTCGTTTATCACATTGCCGTATTGTTGCCACCCCTGATTAAACATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCCATAAt                                                                                                                                 <  1:326072/140‑1 (MQ=255)
          aGGCTGTTTTATATCCCAAACCAACAGTCACTATTTCGTTTATCACATTGCCGTATTGTTGCCACCCCTGATTAAACATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCCATAAt                                                                                                                                 <  1:315557/140‑1 (MQ=255)
             cTGTTTTATATCCCAAACCAACAGTCACTATTTCGTTTATCACATTGCCGTATTGTTGCCACCCCTGATTAAACATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCCATAATCaa                                                                                                                              <  1:879828/140‑1 (MQ=255)
                ttttATATCCCAAACCAACAGTCACTATTTCGTTTATCACATTGCCGTATTGTTGCCACCCCTGATTAAACATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCCATAATCAATTc                                                                                                                           >  1:1242388/1‑140 (MQ=255)
                ttttATATCCCAAACCAACAGTCACTATTTCGTTTATCACATTGCCGTATTGTTGCCACCCCTGATTAAACATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATcc                                                                                                                                           >  1:226663/1‑124 (MQ=255)
                ttttATATCCCAAACCAACAGTCACTATTTCGTTTATCACATTGCCGTATTGTTGCCACCCCTGATTAAACATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATcc                                                                                                                                           <  2:226663/124‑1 (MQ=255)
                  ttATATCCCAAACCAACAGTCACTATTTCGTTTATCACATTGCCGTATTGTTGCCACCCCTGATTAAACATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCCATAATCAATTCaa                                                                                                                         <  1:746213/140‑1 (MQ=255)
                            aaCCAACAGTCACTATTTCGTTTATCACATTGCCGTATTGTTGCCACCCCTGATTAAACATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCCATAATCAATTCAAGAGGACTTTg                                                                                                               >  2:490922/1‑140 (MQ=255)
                            aaCCAACAGTCACTATTTCGTTTATCACATTGCCGTATTGTTGCCACCCCTGATTAAACATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCCATAATCAATTCAAGAGGACTTTg                                                                                                               >  2:624046/1‑140 (MQ=255)
                               cAACAGTCACTATTTCGTTTATCACATTGCCGTATTGTTGCCACCCCTGATTAAACATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCCATAATCAATTCAAGAGGACTTTgcgg                                                                                                            <  1:624046/140‑1 (MQ=255)
                               cAACAGTCACTATTTCGTTTATCACATTGCCGTATTGTTGCCACCCCTGATTAAACATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCCATAATCAATTCAAGAGGACTTTgcgg                                                                                                            <  1:903006/140‑1 (MQ=255)
                                aaCAGTCACTATTTCGTTTATCACATTGCCGTATTGTTGCCACCCCTGATTAAACATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCCATAATCAATTCAAGAGGACTTTgcggc                                                                                                           <  2:1036895/140‑1 (MQ=255)
                                 aCAGTCACTATTTCGTTTATCACATTGCCGTATTGTTGCCACCCCTGATTAAACATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCCATAATCAATTCAAGAGGACTTTgcggcg                                                                                                          <  2:45992/140‑1 (MQ=255)
                                  cAGTCACTATTTCGTTTATCACATTGCCGTATTGTTGCCACCCCTGATTAAACATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCCATAATCAATTCAAGAGGACTTTGCGgcgc                                                                                                         >  1:793991/1‑140 (MQ=255)
                                    gTCACTATTTCGTTTATCACATTGCCGTATTGTTGCCACCCCTGATTAAACATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCCATAATCAATTCAAGAGGACTTTGCGGCGCaa                                                                                                       <  1:812286/140‑1 (MQ=255)
                                       aCTATTTCGTTTATCACATTGCCGTATTGTTGCCACCCCTGATTAAACATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCCATAATCAATTCAAGAGGACTTTGCGGCGCaa                                                                                                       <  2:394554/137‑1 (MQ=255)
                                       aCTATTTCGTTTATCACATTGCCGTATTGTTGCCACCCCTGATTAAACATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCCATAATCAATTCAAGAGGACTTTGCGGCGCaa                                                                                                       >  1:394554/1‑137 (MQ=255)
                                       aCTATTTCGTTTATCACATTGCCGTATTGTTGCCACCCCTGATTAAACATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCCATAATCAATTCAAGAGGACTTTGCGGCGCAatat                                                                                                    >  1:1022245/1‑140 (MQ=255)
                                       aCTATTTCGTTTATCACATTGCCGTATTGTTGCCACCCCTGATTAAACATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCCATAATCAATTCAAGAGGACTTTGCGGCGCAatat                                                                                                    >  1:1223678/1‑140 (MQ=255)
                                         tATTTCGTTTATCACATTGCCGTATTGTTGCCACCCCTGATTAAACATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCCATAATCAATTCAAGAGGACTTTGCGGCGCAATAttt                                                                                                  <  2:1022245/140‑1 (MQ=255)
                                               gTTTATCACATTGCCGTATTGTTGCCACCCCTGATTAAACATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCCATAATCAATTCAAGAGGACTTTGCGGCGCAATATTTCTGAAt                                                                                            <  1:454662/140‑1 (MQ=255)
                                               gTTTATCACATTGCCGTATTGTTGCCACCCCTGATTAAACATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCCATAATCAATTCAAGAGGACTTTGCGGCGCAATATTTCTGAAt                                                                                            <  2:1242388/140‑1 (MQ=255)
                                                     cacaTTGCCGTATTGTTGCCACCCCTGATTAAACATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCCATAATCAATTCAAGAGGACTTTGCGGCGCAATATTTCTGAATTGATTc                                                                                      >  1:183132/1‑140 (MQ=255)
                                                            ccGTATTGTTGCCACCCCTGATTAAACATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCCATAATCAATTCAAGAGGACTTTGCGGCGCAATATTTCTGAATTGATTCACAAccg                                                                               >  1:1279465/1‑140 (MQ=255)
                                                                      gCCACCCCTGATTAAACATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCCATAATCAATTCAAGAGGACTTTGCGGCGCAATATTTCTGAATTGATTCACAACCGCCGCCAcgcg                                                                     >  2:1267133/1‑140 (MQ=255)
                                                                      gCCACCCCTGATTAAACATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCCATAATCAATTCAAGAGGACTTTGCGGCGCAATATTTCTGAATTGATTCACAACCGCCGCCAcgcg                                                                     >  1:642717/1‑140 (MQ=255)
                                                                      gCCACCCCTGATTAAACATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCCATAATCAATTCAAGAGGACTTTGCGGCGCAATATTTCTGAATTGATTCACAACCGCCGCCAcgcg                                                                     >  2:23865/1‑140 (MQ=255)
                                                                      gCCACCCCTGATTAAACATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCCATAATCAATTCAAGAGGACTTTGCGGCGCAATATTTCTGAATTGATTCACAACCGCCGCCAcgcg                                                                     >  2:533219/1‑140 (MQ=255)
                                                                      gCCACCCCTGATTAAACATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCCATAATCAATTCAAGAGGACTTTGCGGCGCAATATTTCTGAATTGATTCACAACCGCCGCCAcgcg                                                                     >  1:413897/1‑140 (MQ=255)
                                                                      gCCACCCCTGATTAAACATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCCATAATCAATTCAAGAGGACTTTGCGGCGCAATATTTCTGAATTGATTCACAACCGCCGCCAcgcg                                                                     >  1:119234/1‑140 (MQ=255)
                                                                         aCCCCTGATTAAACATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCCATAATCAATTCAAGAGGACTTTGCGGCGCAATATTTCTGAATTGATTCACAACCGCCGCCACGCGAtt                                                                  >  2:902873/1‑140 (MQ=255)
                                                                         aCCCCTGATTAAACATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCCATAATCAATTCAAGAGGACTTTGCGGCGCAATATTTCTGAATTGATTCACAACCGCCGCCACGCGAtt                                                                  <  1:490922/140‑1 (MQ=255)
                                                                           cccTGATTAAACATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCCATAATCAATTCAAGAGGACTTTGCGGCGCAATATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTc                                                                >  1:266376/1‑140 (MQ=255)
                                                                           cccTGATTAAACATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCCATAATCAATTCAAGAGGACTTTGCGGCGCAATATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTc                                                                >  2:846528/1‑140 (MQ=255)
                                                                              tGATTAAACATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCCATAATCAATTCAAGAGGACTTTGCGGCGCAATATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCcat                                                             >  1:682612/1‑140 (MQ=255)
                                                                                aTTAAACATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCCATAATCAATTCAAGAGGACTTTGCGGCGCAATATTTCTGAATTGATTCACAACCGCCGCCACGCGa                                                                    >  2:986231/1‑131 (MQ=255)
                                                                                aTTAAACATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCCATAATCAATTCAAGAGGACTTTGCGGCGCAATATTTCTGAATTGATTCACAACCGCCGCCACGCGa                                                                    <  1:986231/131‑1 (MQ=255)
                                                                                    aaCATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCCATAATCAATTCAAGAGGACTTTGCGGCGCAATATTTCTGAATTGATTCACAATCGCCGCCACGCGATTTCcattttcat                                                       <  2:637598/140‑1 (MQ=255)
                                                                                     aCATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCCATAATCAATTCAAGAGGACTTTGCGGCGCAATATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCcattttcatt                                                      <  2:665094/140‑1 (MQ=255)
                                                                                      cATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCCATAATCAATTCAAGAGGACTTTGCGGCGCAATATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCcattttcattt                                                     >  1:1016314/1‑140 (MQ=255)
                                                                                         cAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCCATAATCAATTCAAGAGGACTTTGCGGCGCAATATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGa                                                  >  2:874073/1‑140 (MQ=255)
                                                                                         cAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCCATAATCAATTCAAGAGGACTTTGCGGCGCAATATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGa                                                  >  2:724238/1‑140 (MQ=255)
                                                                                         cAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCCATAATCAATTCAAGAGGACTTTGCGGCGCAATATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGa                                                  >  2:769972/1‑140 (MQ=255)
                                                                                          aaCACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCCATAATCAATTCAAGAGGACTTTGCGGCGCAATATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCcattttcatttttaa                                                 <  2:266376/140‑4 (MQ=255)
                                                                                          aaCACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCCATAATCAATTCAAGAGGACTTTGCGGCGCAATATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGaa                                                 <  2:793991/140‑1 (MQ=255)
                                                                                          aaCACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCCATAATCAATTCAAGAGGACTTTGCGGCGCAATATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGaa                                                 <  1:1197278/140‑1 (MQ=255)
                                                                                          aaCACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCCATAATCAATTCAAGAGGACTTTGCGGCGCAATATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGaa                                                 <  2:119234/140‑1 (MQ=255)
                                                                                            cacaTCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCCATAATCAATTCAAGAGGACTTTGCGGCGCAATATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGaata                                               <  2:682612/140‑1 (MQ=255)
                                                                                               atcatGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCCATAATCAATTCAAGAGGACTTTGCGGCGCAATATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGAATAATc                                            >  2:585515/1‑140 (MQ=255)
                                                                                                 catGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCCATAATCAATTCAAGAGGACTTTGCGGCGCAATATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGAATAATCtt                                          >  1:1171378/1‑140 (MQ=255)
                                                                                                  atGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCCATAATCAATTCAAGAGGACTTTGCGGCGCAATATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGAATAATCTTa                                         <  1:1146216/140‑1 (MQ=255)
                                                                                                       gCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCCATAATCAATTCAAGAGGACTTTGCGGCGCAATATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGAATAATCTTATCTAc                                    <  2:1187685/140‑1 (MQ=255)
                                                                                                         cGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCCATAATCAATTCAAGAGGACTTTGCGGCGCAATATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGAATAATc                                            >  2:1268411/1‑130 (MQ=255)
                                                                                                         cGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCCATAATCAATTCAAGAGGACTTTGCGGCGCAATATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGAATAATc                                            <  1:1268411/130‑1 (MQ=255)
                                                                                                         cGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCCATAATCAATTCAAGAGGACTTTGCGGCGCAATATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGAATAATCTTATCTACca                                  >  2:995028/1‑140 (MQ=255)
                                                                                                         cGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCCATAATCAATTCAAGAGGACTTTGCGGCGCAATATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGAATAATCTTATCTACca                                  >  2:897052/1‑140 (MQ=255)
                                                                                                           ttCAAGAATAGATGCTTGATTATTACGATATCCCCCCCATAATCAATTCAAGAGGACTTTGCGGCGCAATATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGaa                                                 <  2:1262694/123‑1 (MQ=255)
                                                                                                           ttCAAGAATAGATGCTTGATTATTACGATATCCCCCCCATAATCAATTCAAGAGGACTTTGCGGCGCAATATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGaa                                                 >  1:1262694/1‑123 (MQ=255)
                                                                                                           ttCAAGAATAGATGCTTGATTATTACGATATCCCCCCCATAATCAATTCAAGAGGACTTTGCGGCGCAATATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGAATAATCTTATCTACcaca                                >  1:979431/1‑140 (MQ=255)
                                                                                                           ttCAAGAATAGATGCTTGATTATTACGATATCCCCCCCATAATCAATTCAAGAGGACTTTGCGGCGCAATATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGAATAATCTTATCTACcaca                                >  1:917800/1‑140 (MQ=255)
                                                                                                           ttCAAGAATAGATGCTTGATTATTACGATATCCCCCCCATAATCAATTCAAGAGGACTTTGCGGCGCAATATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGAATAATCTTATCTACcaca                                >  2:716717/1‑140 (MQ=255)
                                                                                                           ttCAAGAATAGATGCTTGATTATTACGATATCCCCCCCATAATCAATTCAAGAGGACTTTGCGGCGCAATATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGAATAATCTTATCTACcaca                                >  2:691996/1‑140 (MQ=255)
                                                                                                           ttCAAGAATAGATGCTTGATTATTACGATATCCCCCCATAATCAATTCAAGAGGACTTTGCGGCGCAATATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGAATAATCTTATCTACCACaa                               >  2:256824/1‑140 (MQ=255)
                                                                                                              aaGAATAGATGCTTGATTATTACGATATCCCCCCCATAATCAATTCAAGAGGACTTTGCGGCGCAATATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGAATAATCTTATCTACCACAACt                             <  2:642717/140‑1 (MQ=255)
                                                                                                                gAATAGATGCTTGATTATTACGATATCCCCCCCATAATCAATTCAAGAGGACTTTGCGGCGCAATATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGAATAATCTTATCTACCACAACTCg                           <  2:979431/140‑1 (MQ=255)
                                                                                                                   tAGATGCTTGATTATTACGATATCCCCCCCATAATCAATTCAAGAGGACTTTGCGGCGCAATATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGAATAATCTTATCTACCACAACTCGGCt                        <  2:727732/140‑1 (MQ=255)
                                                                                                                    aGATGCTTGATTATTACGATATCCCCCCCATAATCAATTCAAGAGGACTTTGCGGCGCAATATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGAATAATCTTATCTACCACAACTCGGCTg                       <  2:917800/140‑1 (MQ=255)
                                                                                                                     gATGCTTGATTATTACGATATCCCCCCCATAATCAATTCAAGAGGACTTTGCGGCGCAATATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGAATAATCTTATCTACCACAACTCGGCTGa                      >  2:112886/1‑140 (MQ=255)
                                                                                                                     gATGCTTGATTATTACGATATCCCCCCCATAATCAATTCAAGAGGACTTTGCGGCGCAATATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGAATAATCTTATCTACCACAACTCGGCTGa                      >  1:760893/1‑140 (MQ=255)
                                                                                                                      aTGCTTGATTATTACGATATCCCCCCCATAATCAATTCAAGAGGACTTTGCGGCGCAATATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGAATAATCTTATCTACCACAACTCGGCTGAt                     >  2:516743/1‑140 (MQ=255)
                                                                                                                      aTGCTTGATTATTACGATATCCCCCCCATAATCAATTCAAGAGGACTTTGCGGCGCAATATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGAATAATCTTATCTACCACAACTCGGCTGAt                     >  1:1252096/1‑140 (MQ=255)
                                                                                                                       tGCTTGATTATTACGATATCCCCCCCATAATCAATTCAAGAGGACTTTGCGGCGCAATATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGAATAATCTTATCTACCACAACTCGGCTGAtt                    <  1:1267133/140‑1 (MQ=255)
                                                                                                                        gCTTGATTATTACGATATCCCCCCCATAATCAATTCAAGAGGACTTTGCGGCGCAATATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGAATAATCTTATCTACCACAACTCGGCTGAttt                   <  1:533219/140‑1 (MQ=255)
                                                                                                                        gCTTGATTATTACGATATCCCCCCCATAATCAATTCAAGAGGACTTTGCGGCGCAATATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGAATAATCTTATCTACCACAACTCGGCTGAttt                   <  1:112886/140‑1 (MQ=255)
                                                                                                                        gCTTGATTATTACGATATCCCCCCCATAATCAATTCAAGAGGACTTTGCGGCGCAATATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGAATAATCTTATCTACCACAACTCGGCTGAttt                   <  1:691996/140‑1 (MQ=255)
                                                                                                                        gCTTGATTATTACGATATCCCCCCCATAATCAATTCAAGAGGACTTTGCGGCGCAATATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGAATAATCTTATCTACCACAACTCGGCTGAttt                   >  2:292368/1‑140 (MQ=255)
                                                                                                                          ttGATTATTACGATATCCCCCCCATAATCAATTCAAGAGGACTTTGCGGCGCAATATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGAATAATCTTATCTACCACAACTCGGCTGATTtgt                 <  1:615927/140‑1 (MQ=255)
                                                                                                                           tGATTATTACGATATCCCCCCCATAATCAATTCAAGAGGACTTTGCGGCGCAATATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGAATAATCTTATCTACCACAACTCGGCTGATTtgtg                >  1:764644/1‑140 (MQ=255)
                                                                                                                           tGATTATTACGATATCCCCCCCATAATCAATTCAAGAGGACTTTGCGGCGCAATATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGAATAATCTTATCTACCACAACTCGGCTGATTtgtg                >  1:195283/1‑140 (MQ=255)
                                                                                                                           tGATTATTACGATATCCCCCCCATAATCAATTCAAGAGGACTTTGCGGCGCAATATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGAATAATCTTATCTACCACAACTCGGCTGATTtgtg                >  1:1070088/1‑140 (MQ=255)
                                                                                                                           tGATTATTACGATATCCCCCCAATACTCAATTCAGGAGGCCTTTGCGGCGCAATTTTTCTGAATTGATTCACAACCGCCGCCACGCGTTTTCAATTTTCTTTTGGAATAATCTTATCTACAACAACTTGGCGGATTtgtt                >  2:850255/1‑139 (MQ=255)
                                                                                                                                attaCGATATCCCCCCCATAATCAATTCAAGAGGACTTTGCGGCGCAATATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGAATAATCTTATCTACCACAACTCGGCTGAt                     >  2:309033/1‑130 (MQ=255)
                                                                                                                                attaCGATATCCCCCCCATAATCAATTCAAGAGGACTTTGCGGCGCAATATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGAATAATCTTATCTACCACAACTCGGCTGAt                     <  1:309033/130‑1 (MQ=255)
                                                                                                                                 ttaCGATATCCCCCCCATAATCAATTCAAGAGGACTTTGCGGCGCAATATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGAATAATCTTATCTACCACAACTCGGCTGATTTGTGTCCatg          <  2:760893/140‑1 (MQ=255)
                                                                                                                                  tccGATATCCCCCCCATAATCAATTCAAGAGGACTTTGCGGCGCAATATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGAATAATCTTATCTACCACAACTCGGCTGATTTGTGTCCatga         <  1:48960/138‑1 (MQ=255)
                                                                                                                                  taCGATATCCCCCCCATAATCAATTCAAGAGGACTTTGCGGCGCAATATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGAATAATCTTATCTACCACAACTCGGCTGATTTGTGTCCatga         <  1:874073/140‑1 (MQ=255)
                                                                                                                                  taCGATATCCCCCCCATAATCAATTCAAGAGGACTTTGCGGCGCAATATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGAATAATCTTATCTACCACAACTCGGCTGATTTGTGTCCatga         <  1:585515/140‑1 (MQ=255)
                                                                                                                                        atCCCCCCATAATCAATTCAAGAGGACTTTGCGGCGCAATATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGAATAATCTTATCTACCACAACTCGGCTGATTTGTGTCCATGATGGTGCg  <  1:256824/140‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                           
GCCTTGCGTAAGGCTGTTTTATATCCCAAACCAACAGTCACTATTTCGTTTATCACATTGCCGTATTGTTGCCACCCCTGATTAAACATCAACACATCATGACGCCGTTCAAGAATAGATGCTTGATTATTACGATATCCCCCCATAATCAATTCAAGAGGACTTTGCGGCGCAATATTTCTGAATTGATTCACAACCGCCGCCACGCGATTTCCATTTTCATTTTGAATAATCTTATCTACCACAACTCGGCTGATTTGTGTCCATGATGGTGCG  >  NC_000913/2882955‑2883230

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