Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,638,733 A→C intergenic (‑66/+297) cspI ← / ← rzpQ Qin prophage; cold shock protein/Rz‑like protein, Qin prophage

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,638,7330AC98.7% 536.4 / ‑2.2 155intergenic (‑66/+297)cspI/rzpQQin prophage; cold shock protein/Rz‑like protein, Qin prophage
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base C (81/72);  total (82/73)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 159 total aligned reads displayed.

ATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACCAGAACACCATTAACGCATG  >  NC_000913/1638637‑1638829
                                                                                                |                                                                                                
aTTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGaacaac                                                                                               >  1:481907/1‑100 (MQ=255)
 tttCACTAAACCAGTTATTTTGTTAGATATAATTAATACCTGTTGAATAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAACAAcc                                                                                              <  2:3377980/100‑1 (MQ=255)
   tCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGACAGAACAACcaa                                                                                            >  2:3193427/1‑100 (MQ=255)
    cACTAAACCAGTCATTTTTTTTGACATAATTATTACCATATGAAGAAATTAGCCCTTTGGCAGAATGGTCCGAAAAAAATTATCAGAAAGAAcaaccaac                                                                                           >  2:3906165/1‑100 (MQ=255)
    cACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAcaaccaac                                                                                           >  1:3945108/1‑100 (MQ=255)
     aCTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGATTGGTCCGAAAAAAAATATCGCAGAGAACAACcaaca                                                                                          >  2:1818412/1‑100 (MQ=255)
      cTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAACAACcaacaa                                                                                         >  2:342659/1‑100 (MQ=255)
      cTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAACAACcaacaa                                                                                         >  2:1890387/1‑100 (MQ=255)
        aaaCCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAACAACCAACAAgg                                                                                       <  2:3643825/100‑1 (MQ=255)
        aaaCCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAACAACCAACAAgg                                                                                       <  1:3633712/100‑1 (MQ=255)
           ccAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAACAACCAACAAGGaaa                                                                                    <  2:2965454/100‑1 (MQ=255)
             aGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAACAACCAACAAGGAAAtc                                                                                  <  1:152382/100‑1 (MQ=255)
             aGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAACAACCAACAAGGAAAtc                                                                                  <  1:2633380/100‑1 (MQ=255)
              gTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAACAACCAACAAGGAAAtct                                                                                 >  2:2256563/1‑100 (MQ=255)
              gTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAACAACCAACAAGGAAAtct                                                                                 >  2:1500783/1‑100 (MQ=255)
               tCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAACAACCAACAAGGAAAtctc                                                                                <  1:3721303/100‑1 (MQ=255)
                 aTTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAACAACCAACAAGGAAATCTCaa                                                                              >  1:3767792/1‑100 (MQ=255)
                  ttttttTATACATAATAATTACCTTTTGAAGAAATTAGCCTTTGGGCAGAAAAAAAAAAAAAAAAATATCAGAGAGAACAACCAACAAGGAAATCTCAag                                                                             <  2:2492417/100‑1 (MQ=255)
                   tttGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAACAACCAACAAGGAAATCTCAaga                                                                            <  2:2550547/100‑1 (MQ=255)
                     tGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAACAACCAACAAGGAAATCTCAAGAgg                                                                          <  2:764155/100‑1 (MQ=255)
                      gTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAACAACCAACAAGGAAATCTCAAGAGGt                                                                         >  1:1605917/1‑100 (MQ=255)
                       ttAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTa                                                                        <  1:3594321/100‑1 (MQ=255)
                        tAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAACAACCAACAAGGAAATCTCAAGAGGTAc                                                                       >  1:691386/1‑100 (MQ=255)
                         aGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAACAACCAACAAGGAAATCTCAAGAGGTACa                                                                      >  2:2646191/1‑100 (MQ=255)
                           aCATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAACAACCAACAAGGAAATCTCAAGAGGTACaaa                                                                    <  2:679283/100‑1 (MQ=255)
                           aCATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAACAACCAACAAGGAAATCTCAAGAGGTACaaa                                                                    <  1:1975858/100‑1 (MQ=255)
                           aCATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCAAAAAAAAATATCAGAGAGAACAACCAACAAGGAAATATCAAGAGGTACaaa                                                                    >  1:1821171/1‑100 (MQ=255)
                           aCATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGACAGAAAAAAAATATCAGAGAGAACAACCAACAAGGAAATCTCAAGAGGTACaaa                                                                    <  2:679065/100‑1 (MQ=255)
                            cATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAACAACCAACAAGGAAATCTCAAGAGGTACAaat                                                                   >  1:81308/1‑100 (MQ=255)
                            cATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAACAACCAACAAGGAAATCTCAAGAGGTACAaat                                                                   <  2:2964357/100‑1 (MQ=255)
                            cATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAACAACCAACAAGGAAATCTCAAGAGGTACAaat                                                                   >  2:1955942/1‑100 (MQ=255)
                            cATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAACAACCAACAAGGAAATCTCAAGAGGTACAaat                                                                   >  2:56959/1‑100 (MQ=255)
                             aTAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAACAACCAACAAGGAAATCTCAAGAGGTACAaatn                                                                  <  1:3243748/100‑2 (MQ=255)
                             aTAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAACAACCAACAAGGAAATCTCAAGAGGTACAaata                                                                  >  1:3636119/1‑100 (MQ=255)
                              tAATTATTACCTTTTGAAGAAATTAGCCCTTTGGAAAAATGATCCGAAAAAAAATATCAGAAAGAACAACCAACAAGGAAATCTCAAGAGGTACAaataa                                                                 <  2:3876479/100‑1 (MQ=255)
                              tAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAACAACCAACAAGGAAATCTCAAGAGGTACAaataa                                                                 >  2:3442407/1‑100 (MQ=255)
                              tAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAACAACCAACAAGGAAATCTCAAGAGGTACAaataa                                                                 >  2:2742367/1‑100 (MQ=255)
                              tAATTATTACCATTTGAAGAAAGCAGCCCTTGGGAAGAAAAAAAAAAAAAAAAATATCAGAGAGAACAACCAACAAGGAAATCTCAAGAGGTACAaataa                                                                 <  2:1332119/100‑1 (MQ=255)
                               aaTTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAACAACCAACAAGGAAATCTCAAGAGGTACAaataat                                                                <  1:138708/100‑1 (MQ=255)
                                attattACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAACAACCAACAAGGAAATCTCAAGAGGTACAAataata                                                               <  1:4033869/100‑1 (MQ=255)
                                 ttattaCCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAACAACCAACAAGGAAATCTCAAGAGGTACAAAtaataa                                                              <  2:2393099/100‑1 (MQ=255)
                                  tattaCCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAACAACCAACAAGGAAATCTCAAGAGGTACAAATAATaaa                                                             <  1:1127355/100‑1 (MQ=255)
                                    ttaCCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAACAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAAt                                                           >  1:2983436/1‑100 (MQ=255)
                                    ttaCCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAACAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAAt                                                           <  1:1259603/100‑1 (MQ=255)
                                      aCCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAACAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATta                                                         >  1:3104277/1‑100 (MQ=255)
                                       ccTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAACAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATtat                                                        <  1:1261005/100‑1 (MQ=255)
                                       ccTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAACAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATtat                                                        <  1:3307686/100‑1 (MQ=255)
                                       ccTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATtat                                                        >  2:1817410/1‑100 (MQ=255)
                                        cTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAACAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATtata                                                       >  1:1898234/1‑100 (MQ=255)
                                        cTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAACAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATtata                                                       <  1:194459/100‑1 (MQ=255)
                                           ttGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAACAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACa                                                    <  1:4054318/100‑1 (MQ=255)
                                            tGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAACAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACaa                                                   >  1:3780864/1‑100 (MQ=255)
                                            tGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAACAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACaa                                                   >  2:1435823/1‑100 (MQ=255)
                                             gcagaaACTAGCCCTTGCGGAGAATGGTCCGAAAAAAAATATCAGAGCGAACAACCAACAACGAAATCTCCAGAGGTACAACTAATAAAATTATAACAAt                                                  >  2:323305/3‑100 (MQ=255)
                                             gaagaaATTAGCCCTTGGGCGGAATGGTCCGAAAAAAAATATCAGAGAGAACAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAAt                                                  >  1:1618271/1‑100 (MQ=255)
                                             gaagaaATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAACAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAAt                                                  >  1:1546388/1‑100 (MQ=255)
                                             gaagaaATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAACAACCAACAAGGAAATATCAAGAGGTACAAATAATAAAATTATAACAAt                                                  >  1:3742856/1‑100 (MQ=255)
                                              aagaaATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAACAACCAACAAGGAAATCTCAAGAGGGACAAATAATAAAATTATAACAATg                                                 >  2:856195/1‑100 (MQ=255)
                                              aagaaATGAGCCCTAAAGCAGAAAGGGCAGAAAAAAAATATCAGAGAGAACAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAAtn                                                 <  1:3243154/100‑2 (MQ=255)
                                                gaaATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAACAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGAc                                               <  1:1443646/100‑1 (MQ=255)
                                                 aaaTTAGCCCTTGGGAAGAATGGTACGAAAAAAAATATCAGAGAGAACAACCAACAAGGAAAGCTCAAGTGGTACAAATAATAAAATTATAACAATGACt                                              <  2:2824372/100‑1 (MQ=255)
                                                  aaTTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAACAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTg                                             >  2:19721/1‑100 (MQ=255)
                                                    ttAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAACAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCt                                           <  1:2754713/100‑1 (MQ=255)
                                                     tAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAACAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCtt                                          <  2:1018578/100‑1 (MQ=255)
                                                      aGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAACAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTc                                         <  1:147986/100‑1 (MQ=255)
                                                      aGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAACAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTc                                         <  1:3481116/100‑1 (MQ=255)
                                                       gCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAACAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCa                                        <  1:3042588/100‑1 (MQ=255)
                                                              ggCAGAATGGTCCGAAAAAAAATATCAGAGAGAACAACCAACAAGGAAATCTCAAGAGGTACAAATAATGAAATTATAACAATGACTGCTTCAGATAAAt                                 >  2:739077/1‑100 (MQ=255)
                                                              ggCAGAATGGTCCGAAAAAAAATATCAGAGAGAACAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAAt                                 >  1:1751019/1‑100 (MQ=255)
                                                                cAGAATGGTCCGAAAAAAAATATCAGAGAGAACAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAAttt                               <  2:2886264/100‑1 (MQ=255)
                                                                 aGAATGGTCCGAAAAAAAATATCAGAGAGAACAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTg                              <  1:2843632/100‑1 (MQ=255)
                                                                 aGAATGGTCCGAAAAAAAATATCAGAGAGAACAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTg                              >  2:3702954/1‑100 (MQ=255)
                                                                    aTGGTCCGAAAAAAAATATCAGAGAGAACAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTaa                           >  1:1717327/1‑100 (MQ=255)
                                                                     tGGTCCGAAAAAAAATATCAGAGAGAACAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAAc                          <  1:821687/100‑1 (MQ=255)
                                                                     tGGTCCGAAAAAAAATATCAGAGAGAACAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAAc                          <  1:2291603/100‑1 (MQ=255)
                                                                       gTCCGAAAAAAAATATCAGAGAGAACAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACaa                        <  1:3345750/100‑1 (MQ=255)
                                                                         ccGAAAAAAAATATCAGAGAGAACAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAAc                      >  1:1173527/1‑100 (MQ=255)
                                                                         ccGAAAAAAAATATCAGAGAGAACAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAAc                      >  2:1957512/1‑100 (MQ=255)
                                                                         ccGAAAAAAAATATCAGAGAGAACAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAAc                      >  1:1623857/1‑100 (MQ=255)
                                                                         ccGAAAAAAAATATCAGAGAGAACAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAAc                      >  1:3367892/1‑100 (MQ=255)
                                                                         aaaaaaaaaaaTATCAGAGAGAACAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAAc                      <  2:1001219/97‑1 (MQ=255)
                                                                           gAAAAAAAATATCAGAGAGAACAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACCa                    >  2:3380189/1‑100 (MQ=255)
                                                                            aaaaaaaaTATCAGAGAGAACAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACCAg                   <  2:3479833/100‑1 (MQ=255)
                                                                             aaaaaaaTATCAGAGAGAACAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACCAGa                  <  2:47873/100‑1 (MQ=255)
                                                                              aaaaaaTATCAGAGAGAACAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACCAGaa                 <  1:2021582/100‑1 (MQ=255)
                                                                                aaaaTATCAGAGAGAACAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACCAGAaca               <  2:864821/100‑1 (MQ=255)
                                                                                aaaaTATCAGAGAGAACAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACCAGAaca               >  2:2764260/1‑100 (MQ=255)
                                                                                aaaaTATCAGAGAGAACAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACCAGAaca               <  2:702885/100‑1 (MQ=255)
                                                                                  aaTATCAGAGAGAACAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACCAGAACAcc             <  1:3463626/100‑1 (MQ=255)
                                                                                  aaTATCAGAGAGAACAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACCAGAACAcc             <  1:2891603/100‑1 (MQ=255)
                                                                                  aaTATCAGAGAGAACAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACCAGAACAcc             <  1:2018350/100‑1 (MQ=255)
                                                                                     atCAGAGAGAACAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACCAGAACACCAtt          >  2:1144588/1‑100 (MQ=255)
                                                                                       cAGAGAGAACAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACCAGAACACCATTaa        >  1:572086/1‑100 (MQ=255)
                                                                                        agagagAACAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACCAGAACACCATTAAc       >  2:3862022/1‑100 (MQ=255)
                                                                                        agagagAACAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACCAGAACACCATTAAc       >  1:2906921/1‑100 (MQ=255)
                                                                                         gagagaACAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACCAGAACACCATTAACg      >  2:727863/1‑100 (MQ=255)
                                                                                         gagagaACAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACCAGAACACCATTAACg      <  2:1309177/100‑1 (MQ=255)
                                                                                           gagaACAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACCAGAACACCATTAACGCa    <  1:1842538/100‑1 (MQ=255)
                                                                                            agaACAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACCAGAACACCATTAACGCAt   >  2:2127220/1‑100 (MQ=255)
                                                                                             gaACAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACCAGAACACCATTAACGCATg  >  1:2287783/1‑100 (MQ=255)
                                                                                                |                                                                                                
ATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACCAGAACACCATTAACGCATG  >  NC_000913/1638637‑1638829

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