Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,852,981 A→G 5.5% L9P (CTG→CCG)  ybaB ← nucleoid‑associated broad specificity regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,852,9810AG 5.5% 219.9 / 2.2 73L9P (CTG→CCG) ybaBnucleoid‑associated broad specificity regulator
Reads supporting (aligned to +/- strand):  ref base A (31/38);  new base G (2/2);  total (33/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.30e-01

TCTGGTCGATGCTGACGCGCTTGACGTCGTGACGACCGGTCATCACCACGCTCACCAGACCGCCACCCGACTGACCGGTGACTTCGGCGTTGGCCAGTTCTTCCTGCATCTTCTGCATTTTTTCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGC  >  NC_002947/4852832‑4853129
                                                                                                                                                     |                                                                                                                                                    
tctGGTCGATGCTGACGCGCTTGACGTCGTGACGACCGGTCATCACCACGCTCACCAGACCGCCACCCGACTGACCGGTGACTTCGGCGTTGGCCAGTTCTTCCTGCATCTTCTGCATTTTTTCCTGCATCTGCTGGGCCTGCTtcatca                                                                                                                                                      >  2:1044459/1‑150 (MQ=255)
  tGGTCGATGCTGACGCGCTTGACGTCGTGACGACCGGTCATCACCACGCTCACCAGACCGCCACCCGACTGACCGGTGACTTCGGCGTTGGCCAGTTCTTCCTGCATCTTCTGCATTTTTTCCTGCATCTGCTGGGCCTGCTTCATCAgg                                                                                                                                                    >  2:318261/1‑150 (MQ=255)
     tCGATGCTGACGCGCTTGACGTCGTGACGACCGGTCATCACCACGCTCACCAGACCGCCACCCGACTGACCGGTGACTTCGGCGTTGGCCAGTTCTTCCTGCATCTTCTGCATTTTTTCCTGCATCTGCTGGGCCTGCTTCATCAggccg                                                                                                                                                 >  2:326898/1‑150 (MQ=255)
      cGATGCTGACGCGCTTGACGTCGTGACGACCGGTCATCACCACGCTCACCAGACCGCCACCCGACTGACCGGTGACTTCGGCGTTGGCCAGTTCTTCCTGCATCTTCTGCATTTTTTCCTGCATCTGCTGGGCCTGCTTCATCAggccgg                                                                                                                                                <  2:949619/150‑1 (MQ=255)
          gCTGACGCGCTTGACGTCGTGACGACCGGTCATCACCACGCTCACCAGACCGCCACCCGACTGACCGGTGACTTCGGCGTTGGCCAGTTCTTCCTGCATCTTCTGCATTTTTTCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCAt                                                                                                                                            <  1:1002062/150‑1 (MQ=255)
              aCGCGCTTGACGTCGTGACGACCGGTCATCACCACGCTCACCAGACCGCCACCCGACTGACCGGTGACTTCGGCGTTGGCCAGTTCTTCCTGCATCTTCTGCATTTTTTCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCa                                                                                                                                        <  2:1103151/150‑1 (MQ=255)
                        cgtcgtGACGACCGGTCATCACCACGCTCACCAGACCGCCACCCGACTGACCGGTGACTTCGGCGTTGGCCAGTTCTTCCTGCATCTTCTGCATTTTTTCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTtcatca                                                                                                                              <  1:201760/150‑1 (MQ=255)
                         gtcgtGACGACCGGTCATCACCACGCTCACCAGACCGCCACCCGACTGACCGGTGACTTCGGCGTTGGCCAGTTCTTCCTGCATCTTCTGCATTTTTTCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTcatcat                                                                                                                             <  1:515123/150‑1 (MQ=255)
                         gtcgtGACGACCGGTCATCACCACGCTCACCAGACCGCCACCCGACTGACCGGTGACTTCGGCGTTGGCCAGTTCTTCCTGCATCTTCTGCATTTTTTCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTcatcat                                                                                                                             <  1:191005/150‑1 (MQ=255)
                          tcgtGACGACCGGTCATCACCACGCTCACCAGACCGCCACCCGACTGACCGGTGACTTCGGCGTTGGCCAGTTCTTCCTGCATCTTCTGCATTTTTTCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATg                                                                                                                            <  1:1016096/150‑1 (MQ=255)
                          tcgtGACGACCGGTCATCACCACGCTCACCAGACCGCCACCCGACTGACCGGTGACTTCGGCGTTGGCCAGTTCTTCCTGCATCTTCTGCATTTTTTCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATg                                                                                                                            <  1:1416551/150‑1 (MQ=255)
                           cgtGACGATCTGTCATCACCCCGCTCACCAAACCGCCACCCGACTCACCAGTGACTTCGGCGTTGGCCATGTCTTCCAGCATCTTCTGCATTTTTTCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATgg                                                                                                                           <  1:99938/150‑1 (MQ=255)
                           cgtGACGACCGGTCATCACCACGCTCACCAGACCGCCACCCGACTGACCGGTGACTTCGGCGTTGGCCAGTTCTTCCTGCATCTTCTGCATTTTTTCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATgg                                                                                                                           <  2:980807/150‑1 (MQ=255)
                            gtGACGACCGGTCATCACCACGCTCACCAGACCGCCACCCGACTGACCGGTGACTTCGGCGTTGGCCAGTTCTTCCTGCATCTTCTGCATTTTTTCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATggg                                                                                                                          >  1:652992/1‑150 (MQ=255)
                                cgacCGGTCATCACCACGCTCACCAGACCGCCACCCGACTGACCGGTGACTTCGGCGTTGGCCAGTTCTTCCTGCATCTTCTGCATTTTTTCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGtat                                                                                                                      <  1:379043/150‑1 (MQ=255)
                                    cGGTCATCACCACGCTCACCAGACCGCCACCCGACTGACCGGTGACTTCGGCGTTGGCCAGTTCTTCCTGCATCTTCTGCATTTTTTCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGtatgact                                                                                                                  >  1:499286/1‑146 (MQ=255)
                                       tcatcaCAACGCTCACCAGACCGCCACCCGACTGACCGGTGACTTCGGCGTTGGCCAGTTCTTCCTGCATCTTCTGCATTTTTTCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGa                                                                                                               >  1:837056/1‑150 (MQ=255)
                                         atcaCCACGCTCACCAGACCGCCACCCGACTGACCGGTGACTTCGGCGTTGGCCAGTTCTTCCTGCATCTTCTGCATTTTTTCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCAt                                                                                                                                            <  2:60366/119‑1 (MQ=255)
                                         atcaCCACGCTCACCAGACCGCCACCCGACTGACCGGTGACTTCGGCGTTGGCCAGTTCTTCCTGCATCTTCTGCATTTTTTCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCAt                                                                                                                                            >  1:60366/1‑119 (MQ=255)
                                         atcaCCACGCTCACCAGACCGCCACCCGACTGACCGGTGACTTCGGCGTTGGCCAGTTCTTCCTGCATCTTCTGCATTTTTTCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGta                                                                                                                       >  1:1197100/1‑140 (MQ=255)
                                         atcaCCACGCTCACCAGACCGCCACCCGACTGACCGGTGACTTCGGCGTTGGCCAGTTCTTCCTGCATCTTCTGCATTTTTTCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGta                                                                                                                       <  2:1197100/140‑1 (MQ=255)
                                         atcaCCACGCTCACCAGACCGCCACCCGACTGACCGGTGACTTCGGCGTTGGCCAGTTCTTCCTGCATCTTCTGCATTTTTTCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGAtt                                                                                                             >  2:530164/1‑150 (MQ=255)
                                           caccacGCTCACCAGACCGCCACCCGACTGACCGGTGACTTCGGCGTTGGCCAGTTCTTCCTGCATCTTCTGCATTTTTTCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTgg                                                                                                           <  1:1044459/150‑1 (MQ=255)
                                           caccacGCTCACCAGACCGCCACCCGACTGACCGGTGACTTCGGCGTTGGCCAGTTCTTCCTGCATCTTCTGCATTTTTTCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTgg                                                                                                           >  1:687432/1‑150 (MQ=255)
                                                cGCTCACCAGACCGCCACCCGACTGACCGGTGACTTCGGCGTTGGCCAGTTCTTCCTGCATCTTCTGCATTTTTTCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCAt                                                                                                      <  1:318261/150‑1 (MQ=255)
                                                cGCTCACCAGACCGCCACCCGACTGACCGGTGACTTCGGCGTTGGCCAGTTCTTCCTGCATCTTCTGCATTTTTTCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCAt                                                                                                      <  2:652992/150‑1 (MQ=255)
                                                 gCTCACCAGACCGCCACCCGACTGACCGGTGACTTCGGCGTTGGCCAGTTCTTCCTGCATCTTCTGCATTTTTTCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCAc                                                                                                                                       <  2:1384592/116‑1 (MQ=255)
                                                 gCTCACCAGACCGCCACCCGACTGACCGGTGACTTCGGCGTTGGCCAGTTCTTCCTGCATCTTCTGCATTTTTTCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCAtg                                                                                                     >  2:265281/1‑150 (MQ=255)
                                                 gCTCACCAGACCGCCACCCGACTGACCGGTGACTTCGGCGTTGGCCAGTTCTTCCTGCATCTTCTGCATTTTTTCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCAATGCCAc                                                                                                                                       >  1:1384592/1‑116 (MQ=255)
                                                 ccgcACCAGACCGCCACCCGACTGACCGGTGACTTCGGCGTTGGCCAGTTCTTCCTGCATCTTCTGCATTTTTTCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCAtg                                                                                                     >  1:1342851/4‑150 (MQ=255)
                                                      ccAGACCGCCACCCGACTGACCGGTGACTTCGGCGTTGGCCAGTTCTTCCTGCATCTTCTGCATTTTTTCCTGCATCTGCTGGGCCTGCTTCATCGGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATg                                                                                                >  2:803811/1‑150 (MQ=255)
                                                      ccAGACCGCCACCCGACTGACCGGTGACTTCGGCGTTGGCCAGTTCTTCCTGCATCTTCTGCATTTTTTCCTGCATCTGCTGGGCCTGCTTCATCGGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATg                                                                                                >  2:803793/1‑150 (MQ=255)
                                                         gACCGCCACCCGACTGACCGGTGACTTCGGCGTTGGCCAGTTCTTCCTGCATCTTCTGCATTTTTTCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCgg                                                                                             >  1:579358/1‑150 (MQ=255)
                                                          aCCGCCACCCGACTGACCGGTGACTTCGGCGTTGGCCAGTTCTTCCTGCATCTTCTGCATTTTTTCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGc                                                                                            <  1:326898/150‑1 (MQ=255)
                                                            cGCCACCCGACTGACCGGTGACTTCGGCGTTGGCCAGTTCTTCCTGCATCTTCTGCATTTTTTCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGAt                                                                                                 <  2:444250/143‑1 (MQ=255)
                                                            cGCCACCCGACTGACCGGTGACTTCGGCGTTGGCCAGTTCTTCCTGCATCTTCTGCATTTTTTCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGAt                                                                                                 >  1:444250/1‑143 (MQ=255)
                                                              ccACCCGACTGACCGGTGACTTCGGCGTTGGCCAGTTCTTCCTGCATCTTCTGCATTTTTTCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCgg                                                                                        >  2:826366/1‑150 (MQ=255)
                                                                aCCCGACTGACCGGTGACTTCGGCGTTGGCCAGTTCTTCCTGCATCTTCTGCATTTTTTCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTgg                                                                                                           >  1:709291/1‑129 (MQ=255)
                                                                aCCCGACTGACCGGTGACTTCGGCGTTGGCCAGTTCTTCCTGCATCTTCTGCATTTTTTCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTgg                                                                                                           <  2:709291/129‑1 (MQ=255)
                                                                    gACTGACCGGTGACTTCGGCGTTGGCCAGTTCTTCCTGCATCTTCTGCATTTTTTCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCgg                                                                                  >  1:593296/1‑150 (MQ=255)
                                                                      cTGACCGGTGACTTCGGCGTTGGCCAGTTCTTCCTGCATCTTCTGCATTTTTTCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGTATGCCTCGATTGGGTCATGTGATGCGGCCCGGTGCCGggc                                                                                <  2:499286/150‑1 (MQ=255)
                                                                              tGACTTCGGCGTTGGCCAGTTCTTCCTGCATCTTCTGCATTTTTTCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCa                                                                        <  2:1665/150‑1 (MQ=255)
                                                                               gACTTCGGCGTTGGCCAGTTCTTTCTGCATCTTCTGCATTTTTTCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCAt                                                                       <  2:682878/150‑1 (MQ=255)
                                                                                  ttCGGCTTTGGCCAGTTCTTCCTGCATCTTCTGCATTTTTTCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGggccggc                                                                            <  2:1370491/142‑1 (MQ=255)
                                                                                  ttCGGCGTTGGCCAGTTCTTCCTGCATCTTCTGCATTTTTTCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGggccggc                                                                            >  1:1370491/1‑142 (MQ=255)
                                                                                  ttCGGCGTTGGCCAGTTCTTCCTGCATCTTCTGCATTTTTTCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGt                                                                    >  1:710750/1‑150 (MQ=255)
                                                                                  ttCGGCGTTGGCCAGTTCTTCCTGCATCTTCTGCATTTTTTCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGt                                                                    >  2:914341/1‑150 (MQ=255)
                                                                                  ttCGGCGTTGGCCAGTTCTTCCTGCATCTTCTGAATTTTTTCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGt                                                                    >  1:1146332/1‑150 (MQ=255)
                                                                                   tCGGCGTTGGCCAGTTCTTCCTGCATCTTCTGCATTTTTTCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGtt                                                                   <  2:1083211/150‑1 (MQ=255)
                                                                                         ttGGCCAGTTCTTCCTGCATCTTCTGCATTTTTTCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCAt                                                             <  2:579358/150‑1 (MQ=255)
                                                                                         ttGGCCAGTTCTTCCTGCATCTTCTGCATTTTTTCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCAt                                                             <  1:530164/150‑1 (MQ=255)
                                                                                              cAGTTCTTCCTGCATCTTCTGCATTTTTTCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTg                                                        >  1:567854/1‑150 (MQ=255)
                                                                                              cAGTTCTTCCTGCATCTTCTGCATTTTTTCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTg                                                        <  1:557254/150‑1 (MQ=255)
                                                                                                     tcCTGCATCTTCTGCATTTTTTCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGa                                                 <  2:236676/150‑1 (MQ=255)
                                                                                                       cTGCATCTTCTGCATTTTTTCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATg                                                                                                                            <  2:248793/73‑1 (MQ=255)
                                                                                                       cTGCATCTTCTGCATTTTTTCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATg                                                                                                                            >  1:248793/1‑73 (MQ=255)
                                                                                                       cTGCATCTTCTGCATTTTTTCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGtt                                                                   <  2:478279/130‑1 (MQ=255)
                                                                                                       cTGCATCTTCTGCATTTTTTCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGtt                                                                   >  1:478279/1‑130 (MQ=255)
                                                                                                       cTGCATCTTCTGCATTTTTTCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACt                                               >  2:100298/1‑150 (MQ=255)
                                                                                                             cttctGCATTTTTTCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCATGACTAGCCAg                                         <  2:1342851/150‑1 (MQ=255)
                                                                                                                ctGCATTTTTTCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGc                                      >  2:903232/1‑150 (MQ=255)
                                                                                                                   cATTTTTTCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAg                                         <  2:1426437/144‑1 (MQ=255)
                                                                                                                   cATTTTTTCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAg                                         >  1:1426437/1‑144 (MQ=255)
                                                                                                                   aaTTTTTTCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTc                                   <  2:1146332/149‑1 (MQ=255)
                                                                                                                        tttCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAg                              <  1:433592/150‑1 (MQ=255)
                                                                                                                                aTCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATa                      <  2:567854/150‑1 (MQ=255)
                                                                                                                                      tgGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATg                                                                                                                            <  2:626701/42‑1 (MQ=255)
                                                                                                                                      tgGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATg                                                                                                                            >  1:626701/1‑42 (MQ=255)
                                                                                                                                        ggCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTg              <  2:710750/150‑1 (MQ=255)
                                                                                                                                          ccTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGcc            >  2:121861/1‑150 (MQ=255)
                                                                                                                                           cTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCg           >  2:1112670/1‑150 (MQ=255)
                                                                                                                                            tGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGc          >  2:362606/1‑150 (MQ=255)
                                                                                                                                              cTTCATCGGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCAc        <  1:803793/150‑1 (MQ=255)
                                                                                                                                              cTTCATCGGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCAc        <  1:803811/150‑1 (MQ=255)
                                                                                                                                                tcatcaGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCt      >  1:801150/1‑150 (MQ=255)
                                                                                                                                                   tcaGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCTATTGGGTCATGTGATGCGTCCCGGCACCGGGCCGGCCGCATTGTTATCCGCCGAGGTCGCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTgg   <  2:837056/150‑1 (MQ=255)
                                                                                                                                                   tcaGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTgg   >  1:1210738/1‑150 (MQ=255)
                                                                                                                                                    caGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGgc  <  1:265281/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TCTGGTCGATGCTGACGCGCTTGACGTCGTGACGACCGGTCATCACCACGCTCACCAGACCGCCACCCGACTGACCGGTGACTTCGGCGTTGGCCAGTTCTTCCTGCATCTTCTGCATTTTTTCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGC  >  NC_002947/4852832‑4853129

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