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.C97.4% 1118.0 / 16.6 271intergenic (‑77/+180)iraM/ymgKanti‑adaptor protein IraM, inhibitor of sigma(S) proteolysis/protein YmgK
Reads supporting (aligned to +/- strand):  ref base . (1/6);  new base C (129/135);  total (130/141)
Fisher's exact test for biased strand distribution p-value = 1.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 321 total aligned reads displayed.

CGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGAT  >  NC_000913/1211980‑1212175
                                                                                                      |                                                                                               
  tgtCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAcc                                                                                               >  2:1783764/1‑101 (MQ=255)
  tgtCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAcc                                                                                               <  1:1783359/101‑1 (MQ=255)
    tCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAcccc                                                                                             <  1:1782918/101‑1 (MQ=255)
    tCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGCATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAcccc                                                                                             >  1:1783768/1‑101 (MQ=255)
     cAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAccccc                                                                                            >  2:1783769/1‑101 (MQ=255)
       aTTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAccccccc                                                                                          <  1:1783259/101‑1 (MQ=255)
         tACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAccccccccc                                                                                        >  2:1783774/1‑100 (MQ=255)
          aCTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAccccccccct                                                                                       <  2:1783190/101‑3 (MQ=255)
           cTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAccccccccctc                                                                                      <  2:1783322/101‑4 (MQ=255)
           cTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAccccccccctc                                                                                      >  1:1783777/1‑98 (MQ=255)
           cTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAccccccccct                                                                                       <  2:1783368/100‑3 (MQ=255)
            tATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAccccccccctca                                                                                     <  2:1783434/101‑5 (MQ=255)
            tATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAccccccccctca                                                                                     <  2:1783172/101‑5 (MQ=255)
            tATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAccccccccctca                                                                                     >  2:1783781/1‑97 (MQ=255)
             aTCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAccccccccctcaa                                                                                    >  1:1783782/1‑96 (MQ=255)
               ccACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCaaaa                                                                                   <  1:1783380/101‑1 (MQ=255)
                cACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCaaaaa                                                                                  <  1:1783270/101‑1 (MQ=255)
                cACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCaaaaa                                                                                  >  2:1783786/1‑101 (MQ=255)
                  cTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCaaaaaaa                                                                                >  1:1783789/1‑101 (MQ=255)
                     cATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAaaataa                                                                             >  2:1783794/1‑101 (MQ=255)
                      aTTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAaaataaa                                                                            <  2:1783089/101‑1 (MQ=255)
                       tttGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTCCCCCCCCCCTCAAAAaaataaat                                                                           <  1:1783285/101‑1 (MQ=255)
                        ttGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAaatt                                                                          <  2:1783141/101‑1 (MQ=255)
                        ttGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAaatt                                                                          <  2:1783220/101‑1 (MQ=255)
                        ttGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAaatt                                                                          >  2:1783799/1‑101 (MQ=255)
                          gTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAaattaa                                                                        >  2:1783802/1‑101 (MQ=255)
                            tATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAaattaatt                                                                      >  2:1783805/1‑101 (MQ=255)
                             aTGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAccccccccctc                                                                                      >  2:1783808/1‑80 (MQ=255)
                              tGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAAttaattaa                                                                    >  1:1783809/1‑101 (MQ=255)
                               gTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTaaa                                                                   <  1:1783115/101‑1 (MQ=255)
                                tCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTaaaa                                                                  >  1:1783812/1‑101 (MQ=255)
                                 cTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAat                                                                 <  1:1783219/101‑1 (MQ=255)
                                 cTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAat                                                                 <  1:1783192/101‑1 (MQ=255)
                                  ttATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAatg                                                                >  2:1783816/1‑101 (MQ=255)
                                  ttATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAatg                                                                <  1:1783147/101‑1 (MQ=255)
                                    aTCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAatgat                                                              >  2:1783822/1‑101 (MQ=255)
                                    aTCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAatgat                                                              >  1:1783821/1‑101 (MQ=255)
                                     tCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTaaaa                                                                  >  1:1783826/1‑96 (MQ=255)
                                     tCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAatgatg                                                             <  1:1783331/101‑1 (MQ=255)
                                       ccACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGc                                                           >  1:1783831/1‑101 (MQ=255)
                                       ccACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCt                                                          >  2:1783830/1‑101 (MQ=255)
                                        cACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCt                                                          <  2:1783330/101‑1 (MQ=255)
                                         aCGGTATTTTATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATCAATTAATTAAAATGATGGCtt                                                         <  1:1783588/101‑1 (MQ=255)
                                          cGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTta                                                        >  2:1783834/1‑101 (MQ=255)
                                            gTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTtata                                                      <  2:1783237/101‑1 (MQ=255)
                                            gTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTtata                                                      <  1:1783133/101‑1 (MQ=255)
                                              aTTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTtatata                                                    <  1:1783441/101‑1 (MQ=255)
                                               tttAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATataa                                                   >  1:1783839/1‑101 (MQ=255)
                                                ttAATATGGGTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATataaa                                                  <  1:1783206/101‑1 (MQ=255)
                                                 tAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAatga                                                               >  2:1783846/1‑87 (MQ=255)
                                                 tAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATataaaa                                                 <  2:1783421/101‑1 (MQ=255)
                                                 tAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATataaaa                                                 >  2:1783843/1‑101 (MQ=255)
                                                  aaTATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATataaaat                                                >  1:1783849/1‑101 (MQ=255)
                                                    tatGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATataaaataa                                              >  1:1783853/1‑101 (MQ=255)
                                                    tatGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATataaaataa                                              <  2:1783361/101‑1 (MQ=255)
                                                      tGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATAtaaaataaaa                                            <  1:1783501/101‑1 (MQ=255)
                                                       ggTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATaaaataaaat                                           <  1:1783393/101‑1 (MQ=255)
                                                       ggTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATaaaataaaat                                           <  2:1783497/101‑1 (MQ=255)
                                                        gTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAAtt                                          <  1:1783348/101‑1 (MQ=255)
                                                        gTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAAtt                                          <  2:1783260/101‑1 (MQ=255)
                                                        gTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAAtt                                          >  1:1783856/1‑101 (MQ=255)
                                                         ttCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAAttt                                         <  1:1783229/101‑1 (MQ=255)
                                                            aTTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTaaa                                      >  2:1783861/1‑101 (MQ=255)
                                                             ttAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAg                                     <  1:1783371/101‑1 (MQ=255)
                                                             ttAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAg                                     >  1:1783864/1‑101 (MQ=255)
                                                                ggATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCaa                                  <  2:1783403/101‑1 (MQ=255)
                                                                 gATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAg                                 >  2:1783870/1‑101 (MQ=255)
                                                                 gATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCCCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAg                                 >  1:1783869/1‑101 (MQ=255)
                                                                    gTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGaa                              >  2:1783875/1‑101 (MQ=255)
                                                                    gTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGaa                              <  1:1783045/101‑1 (MQ=255)
                                                                     tttATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAAt                             <  2:1783723/101‑1 (MQ=255)
                                                                     tttATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAAt                             >  1:1783876/1‑101 (MQ=255)
                                                                      ttATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAAt                             >  2:1783881/1‑100 (MQ=255)
                                                                       tATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAAtct                           <  2:1783202/101‑1 (MQ=255)
                                                                       tATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAAtct                           >  1:1783882/1‑101 (MQ=255)
                                                                        aTTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAAtctc                          <  2:1783308/101‑1 (MQ=255)
                                                                            cTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATg                      >  2:1783887/1‑101 (MQ=255)
                                                                             ttGATTTTGCATATGAGTATATTCCCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATataaaata                                               <  2:1783041/75‑1 (MQ=255)
                                                                             ttGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATgg                     <  1:1783559/101‑1 (MQ=255)
                                                                              tGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGa                    >  1:1783888/1‑101 (MQ=255)
                                                                                 ttttGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGt                 <  1:1783161/101‑1 (MQ=255)
                                                                                  tttGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGtt                <  2:1783280/101‑1 (MQ=255)
                                                                                  tttGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGtt                >  2:1783894/1‑101 (MQ=255)
                                                                                   ttGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTa               <  1:1783454/101‑1 (MQ=255)
                                                                                    tGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTaa              <  2:1783295/101‑1 (MQ=255)
                                                                                     gCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTaaa             <  1:1783405/101‑1 (MQ=255)
                                                                                      cATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAAc            >  2:1783898/1‑101 (MQ=255)
                                                                                      cATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAAc            >  1:1783900/1‑101 (MQ=255)
                                                                                       atatGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACa           >  2:1783904/1‑101 (MQ=255)
                                                                                        tatGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACaa          >  1:1783905/1‑101 (MQ=255)
                                                                                         atGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACaaa         <  1:1783634/101‑1 (MQ=255)
                                                                                         atGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACaaa         <  1:1783473/101‑1 (MQ=255)
                                                                                          tGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAatgat                                                              >  2:1783911/1‑47 (MQ=39)
                                                                                          tGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACaaaa        >  1:1783909/1‑101 (MQ=255)
                                                                                          tGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACaaaa        <  1:1783424/101‑1 (MQ=255)
                                                                                              tataTTACCCCCCCCCTCAAAAAAATAAATTAATTAAAatgatg                                                             <  1:1783312/44‑1 (MQ=38)
                                                                                              tataTTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATgag    <  1:1783520/101‑1 (MQ=255)
                                                                                               atatTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATgaga   >  1:1783917/1‑101 (MQ=255)
                                                                                                tatTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATg                      <  2:1783247/81‑1 (MQ=255)
                                                                                                tatTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGAt  <  2:1783586/101‑1 (MQ=255)
                                                                                                      |                                                                                               
CGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGAT  >  NC_000913/1211980‑1212175

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: