Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,212,088 (C)8→9 intergenic (‑85/+615) iraM ← / ← ycgX RpoS stabilzer during Mg starvation, anti‑RssB factor/DUF1398 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,212,0801.C98.4% 209.1 / ‑3.4 61intergenic (‑77/+623)iraM/ycgXRpoS stabilzer during Mg starvation, anti‑RssB factor/DUF1398 family protein
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (27/33);  total (27/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

ATATTACCCCATATGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTA‑‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTAT  >  NC_000913/1211930‑1212224
                                                                                                                                                       |                                                                                                                                                 
atatTACCCCATATGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTa                                                                                                                                                    <  2:532873/151‑1 (MQ=255)
   ttACCCCATATGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTA‑‑ccc                                                                                                                                               >  2:366323/1‑151 (MQ=255)
     aCCCCATATGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTA‑‑ccccc                                                                                                                                             <  1:110578/151‑1 (MQ=255)
         cATATGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTA‑‑ccccccccc                                                                                                                                         >  2:884073/1‑150 (MQ=255)
              ggCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTA‑‑ccccccccctcaa                                                                                                                                     <  2:450090/150‑6 (MQ=255)
                      aaGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCCCCCCCTCAAAAaaataa                                                                                                                              >  1:984268/1‑150 (MQ=255)
                         atatACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCCCCCCCTCAAAAAAATAaatta                                                                                                                          <  1:371429/151‑1 (MQ=255)
                           ataCCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCCCCCCCTCAAAAAAATAaattaa                                                                                                                         <  2:12848/150‑1 (MQ=255)
                            ttccACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCCCCCCCTCAAAAAAATAaattaat                                                                                                                        <  1:609675/148‑1 (MQ=255)
                             aCCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCCCCCCCTCAAAAAAATAaattaatt                                                                                                                       <  1:695102/150‑1 (MQ=255)
                             aCCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCCCCCCCTCAAAAAAATAaattaatt                                                                                                                       >  1:663770/1‑150 (MQ=255)
                                acaTGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCCCCCCCTCAAAAAAATAAATTAATTaaa                                                                                                                    >  1:700595/1‑150 (MQ=255)
                                 caTGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCCCCCCCTCAAAAAAATAAAttaattaa                                                                                                                     <  1:981569/148‑1 (MQ=255)
                                  aTGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCCCCCCCTCAAAAAAATAAATTAATTAAAatg                                                                                                                 >  2:3031/1‑151 (MQ=255)
                                   tGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCCCCCCCTCAAAAAAATAAATTAATTAAAatg                                                                                                                 >  1:915818/1‑150 (MQ=255)
                                   tGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCCCCCCCTCAAAAAAATAAATTAATTAAAatg                                                                                                                 >  1:741725/1‑150 (MQ=255)
                                     tAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCCCCCCCTCAAAAAAATAAATTAATTAAAat                                                                                                                  <  1:58292/147‑1 (MQ=255)
                                       ggTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCCCCCCCTCAAAAAAATAAATTAATTAAAatgatg                                                                                                              >  1:369624/1‑149 (MQ=255)
                                           gAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCCTCAAAAAAATAAATTTATTTAAatgaatggtt                                                                                                          >  2:843627/1‑144 (MQ=255)
                                           gAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTta                                                                                                         <  1:303816/150‑1 (MQ=255)
                                             aTTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTtatat                                                                                                      >  2:979502/1‑151 (MQ=255)
                                             aTTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTtatat                                                                                                      <  1:755840/151‑1 (MQ=255)
                                               tACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTtatata                                                                                                     >  1:379450/1‑150 (MQ=255)
                                                aCCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATataa                                                                                                    >  1:145323/1‑150 (MQ=255)
                                                 ccGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATataaaa                                                                                                  <  1:549597/151‑1 (MQ=255)
                                                       cAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATAtaaaataaaa                                                                                             <  1:657484/150‑1 (MQ=255)
                                                       cAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATaaaataaaat                                                                                            <  1:144621/151‑1 (MQ=255)
                                                              tATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCCCCCCCTCAAAAATATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTaa                                                                                        >  2:791655/1‑148 (MQ=255)
                                                              tATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGc                                                                                     >  1:616849/1‑151 (MQ=255)
                                                                 ccACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCaa                                                                                   <  1:329844/151‑1 (MQ=255)
                                                                  cACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAgg                                                                                 >  1:617022/1‑151 (MQ=255)
                                                                   aCTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGa                                                                                <  1:832510/151‑1 (MQ=255)
                                                                      tCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAAt                                                                              <  2:800538/150‑1 (MQ=255)
                                                                         tttGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAAtctc                                                                           >  1:368681/1‑150 (MQ=255)
                                                                          ttGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCa                                                                          <  1:1000232/150‑1 (MQ=255)
                                                                          ttGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCa                                                                          <  2:616595/150‑1 (MQ=255)
                                                                                  tCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGt                                                                  <  1:520894/150‑1 (MQ=255)
                                                                                     tATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTaa                                                               >  2:821909/1‑150 (MQ=255)
                                                                                     tATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTaa                                                               >  1:817515/1‑150 (MQ=255)
                                                                                       tCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAAc                                                             >  1:558262/1‑150 (MQ=255)
                                                                                       tCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAAc                                                             >  1:413832/1‑150 (MQ=255)
                                                                                       tCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAAc                                                             >  2:422620/1‑150 (MQ=255)
                                                                                             ggTATTTAATATGGTTCATTAGGATGTTTATTTATTGATTTTGCATATGAGTATATTAC‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATga                                                      <  2:65747/151‑1 (MQ=255)
                                                                                              gTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATg                                                       >  2:259373/1‑149 (MQ=255)
                                                                                                   tAATATGGTTCATTAGGATGTTTAATTCTTGATTTTGCATATGAGTATATTAC‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGAttt                                                 <  2:754746/150‑1 (MQ=255)
                                                                                                    aaTATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGAtttt                                                <  2:498032/150‑1 (MQ=255)
                                                                                                       atGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGa                                            <  2:547264/151‑1 (MQ=255)
                                                                                                        tGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTtgtg                                             <  2:741643/149‑1 (MQ=255)
                                                                                                         ggTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTA‑‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGaaa                                          <  1:75213/150‑1 (MQ=255)
                                                                                                          gTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAg                                         <  1:152484/151‑1 (MQ=255)
                                                                                                            tCATTAGGATGTTTATTTCTTTATTTTGCATATGAGTATATTAC‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGTATGTTAAACAAAATGAGATTTTGTGAAAg                                         <  2:567063/149‑1 (MQ=255)
                                                                                                             cATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCaa                                      <  2:511503/151‑1 (MQ=255)
                                                                                                             cATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCa                                       >  2:238830/1‑150 (MQ=255)
                                                                                                                      gTTTATTTCTTGATTTTGCATATGAGTATATTAC‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTCGTGAAAGCAATAAattat                              <  1:909816/150‑1 (MQ=255)
                                                                                                                          aTTTCTTGATTTTGCATATGAGTATATTAC‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGAc                          <  1:926499/150‑1 (MQ=255)
                                                                                                                            ttCTTGATTTTGCATATGAGTATATTAC‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACtt                        >  2:265349/1‑150 (MQ=255)
                                                                                                                              cTTGATTTTGCATATGAGTATATTAC‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGt                     >  1:847884/1‑151 (MQ=255)
                                                                                                                                 gATTTTGCATATGAGTATATTAC‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGttt                   <  1:286756/150‑1 (MQ=255)
                                                                                                                                  aTTTTGCATATGAGTATATTAC‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGtttt                  >  1:789272/1‑150 (MQ=255)
                                                                                                                                       gCATATGAGTATATTAC‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGAttt            <  1:41063/151‑1 (MQ=255)
                                                                                                                                        cATATGAGTATATTAC‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATATAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGAttt            >  1:538071/1‑150 (MQ=255)
                                                                                                                                        cATATGAGTATATTAC‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGAttt            >  1:118812/1‑150 (MQ=255)
                                                                                                                                             gAGTATATTAC‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTa       <  2:616960/150‑1 (MQ=255)
                                                                                                                                             gAGTATATTAC‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAg      <  1:718567/151‑1 (MQ=255)
                                                                                                                                                 atatTAC‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCtat  <  1:557460/151‑1 (MQ=255)
                                                                                                                                                 atatTAC‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCta   <  1:420541/150‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                 
ATATTACCCCATATGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTA‑‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTAT  >  NC_000913/1211930‑1212224

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