Predicted mutation
evidence seq id position mutation annotation gene description
RA CP000730 1,440,007 Δ1 bp coding (34/201 nt) cspA ← cold shock protein CspA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,440,0070C.93.4% 326.5 / 14.7 76coding (34/201 nt)cspAcold shock protein CspA
Reads supporting (aligned to +/- strand):  ref base C (2/3);  new base . (43/28);  total (45/31)
Fisher's exact test for biased strand distribution p-value = 3.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGGTCGCCTTCAACTACTTCAAACTCAACAGCTTGACCTTCTTCTAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATA  >  CP000730/1439873‑1440141
                                                                                                                                      |                                                                                                                                      
gcgGACGCCTTCAACTACTTCAAACTCAACAGCTTGACCTTCTTCTAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTCAGCGtt                                                                                                                                  >  2:401346/1‑141 (MQ=255)
  ggTCGCCTTCAACTACTTCAAACTCAACAGCTTGACCTTCTTCTAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTaaa                                                                                                                               >  1:301135/1‑141 (MQ=255)
  ggTCGCCTTCAACTACTTCAAACTCAACAGCTTGACCTTCTTCTAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTaaa                                                                                                                               >  1:487023/1‑141 (MQ=255)
       ccTTCAACTACTTCAAACTCAACAGCTTGACCTTCTTCTAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAACCCGAATCCTTTTT‑AGCGTTAAACCAtt                                                                                                                          >  1:171801/1‑141 (MQ=255)
               tactTCAAACTCAACAGCTTGACCTTCTTCTAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTa                                                                                                                  <  2:568175/141‑1 (MQ=255)
                   tCAAACTCAACAGCTTGACCTTCTTCTAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTCAGCGTTAAACCATTTAACTGTACCt                                                                                                               <  1:698739/141‑1 (MQ=255)
                     aaaCTCAACAGCTTGACCTTCTTCTAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGt                                                                                                            >  1:26388/1‑141 (MQ=255)
                      aaCTCAACAGCTTGACCTTCTTCTAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGtt                                                                                                           <  2:244137/141‑1 (MQ=255)
                       aCTCAACAGCTTGACCTTCTTCTAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGttt                                                                                                          >  2:694227/1‑141 (MQ=255)
                         tCAACAGCTTGACCTTCTTCTAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCa                                                                                                        <  2:526055/141‑1 (MQ=255)
                         tCAACAGCTTGACCTTCTTCTAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCAGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCa                                                                                                        <  1:268492/141‑1 (MQ=255)
                           aaCAGCTTGACCTTCTTCTAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCAATAAATCCAAAACCGTTTTCAGCGTGGAACAATTTAAATGTACCTTGTTTCAt                                                                                                       >  1:106215/1‑141 (MQ=38)
                              aGCTTGACCTTCTTCTAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATc                                                                                                   <  1:278197/141‑1 (MQ=255)
                               gCTTGACCTTCTTCTAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCt                                                                                                  >  1:148696/1‑141 (MQ=255)
                               gCTTGACCTTCTTCTAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCt                                                                                                  >  2:367994/1‑141 (MQ=255)
                                cTTGACCTTCTTCTAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTg                                                                                                 <  1:383832/141‑1 (MQ=255)
                                cTTGACCTTCTTCTAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTg                                                                                                 >  2:123110/1‑141 (MQ=255)
                                 ttGACCTTCTTCTAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGa                                                                                                <  2:242066/141‑1 (MQ=255)
                                 ttGACCTTCTTCTAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGa                                                                                                <  1:401735/141‑1 (MQ=255)
                                 ttGACCTTCTTCTAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGa                                                                                                <  2:491604/141‑1 (MQ=255)
                                  tGACCTTCTTCTAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGaa                                                                                               >  1:160815/1‑141 (MQ=255)
                                         cttctAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCa                                                                                        >  1:516876/1‑141 (MQ=255)
                                               aTGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCGACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTa                                                                                  >  1:465240/1‑141 (MQ=255)
                                                tGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTaa                                                                                 >  1:461630/1‑141 (MQ=255)
                                                tGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTaa                                                                                 <  2:387852/141‑1 (MQ=255)
                                                  aTTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTaaaa                                                                               <  2:171801/141‑1 (MQ=255)
                                                     tGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGa                                                                                                >  1:608567/1‑121 (MQ=255)
                                                     tGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGa                                                                                                <  2:608567/121‑1 (MQ=255)
                                                      gTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAattca                                                                           >  2:14810/1‑141 (MQ=255)
                                                            aTCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAat                                                                     <  1:525538/141‑1 (MQ=255)
                                                            aTCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAat                                                                     >  2:640590/1‑141 (MQ=255)
                                                             tCTTGGTTAATTGCTGAAAAATGTACGAATTCGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAata                                                                    <  1:156677/141‑1 (MQ=255)
                                                                tGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGc                                                                 >  1:590933/1‑141 (MQ=255)
                                                                  gTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCtt                                                               <  1:298460/141‑1 (MQ=255)
                                                                   ttAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTa                                                              >  1:571338/1‑141 (MQ=255)
                                                                         gCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCa                                                        >  2:226512/1‑141 (MQ=255)
                                                                         gCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATACTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTACAAGACTAAAATTCATTCAATATGCTTATTCGCa                                                        >  2:510581/1‑141 (MQ=255)
                                                                             aaaaaTGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATAtt                                                    <  1:343077/141‑1 (MQ=255)
                                                                                aaTGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACa                                                 >  2:518129/1‑141 (MQ=255)
                                                                                     aCGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCagatta                                                                                                   >  2:577988/1‑81 (MQ=255)
                                                                                     aCGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCagatta                                                                                                   <  1:577988/86‑6 (MQ=255)
                                                                                     aCGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACaaaaaa                                            >  2:231336/1‑141 (MQ=255)
                                                                                      cGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTATAAAAAAt                                           >  2:623522/1‑141 (MQ=255)
                                                                                           aCGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTC‑GCTTTAAACCATTTAGCTGTACCTTGTTGCATAATCTAAACCCTCCAAGACTAAAATTCATTCACTATTCTTATTCGTATATTACAAAAAATACATg                                      <  2:261127/141‑1 (MQ=255)
                                                                                           aCGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATg                                      >  1:177318/1‑141 (MQ=255)
                                                                                            cGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTaaaa                                                                               <  2:253288/99‑1 (MQ=255)
                                                                                            cGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTaaaa                                                                               >  1:253288/1‑99 (MQ=255)
                                                                                            aatCAATTTCTCCTTCAACTTCGATAAAGCCGAATGCTTTTTCAGCTTTGAACCATTTAACTGTTCCTTGTTTCATAATCAGAAAACTCAAAAACAAAAATTCATTCAATATGCTTATTCACATATTACAAAAAATACATg                                      <  2:106215/139‑1 (MQ=255)
                                                                                             gTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTAGAAAAATACAAAAAATACATGCa                                    <  2:306022/141‑1 (MQ=255)
                                                                                              tCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCaa                                   <  1:195096/141‑1 (MQ=255)
                                                                                                aTTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAAt                                 >  2:256311/1‑141 (MQ=255)
                                                                                                aTTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAAt                                 >  1:222074/1‑141 (MQ=255)
                                                                                                aTTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAAt                                 >  2:256339/1‑141 (MQ=255)
                                                                                                      tcCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACt                           >  1:348589/1‑141 (MQ=255)
                                                                                                        cTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCGAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTCa                         >  2:660678/1‑141 (MQ=255)
                                                                                                          tCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCAtt                                                                                                                          >  1:79393/1‑42 (MQ=255)
                                                                                                          tCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCAtt                                                                                                                          <  2:79393/42‑1 (MQ=255)
                                                                                                               ttCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTGTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTCAAATATTc                  >  2:418632/1‑141 (MQ=255)
                                                                                                                tCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCt                 >  1:572919/1‑141 (MQ=255)
                                                                                                                 cGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAata                                                                    >  1:219093/1‑89 (MQ=255)
                                                                                                                 cGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAata                                                                    <  2:219093/89‑1 (MQ=255)
                                                                                                                 cGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCtt                <  2:345634/141‑1 (MQ=255)
                                                                                                                  gATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTa                                                                                                                  <  1:237196/42‑1 (MQ=37)
                                                                                                                  gATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTa                                                                                                                  >  2:237196/1‑42 (MQ=37)
                                                                                                                  gATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCttt               >  1:326190/1‑141 (MQ=255)
                                                                                                                  gATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCttt               >  2:287328/1‑141 (MQ=255)
                                                                                                                        gCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAata         >  1:557989/1‑141 (MQ=255)
                                                                                                                        gCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAata         >  2:636365/1‑141 (MQ=255)
                                                                                                                        gCCGAATCCTTTTT‑AGCATTAAACCATCTAACTGTACCTTGATTCATAATCTTTAACCTCCAAGACAAAAATTCATTCAATATGCTTATTCGCCTATTCCCAAAATTACATCCAAATTTCACACAAatcctgtttacataa         >  2:213619/1‑128 (MQ=255)
                                                                                                                         ccGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAatat        <  2:368272/141‑1 (MQ=255)
                                                                                                                           gAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTa                                                   >  1:195882/1‑96 (MQ=255)
                                                                                                                           gAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTa                                                   <  2:195882/96‑1 (MQ=255)
                                                                                                                           gAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATgg      <  1:138674/141‑1 (MQ=255)
                                                                                                                             aTCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGaa    >  2:240698/1‑141 (MQ=255)
                                                                                                                              tCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGaat   <  1:287328/141‑1 (MQ=255)
                                                                                                                               ccTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATGCATTCAATAGGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGaata  >  1:140354/1‑141 (MQ=255)
                                                                                                                                      |                                                                                                                                      
GCGGTCGCCTTCAACTACTTCAAACTCAACAGCTTGACCTTCTTCTAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATA  >  CP000730/1439873‑1440141

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