Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,212,088 (C)8→9 intergenic (‑85/+172) iraM ← / ← ymgK anti‑adaptor protein IraM, inhibitor of sigma(S) proteolysis/protein YmgK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,212,0801.C88.6% 1120.1 / 134.3 308intergenic (‑77/+180)iraM/ymgKanti‑adaptor protein IraM, inhibitor of sigma(S) proteolysis/protein YmgK
Reads supporting (aligned to +/- strand):  ref base . (17/18);  new base C (138/135);  total (155/153)
Fisher's exact test for biased strand distribution p-value = 8.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.12e-01

Only 100 of 359 total aligned reads displayed.

CGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTA‑‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTT  >  NC_000913/1211980‑1212177
                                                                                                     |                                                                                                  
cGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTa                                                                                                     >  1:1973107/1‑101 (MQ=255)
 gtgtCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTA‑‑c                                                                                                  <  1:1972292/101‑1 (MQ=255)
 gtgtCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTA‑‑c                                                                                                  <  1:1972964/101‑1 (MQ=255)
  tgtCAATTACTATCCACTACATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTA‑‑cc                                                                                                 <  1:1972605/101‑1 (MQ=255)
    tCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTA‑‑cccc                                                                                               >  2:1973112/1‑101 (MQ=255)
     cAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTA‑‑ccccc                                                                                              <  2:1972496/101‑1 (MQ=255)
     cAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTA‑‑ccccc                                                                                              <  1:1972463/101‑1 (MQ=255)
     cAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTA‑‑ccccc                                                                                              >  1:1973114/1‑101 (MQ=255)
           cTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTA‑‑ccccccccctc                                                                                        >  1:1973119/1‑98 (MQ=255)
            tATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTA‑‑ccccccccctca                                                                                       >  2:1973120/1‑97 (MQ=255)
            tATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTA‑‑ccccccccctca                                                                                       <  1:1972638/101‑5 (MQ=255)
              tCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCaaa                                                                                      <  2:1972426/101‑1 (MQ=255)
              tCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCaaa                                                                                      <  1:1972534/101‑1 (MQ=255)
               ccACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCaaaa                                                                                     <  1:1972658/101‑1 (MQ=255)
               ccACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCaaaa                                                                                     >  2:1973126/1‑101 (MQ=255)
                cACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTA‑‑CCCCCCCCTCaaaaaa                                                                                   <  2:1972560/101‑1 (MQ=255)
                 aCTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTA‑‑ccccccccctc                                                                                        >  1:1973129/1‑92 (MQ=255)
                 aCTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTA‑‑CCCCCCCCTCgaaaaaa                                                                                  <  1:1972446/101‑1 (MQ=255)
                       tttGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAaaataaat                                                                             <  1:1972376/101‑1 (MQ=255)
                       tttGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAaaataaat                                                                             <  2:1972697/101‑1 (MQ=255)
                       tttGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAaaataaat                                                                             <  2:1972590/101‑1 (MQ=255)
                       tttGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTA‑‑CCCCCCCCTCAAAAAAATAaatt                                                                            >  2:1973135/1‑101 (MQ=255)
                        ttGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAaaat                                                                                 <  1:1972515/96‑1 (MQ=255)
                          gTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAaattaa                                                                          >  1:1973139/1‑101 (MQ=255)
                           ttATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAaattaat                                                                         >  2:1973142/1‑101 (MQ=255)
                            tATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTA‑‑CCCCCCCCTCAAAAAAATAAattaatta                                                                       <  1:1972629/101‑1 (MQ=255)
                             aTGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAattaatta                                                                       <  2:1972456/101‑1 (MQ=255)
                              tGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAAttaattaa                                                                      >  2:1973147/1‑101 (MQ=255)
                              tGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAAttaattaa                                                                      <  1:1972618/101‑1 (MQ=255)
                              tGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAAttaattaa                                                                      <  2:1972611/101‑1 (MQ=255)
                               gTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTaaa                                                                     >  1:1973149/1‑101 (MQ=255)
                                 cTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAat                                                                   <  2:1972626/101‑1 (MQ=255)
                                 cTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAat                                                                   >  2:1973151/1‑101 (MQ=255)
                                  ttATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAatg                                                                  <  1:1972394/101‑1 (MQ=255)
                                    aTCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAatgat                                                                <  1:1972520/101‑1 (MQ=255)
                                    aTCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAatgat                                                                >  2:1973158/1‑101 (MQ=255)
                                     tCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAatgatg                                                               <  1:1972886/101‑1 (MQ=255)
                                     tCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAatgatg                                                               <  1:1972225/101‑1 (MQ=255)
                                      cccACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTA‑‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGc                                                             <  1:1972499/101‑1 (MQ=255)
                                       ccACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGc                                                             >  2:1973161/1‑101 (MQ=255)
                                       ccACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGc                                                             <  2:1972751/101‑1 (MQ=255)
                                        cACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCt                                                            >  1:1973164/1‑101 (MQ=255)
                                         aCGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCtt                                                           >  2:1973166/1‑101 (MQ=255)
                                         aCGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTA‑‑ccccccccct                                                                                         >  1:1973169/1‑68 (MQ=255)
                                             tATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTtatat                                                       >  1:1973174/1‑101 (MQ=255)
                                              aTTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTtatata                                                      >  1:1973180/1‑101 (MQ=255)
                                              aTTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTtatata                                                      >  2:1973177/1‑101 (MQ=255)
                                                 tAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATataaaa                                                   <  2:1972734/101‑1 (MQ=255)
                                                 tAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATataaaa                                                   >  1:1973186/1‑101 (MQ=255)
                                                  aaTATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTA‑‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATataaaata                                                 >  2:1973187/1‑101 (MQ=255)
                                                    tatGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATataaaataa                                                <  2:1972299/101‑1 (MQ=255)
                                                    tatGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATataaaataa                                                >  1:1973193/1‑101 (MQ=255)
                                                     atGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATataaaataaa                                               >  2:1973198/1‑101 (MQ=255)
                                                     atGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATataaaataaa                                               >  1:1973197/1‑101 (MQ=255)
                                                       ggTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATaaaataaaat                                             <  1:1972681/101‑1 (MQ=255)
                                                       ggTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATaaaataaaat                                             <  2:1972577/101‑1 (MQ=255)
                                                        gTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAAtt                                            <  2:1972983/101‑1 (MQ=255)
                                                         ttCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATAtaaaataaaa                                              >  2:1973201/1‑98 (MQ=255)
                                                          tCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTa                                          <  2:1972847/101‑1 (MQ=255)
                                                           cATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATgg                                                              >  1:1973207/1‑80 (MQ=255)
                                                            aTTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATataaaataa                                                <  2:1972547/93‑1 (MQ=255)
                                                            aTTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTaaa                                        >  2:1973209/1‑101 (MQ=255)
                                                            aTTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTaaa                                        <  2:1972519/101‑1 (MQ=255)
                                                             ttAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAg                                       <  2:1972662/101‑1 (MQ=255)
                                                              tAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGc                                      <  2:1972372/101‑1 (MQ=255)
                                                              tAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGc                                      <  2:1972708/101‑1 (MQ=255)
                                                              tAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGc                                      >  2:1973215/1‑101 (MQ=255)
                                                               aGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAatga                                                                 >  2:1973220/1‑73 (MQ=255)
                                                               aGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCa                                     <  1:1972782/101‑1 (MQ=255)
                                                               aGGATGTTTATTTCTTGATTTTGCATATGAGTATATTA‑‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCaa                                    <  1:1972564/101‑1 (MQ=255)
                                                               aGGATGTTTATTTCTTGATTTTGCATATGAGTATATTA‑‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCaa                                    >  1:1973217/1‑101 (MQ=255)
                                                                ggATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCaa                                    <  1:1972819/101‑1 (MQ=255)
                                                                  aTGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAgg                                  >  1:1973224/1‑101 (MQ=255)
                                                                      ttATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAAtc                              >  2:1973228/1‑101 (MQ=255)
                                                                       tATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAAtct                             <  1:1972762/101‑1 (MQ=255)
                                                                        aTTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAAtctc                            >  1:1973232/1‑101 (MQ=255)
                                                                         tttCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCa                           >  2:1973236/1‑101 (MQ=255)
                                                                          ttCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCaa                          <  2:1973052/101‑1 (MQ=255)
                                                                          ttCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCaa                          >  2:1973240/1‑101 (MQ=255)
                                                                          ttCTTGATTTTGCATATGAGTATATTA‑‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAAt                         <  1:1972792/101‑1 (MQ=255)
                                                                           tCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAAt                         <  1:1972726/101‑1 (MQ=255)
                                                                             ttGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATgg                       <  1:1972704/101‑1 (MQ=255)
                                                                              tGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGa                      >  2:1973245/1‑101 (MQ=255)
                                                                               gATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGAt                     <  2:1972597/101‑1 (MQ=255)
                                                                               gATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGAt                     >  1:1973246/1‑101 (MQ=255)
                                                                                aTTTTGCATATGAGTATATTACCCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAACAAGGAATCTCAATGGAt                     >  1:1973252/1‑101 (MQ=255)
                                                                                aTTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAAAAAaataatta                                                                       >  2:1973253/1‑50 (MQ=255)
                                                                                  tttGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGtt                  >  1:1973259/1‑101 (MQ=255)
                                                                                    tGCATATGAGTATATTA‑‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTaaa               <  1:1972575/101‑1 (MQ=255)
                                                                                      cATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAAc              <  2:1972785/101‑1 (MQ=255)
                                                                                       atatGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACa             >  2:1973262/1‑101 (MQ=255)
                                                                                       atatGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACa             >  1:1973264/1‑101 (MQ=255)
                                                                                         atGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACaaa           <  2:1972675/101‑1 (MQ=255)
                                                                                         atGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACaaa           >  2:1973268/1‑101 (MQ=255)
                                                                                          tGAGTATATTACCCCCCCCCTCAAAAAAAAAAATTAATTAAAatga                                                                 >  1:1973271/1‑46 (MQ=255)
                                                                                            aGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATg        <  1:1972742/101‑1 (MQ=255)
                                                                                              tataTTA‑‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATgaga     >  2:1973276/1‑101 (MQ=255)
                                                                                               atatTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATgaga     <  2:1972832/101‑1 (MQ=255)
                                                                                                tatTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGAt    >  1:1973279/1‑101 (MQ=255)
                                                                                                tatTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGAt    <  1:1972729/101‑1 (MQ=255)
                                                                                                     |                                                                                                  
CGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTA‑‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTT  >  NC_000913/1211980‑1212177

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