Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 252,258 A→G 23.1% N85S (AAT→AGT)  yafN → antitoxin of the YafO‑YafN toxin‑antitoxin system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913252,2580AG23.1% 181.9 / 60.7 104N85S (AAT→AGT) yafNantitoxin of the YafO‑YafN toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base A (39/41);  new base G (14/10);  total (53/51)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.05e-01

Only 100 of 105 total aligned reads displayed.

CTTGCTGAACAAGAGGAGAAAAAGCCCATAAAGGCGCGCTTCCGTCCAAGTGCTGCAAGATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCAGAGGAACTTGA  >  NC_000913/252161‑252356
                                                                                                 |                                                                                                  
cttgctGAACAAGAGGAGAAAAAGCCCATAAAGGCGCGCTTCCGTCCAAGTGCTGCAAGATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAATg                                                                                                  <  1:2009040/100‑1 (MQ=255)
 ttgctGAACAAGAGGAGAAAAAGCCCATAAAGGCGCGCTTCCGTCCAAGTGCTGCAAGATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAATGa                                                                                                 <  2:640209/100‑1 (MQ=255)
 ttgctGAACAAGAGGAGAAAAAGCCCATAAAGGCGCGCTTCCGTCCAAGTGCTGCAAGATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAATGa                                                                                                 <  1:84896/100‑1 (MQ=255)
   gctGAACAAGAGGAGAAAAAGCCCATAAAGGCGCGCTTCCGTCCAAGTGCTGCAAGATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAGTGata                                                                                               >  2:665388/1‑100 (MQ=255)
   gctGAACAAGAGGAGAAAAAGCCCATAAAGGCGCGCTTCCGTCCAAGTGCTGCAAGATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAATGata                                                                                               <  2:2056871/100‑1 (MQ=255)
    ctGAACAAGAGGAGAAAAAGCCCATAAAGGCGCGCTTCCGTCCAAGTGCTGCAAGATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAATGatat                                                                                              >  2:533408/1‑100 (MQ=255)
     tGAACAATAGGAGAAAAAGCCCATAAAGGCGCGCTTCCGTCAAAGTGCTGCAAGATTAGAGGAAATTACACGCCGCGCTGAGCAATATCTTAATGatat                                                                                              >  2:1055254/1‑99 (MQ=255)
     tGAACAAGAGGAGAAAAAGCCCATAAAGGCGCGCTTCCGTCCAAGTGCTGCAAGATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAGTGATATg                                                                                             >  2:2086203/1‑100 (MQ=255)
     tGAACAAGAGGAGAAAAAGCCCATAAAGGCGCGCTTCCGTCCAAGTGCTGCAAGATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAATGATATg                                                                                             >  1:2077449/1‑100 (MQ=255)
      gAACAAGAGGAGAAAAAGCCCATAAAGGCGCGCTTCCGTCCAAGTGCTGCAAGATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAATGATATGa                                                                                            >  2:232582/1‑100 (MQ=255)
      gAACAAGAGGAGAAAAAGCCCATAAAGGCGCGCTTCCGTCCAAGTGCTGCAAGATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAATGATATGa                                                                                            >  2:1847041/1‑100 (MQ=255)
      gAACAAGAGGAGAAAAAGCCCATAAAGGCGCGCTTCCGTCCAAGTGCTGCAAGATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAATGATATGa                                                                                            <  2:1235645/100‑1 (MQ=255)
      gAACAAGAGGAGAAAAAGCCCATAAAGGCGCGCTTCCGTCCAAGTGCTGCAAGATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAATGATATGa                                                                                            >  1:518189/1‑100 (MQ=255)
        aCAAGAGGAGAAAAAGCCCATAAAGGCGCTCTTCCGTCCAAGTGCTGCAAGATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAATGATATGACg                                                                                          <  1:193459/100‑1 (MQ=255)
        aCAAGAGGAGAAAAAGCCCATAAAGGCGCGCTTCCGTCCAAGTGCTGCAAGATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAATGATATGAc                                                                                           >  1:669081/1‑99 (MQ=255)
         cAAGAGGAGAAAAAGCCCATAAAGGCGCGCTTCCGTCCAAGTGCTGCAAGATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAATGATATGACg                                                                                          >  1:619006/1‑99 (MQ=255)
            gaggagAAAAAGCCCATAAAGGCGCGCTTCCGTCCAAGTGCTGCAAGATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAGTGATATGACGgatg                                                                                      <  2:1374187/100‑1 (MQ=255)
              ggagAAAAAGCCCATAAAGGCGCGCTTCCGTCCAAGTGCTGCAAGATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAGTGATATGACGgatgat                                                                                    <  1:2220072/100‑1 (MQ=255)
              ggagAAAAAGCCCATAAAGGCGCGCTTCCGTCCAAGTGCTGCAAGATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAATGATATGACGgatgat                                                                                    <  2:1055556/100‑1 (MQ=255)
                 gaAAAAGCCCATAAAGGTGCGCTTGCGTCCAGGTTGTGAAAGATTAGAGTAGATTACACGCCGCGCTGAACAATATCTTAATTATATGACGgatgatga                                                                                  <  2:900010/99‑1 (MQ=255)
                 gaAAAAGCCCATAAAGGCGCGCTTCCGTCCAAGTGCTGCAAGATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAGTGATATGACGgatgatgat                                                                                 >  1:1383208/1‑100 (MQ=255)
                 gaAAAAGCCCATAAAGGCGCGCTTCCGTCCAAGTGCTGCAAGATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAATGATATGACGgatgatgat                                                                                 <  1:556927/100‑1 (MQ=255)
                  aaaaaGCCCATAAAGGCGCGCTTCCGTCCAAGTGCTGCAAGATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAGTGATATGACGGATGATGAtt                                                                                <  1:1246086/100‑1 (MQ=255)
                  aaaaaGCCCATAAAGGCGCGCTTCCGTCCAAGTGCTGCAAGATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAGTGATATGACGGATGATGAtt                                                                                <  2:1635343/100‑1 (MQ=255)
                    aaaGCCCATAAAGGCGCGCTTCCGTCCAAGTGCTGCAAGATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAATGATATGACGGATGATGATTTc                                                                              <  2:1156789/100‑1 (MQ=255)
                      atcccATAAAGGCGCGCTTCCGTCCAAGTGCTACAAGATTAGAGGAAATTACACGACGCGCTGAACAATATCTTAATGAAATGACGGATGATGATTTCaa                                                                            <  2:521615/98‑1 (MQ=255)
                       gCCCATAAAGGCGCGCTTCCGTCCAAGTGCTGCAAGATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAGTGATATGACGGATGATGATTTCAAt                                                                           <  2:1427992/100‑1 (MQ=255)
                       gCCCATAAAGGCGCGCTTCCGTCCAAGTGCTGCAAGATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAGTGATATGACGGATGATGATTTCAAt                                                                           >  1:226183/1‑100 (MQ=255)
                         ccATAAAGGCGCGCTTCCGTCCAAGTGCTGCAAGATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAATGATATGACGGATGATGATTTCAATGa                                                                         >  2:1778742/1‑100 (MQ=255)
                           aTAAAGGCGCGCTTCCGTCCAAGTGCTGCAAGATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAGTGATATGACGGATGATGATTTCAATGACt                                                                       >  2:1905836/1‑100 (MQ=255)
                            tAAAGGCGCGCTTCCGTCCAAGTGCTGCAAGATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAGTGATATGACGGATGATGATTTCAATGACtt                                                                      >  2:1924121/1‑100 (MQ=255)
                            tAAAGGCGCGCTTCCGTCCAAGTGCTGCAAGATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAATGATATGACGGATGATGATTTCAATGACtt                                                                      <  2:461521/100‑1 (MQ=255)
                            tAAAGGCGCGCTTCCGTCCAAGTGCTGCAAGATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAATGATATGACGGATGATGATTTCAATGACtt                                                                      <  1:351586/100‑1 (MQ=255)
                              aaGGCGCGCTTCCGTCCAAGTGCTGCAAGATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTa                                                                    >  1:450054/1‑100 (MQ=255)
                                  cgcgcTTCCGTCCAAGTGCTGCAAGATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGa                                                                >  2:296328/1‑100 (MQ=255)
                                   gcgcTTCCGTCCAAGTGCTGCAAGATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGa                                                                <  2:1339724/99‑1 (MQ=255)
                                     gcTTCCGTCCAAGTGCTGCAAGATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATa                                                             >  2:1599897/1‑100 (MQ=255)
                                      cTTCCGTCCAAGTGCTGCAAGATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAGTGATATGACGGATGATGATTTCAATGACTTTAAGGAATaa                                                            >  2:1265252/1‑100 (MQ=255)
                                       ttCCGTCCAAGTGCTGCAAGATTAGAGGAAATTACACGCCGCGCTGAACAATACCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAg                                                           <  2:1239485/100‑1 (MQ=255)
                                        tCCGTCCAAGTGCTGCAAGATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAgg                                                          >  1:1823798/1‑100 (MQ=255)
                                          cGTCCAAGTGCTGCAAGATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGAt                                                        >  1:1392009/1‑100 (MQ=255)
                                          cGTCCAAGTGCTGCAAGATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGAt                                                        >  2:1762654/1‑100 (MQ=255)
                                           gTCCAAGTGCTGCAAGATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAGTGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGAt                                                        >  2:1831813/1‑99 (MQ=255)
                                            tCCAAGTGCTGCAAGATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAGTGATATGACGGATGATGATTTCAATGACTTTAAGCAATAAGGATGc                                                      <  2:1172356/100‑1 (MQ=255)
                                            tCCAAGTGCTGCAAGATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGc                                                      >  1:1219749/1‑100 (MQ=255)
                                             ccAAGTGCTGCAAGATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAATGATATGACGGATGATGATTTCACTGACTTTAAGGAATAAGGATGCg                                                     >  2:748314/1‑100 (MQ=255)
                                             ccAAGTGCTGCAAGATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCg                                                     <  1:1447548/100‑1 (MQ=255)
                                               aaGTGCTGCAAGATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCggg                                                   <  1:854668/100‑1 (MQ=255)
                                               aaGTGCTGCAAGATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCggg                                                   >  1:1504296/1‑100 (MQ=255)
                                               aaGTGCTGCAAGATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCggg                                                   <  2:834947/100‑1 (MQ=255)
                                                    ctgcAAGATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAGTGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTc                                              <  1:1265258/100‑1 (MQ=255)
                                                    ctgcAAGATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTAtt                                               <  1:1916370/99‑1 (MQ=255)
                                                    ctgcAAGATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTc                                              <  1:1896612/100‑1 (MQ=255)
                                                     tgcAAGATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTca                                             >  2:369499/1‑100 (MQ=255)
                                                     tgcAAGATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTca                                             >  1:726495/1‑100 (MQ=255)
                                                     tgcAAGATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTc                                              >  2:1609326/1‑99 (MQ=255)
                                                      gcAAGATTAGAGGAAATTACGCGCCGCGCTGAACAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTcaa                                            >  2:1368360/1‑100 (MQ=255)
                                                      gcAAGATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTcaa                                            >  2:505841/1‑100 (MQ=255)
                                                      gcAAGATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTcaa                                            <  2:1035827/100‑1 (MQ=255)
                                                      gcAAGAATAGATGAAATTACACGCCGCGCTGAACAATATCTGAATGATATGACGGATGATGATTTAAATGACTTTAAGGAATAAGGATGCGAGTATTaaa                                            <  2:1035717/100‑1 (MQ=255)
                                                       cAAGATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAGTGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTcaaa                                           >  2:22941/1‑100 (MQ=255)
                                                        ttccgatctggAAATTACACGCCGCGCTGAACAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTcaaaa                                          <  2:1967965/91‑1 (MQ=255)
                                                        aaGATTAGAGGAAATTACACGCCGCGCTGAACAATATATTAGTGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTcaaaa                                          >  2:1276284/1‑100 (MQ=255)
                                                           aTTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAGTGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTcaaaacaa                                       <  1:2086217/100‑1 (MQ=255)
                                                           aTTAGAAGAAATTACACGCCGCGCTGAACAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTcaaaacaa                                       <  2:1671930/100‑1 (MQ=255)
                                                            ttAGAGGAAATTACACGCCGCGCTGAACAATATCTTAGTGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTcaaaacaaa                                      >  2:761282/1‑100 (MQ=255)
                                                            ttAGAGGAAATTACACGCCGCGCTGAACAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTcaaaacaaa                                      <  1:1927478/100‑1 (MQ=255)
                                                              agagGAAATTACACGCCGCGCTGAACAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCaaaacaaaac                                    >  2:1412224/1‑100 (MQ=255)
                                                              agagGAAATTACACGCCGCGCTGAACAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCaaaacaaaac                                    <  1:877214/100‑1 (MQ=255)
                                                               gagTAAATTACACGCCGCGCTGAACAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACt                                   <  1:97365/100‑1 (MQ=255)
                                                               gagGAAATTACACGCCGCGCTGAACAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACt                                   >  1:615158/1‑100 (MQ=255)
                                                                 ggAAATTACACGCCGCGCTGAACAATATCTTAGAGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTa                                 >  2:1459207/1‑100 (MQ=255)
                                                                 ggAAATTACACGCCGCGCTGAACAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTcaaaaagatcggaa                                 >  1:1967978/1‑91 (MQ=255)
                                                                 ggAAATTACACGCCGCGCTGAACAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTa                                 >  1:1284667/1‑100 (MQ=255)
                                                                 ggAAATTACACGCCGCGCTGAACAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTa                                 >  2:854532/1‑100 (MQ=255)
                                                                    aaTTACACGCCGCGCTGAACAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTc                              <  1:895680/100‑1 (MQ=255)
                                                                     aTTACACGCCGCGCTGAACAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCg                             <  1:1362616/100‑1 (MQ=255)
                                                                       tACACGCCGCGCTGAACAATATCTTAGTGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCGcc                           <  1:2093999/100‑1 (MQ=255)
                                                                           cgccgcGCTGAACAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCa                       <  2:1902704/100‑1 (MQ=255)
                                                                           cgccgcGCTGAACAATATCTTAATGATATGAAAGATGATGATTTCAATGACTTTAAGGAATAAAGATGCGGGTATTCAAAACAAAACTTATTCGACTGCa                       >  2:2145293/1‑100 (MQ=255)
                                                                             ccgcGCTGAACAATATCTTAGTGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAAc                     >  1:761260/1‑100 (MQ=255)
                                                                              cgcgCTGAACAATATCTTAATGATATGACGGATGATGACTTCAATGACTTTAAGGAATTAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACt                    <  2:1219745/100‑1 (MQ=255)
                                                                                cgcTGAACAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTa                  <  1:1587166/100‑1 (MQ=255)
                                                                                 gcTGAACAATATCTTAGTGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTAc                 <  2:226183/100‑1 (MQ=255)
                                                                                 gcTGAACAATATCTTAGTGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTAc                 >  2:672424/1‑100 (MQ=255)
                                                                                   tGAACAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTAcag               <  1:448123/100‑1 (MQ=255)
                                                                                    gAACAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTAcagc              <  1:232582/100‑1 (MQ=255)
                                                                                     aaCAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTAcagca             >  1:192125/1‑100 (MQ=255)
                                                                                      aCAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAATCTTATTCGCCTGCAACTTAcagc              <  2:372138/98‑1 (MQ=255)
                                                                                       cAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCaga           >  1:2120166/1‑100 (MQ=255)
                                                                                       cAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCaga           >  1:156943/1‑100 (MQ=255)
                                                                                       cAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCaga           <  1:533407/100‑1 (MQ=255)
                                                                                        aaTATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCagag          <  2:503202/100‑1 (MQ=255)
                                                                                        aaTATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCagag          <  1:369498/100‑1 (MQ=255)
                                                                                        aaTATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCagag          <  1:444583/100‑1 (MQ=255)
                                                                                        aaTATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCagag          <  1:667621/100‑1 (MQ=255)
                                                                                         ataaCTTAATCATATAACGGAAGATGATTTCAATAACTTTAAGGAATAAGGAAGCAGGTATTTAAAACAAAACTTATTCACCTGACACTTACagaaga           >  2:1002781/1‑98 (MQ=255)
                                                                                               tAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCAGAGGAACTTg   >  1:181434/1‑100 (MQ=255)
                                                                                               tAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCGACTGCAACTTACAACAGAGGAACTTg   >  1:213317/1‑100 (MQ=255)
                                                                                                aaTGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCAGAGGAACTTGa  <  2:291188/100‑1 (MQ=255)
                                                                                                 |                                                                                                  
CTTGCTGAACAAGAGGAGAAAAAGCCCATAAAGGCGCGCTTCCGTCCAAGTGCTGCAAGATTAGAGGAAATTACACGCCGCGCTGAACAATATCTTAATGATATGACGGATGATGATTTCAATGACTTTAAGGAATAAGGATGCGGGTATTCAAAACAAAACTTATTCGCCTGCAACTTACAGCAGAGGAACTTGA  >  NC_000913/252161‑252356

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: