Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,642,793 G→T 5.4% R18S (CGC→AGC)  rodZ ← MreB assembly cytoskeletal protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,642,7930GT 5.4% 280.2 / 3.7 83R18S (CGC→AGC) rodZMreB assembly cytoskeletal protein
Reads supporting (aligned to +/- strand):  ref base G (37/41);  new base T (2/3);  total (39/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.26e-01

GATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGTACGTTTAAATGAGAATTGAGATGCCGGTGAACCATCATGCTCACCCACACCGCGAAACATCCGTTAAGTTAACCGTTATCATACAG  >  NC_000913/2642650‑2642941
                                                                                                                                               |                                                                                                                                                    
gATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGa                                                                                                                                                     >  8:84365/1‑145 (MQ=255)
gATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGa                                                                                                                                                     <  7:84365/145‑1 (MQ=255)
gATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAgcgc                                                                                                                                                 <  5:356277/149‑1 (MQ=255)
 atatatCCGCGCAGGAATGGTGAGGCAAGCTCGGCGGGTGCCTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCACCGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGcc                                                                                                                                                <  7:295676/149‑1 (MQ=255)
    tatCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCgg                                                                                                                                              <  8:23772/148‑1 (MQ=255)
      tCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTATTTGTTCGCGAGCATTACGCAGGCGAGCGCCGgta                                                                                                                                            <  2:152793/148‑1 (MQ=255)
      tCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGgtag                                                                                                                                           >  4:79727/1‑149 (MQ=255)
       ccGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGgtag                                                                                                                                           >  1:95172/1‑148 (MQ=255)
        cgcgCAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGgtagt                                                                                                                                          <  3:79727/148‑1 (MQ=255)
             aGGAATGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTCCGCAGGCGAGCGCCGGTAGTAAGTGc                                                                                                                                    >  1:63547/1‑149 (MQ=255)
             aGGAATGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGc                                                                                                                                    <  8:132360/149‑1 (MQ=255)
             aGGAATGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGc                                                                                                                                    <  7:268938/149‑1 (MQ=255)
             aGGAATGTTGAAGAAAGATCTGCGGGTGCCTTATCTTCTTCAATGTAGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGc                                                                                                                                    <  6:188251/149‑1 (MQ=255)
                       aaGCAAGATCGGCGGGTGCCTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTg                                                                                                                           <  1:237259/148‑1 (MQ=255)
                       aaGCAAGATCGGCGGGTGCCTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCAGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTgg                                                                                                                          >  1:139107/1‑149 (MQ=255)
                         gCAAGATCGGCGGGTGCCTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTc                                                                                                                        >  5:361310/1‑149 (MQ=255)
                          cAAGATCGGCGGGTGCCTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTcg                                                                                                                       >  1:239920/1‑149 (MQ=255)
                                         ccTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGgtagt                                                                                                                                          <  2:314800/115‑1 (MQ=255)
                                         ccTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGgtagt                                                                                                                                          >  1:314800/1‑115 (MQ=255)
                                          cTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTgg                                                                                                                          <  7:123441/130‑1 (MQ=255)
                                          cTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTgg                                                                                                                          >  8:123441/1‑130 (MQ=255)
                                                 gcttcAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTat                                                                                                       <  5:104110/141‑1 (MQ=255)
                                                  cttcAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTattcattc                                                                                                 <  7:226654/147‑1 (MQ=255)
                                                   ttcAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCt                                                                                              >  5:157626/1‑149 (MQ=255)
                                                     cAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGc                                                                                               <  2:95172/146‑1 (MQ=255)
                                                     cAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTg                                                                                             >  7:288605/1‑148 (MQ=255)
                                                     cAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGGCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTg                                                                                             <  5:180956/148‑1 (MQ=255)
                                                        tGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCTAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAgg                                                                                         >  3:368399/1‑149 (MQ=255)
                                                        tGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCTAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAgg                                                                                         >  3:362661/1‑149 (MQ=255)
                                                            gcgTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTg                                                                                                                                     <  1:204630/101‑1 (MQ=255)
                                                            gcgTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTg                                                                                                                                     >  2:204630/1‑101 (MQ=255)
                                                                  cGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGTACGTTTaaa                                                                               >  2:326635/1‑149 (MQ=255)
                                                                    tGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGTACGTTTAAATg                                                                             <  5:197008/149‑1 (MQ=255)
                                                                    tGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATGTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGTACGTTTAAAt                                                                              >  8:96335/1‑148 (MQ=255)
                                                                     ggAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGTACGTTTAAAt                                                                              >  4:257382/1‑147 (MQ=255)
                                                                      gAATCCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGTTCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGTACGTTTTAATga                                                                            <  4:295249/148‑1 (MQ=255)
                                                                       aaaCCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGTACGTTTAAATgag                                                                           >  3:379071/1‑148 (MQ=255)
                                                                        aaCCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGTACGTTTAAATGAGaa                                                                         <  3:355179/149‑1 (MQ=255)
                                                                          ccTTCAGGCAAAGGCGCTCGGCATCGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGTACGTTTAAATGAGAAtt                                                                       <  2:367701/149‑1 (MQ=255)
                                                                           cTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCgg                                                                                                                                              <  6:310497/77‑1 (MQ=255)
                                                                           cTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCgg                                                                                                                                              >  5:310497/1‑77 (MQ=255)
                                                                            ttCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGTATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGTACGTTTAAATGAGAATTga                                                                     >  6:85499/1‑149 (MQ=255)
                                                                            ttCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGTACGTTTAAATGAGAAtt                                                                       >  7:320426/1‑147 (MQ=255)
                                                                              cAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTGGTGCGTGGCTTCAGTATTCATTCGCTGCAGGTACGTTTAAATGAGAATTgag                                                                    <  5:99114/148‑1 (MQ=255)
                                                                               aGGCAAAGTCGCTAGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGTACGTTTAAATGAGAATTgaga                                                                   <  8:167267/148‑1 (MQ=255)
                                                                                 gCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGCGCGCCGGTAGTAAGTGCTTCATTTAGGTCGTGCGTGGCTTCAGTATTCATTCTCTGCAGGTACGTGTAAATGAGAATTGAGATGc                                                                >  8:254402/1‑149 (MQ=255)
                                                                                 gCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTat                                                                                                       <  7:82670/110‑1 (MQ=255)
                                                                                 gCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTat                                                                                                       >  8:82670/1‑110 (MQ=255)
                                                                                   aaaGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGTACGTTTAAATGAGAATTGAGATGCCg                                                              >  3:7959/1‑149 (MQ=255)
                                                                                      ttCGCTCGGATACGGCCTTCTGAGTAAGTGCTAGTTTTTCGCGAGCTTTATGCTGGCTAGCGCCGGTATTAAGTGCTTCATTTTGGTCGTTGGTGGTTTCAGTATTCATTCGCTGCAGGTTCGTTTAAATGAGAATTTAGATGACGGt                                                            <  8:320426/147‑1 (MQ=255)
                                                                                      gTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGTACGTTTAAATGAGAATTGAGATGCCGGt                                                            >  8:399401/1‑148 (MQ=255)
                                                                                      gTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGTACGTTTAAATGAGAATTGAGATGCCGGt                                                            >  7:146495/1‑148 (MQ=255)
                                                                                           tCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGTACGTTTAAATGAGAATTGAGATGCCGGTGAACca                                                      >  6:134964/1‑149 (MQ=255)
                                                                                             ggCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGTACGTTTAAATGAGAATTGAGATGCCGGTGAACcatc                                                    <  2:253896/149‑1 (MQ=255)
                                                                                             ggCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGTACGTTTAAATGAGAATTGAGATGCCGGTGAACcatc                                                    <  2:323662/149‑1 (MQ=255)
                                                                                              gCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCTAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGTACGTTTAAATGAGAATTGAGATGCCGGTGAACcatca                                                   <  4:368399/149‑1 (MQ=255)
                                                                                              gCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCTAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGTACGTTTAAATGAGAATTGAGATGCCGGTGAACcatca                                                   <  4:362661/149‑1 (MQ=255)
                                                                                                aaCGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTATGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGTACGTTTAAATGATAATTGAGATGCCGGTGAACCATCATg                                                 <  2:139107/149‑1 (MQ=255)
                                                                                                    gCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTACATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGTACGTTTAAATGAGAATTGAGATGCCGGTGACCAATCATGCTCa                                             >  6:95797/1‑149 (MQ=255)
                                                                                                          tgACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGAAGGTACGTTTAAATGAGAATTGAGATGCCGGTGAACCATCATGCTCACCccca                                       >  6:337056/1‑149 (MQ=255)
                                                                                                              tAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTc                                                                                                                        <  4:180567/64‑1 (MQ=255)
                                                                                                              tAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTc                                                                                                                        >  3:180567/1‑64 (MQ=255)
                                                                                                              tAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGTACGTTTAAATGAGAATTGAGATGCCGGTGAACCATCATGCTCACCCACACcgc                                   >  3:360296/1‑149 (MQ=255)
                                                                                                              tAAGACCTAGTTGTTCGCGAGCATAACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGTACGTTTAAATGAGAATTGAGATGCCGGAGAACCATCATGCTCACCCCCACcgc                                   >  8:126827/1‑149 (MQ=255)
                                                                                                                aGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGTACGTTTAAATGAGAATTGAGATGCCGGTGAACCATCATGCTCACCCACACcgcg                                  <  4:7959/148‑1 (MQ=255)
                                                                                                                 ggccTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTTGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGTACGTTTAAATGAGAATTTAGATGCCGGTGAAGCATCATGCTCACCCACACCGCGaa                                <  5:95797/147‑1 (MQ=255)
                                                                                                                 gTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGTACGTTTAAATGAGAATTGAGATGCCGGTGAACCATCATGCTCACCCACACCGCGa                                 >  3:176065/1‑148 (MQ=255)
                                                                                                                  tCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGTa                                                                                       >  1:250219/1‑93 (MQ=255)
                                                                                                                  tCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGTa                                                                                       <  2:250219/93‑1 (MQ=255)
                                                                                                                   ccTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGTACGTTTAAATGAGAATTGAGATGCCGGTGAACCATCATGCTCACCCACACCGCGaaa                               >  2:103090/1‑148 (MQ=255)
                                                                                                                    cTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGTACGTTTAAATga                                                                            <  6:381345/102‑1 (MQ=255)
                                                                                                                    cTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGTACGTTTAAATga                                                                            >  5:381345/1‑102 (MQ=255)
                                                                                                                    cTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGGCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGTACGTTTAAATGAGAATTGAGATGCCGGTGAACCATCATGCTCACCCACACCGCGAAAc                              <  7:399401/148‑1 (MQ=255)
                                                                                                                          gttCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGTACGTTTAAATGAGAATTGAGATGCCGGTGAACCATCATGCTCACCCACACCGCGAAACATCCgt                        <  3:257382/148‑1 (MQ=255)
                                                                                                                           ttCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGTACGTTTAAATGAGAATTGAGATGCCGGTGAACCATCATGCTCACCCACACCGCGAAACATCCgtt                       <  2:239920/148‑1 (MQ=255)
                                                                                                                              gcgAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGTACGTTTAAATGAGAATTGAGATGCCGGTGAACCATCATGCTCACCCACACCGCGAAACATCCgttaagt                   >  2:70303/1‑149 (MQ=255)
                                                                                                                                 aGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGTACGTTTAAATGAGAATTGAGATGCCGGTGAACCATCATGCTCACCCACACCGCGAAACATCCgttaagttaa                <  1:326635/149‑1 (MQ=255)
                                                                                                                                         gCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGTACGTTTAAATGAGAATTGAGATGCCGGTGAACCATCATGCTCACCCACACCGCGAAACATCCgttaa                     <  8:245446/136‑1 (MQ=255)
                                                                                                                                         gCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGTACGTTTAAATGAGAATTGAGATGCCGGTGAACCATCATGCTCACCCACACCGCGAAACATCCgttaa                     >  7:245446/1‑136 (MQ=255)
                                                                                                                                            ggCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGTACGTTTAAATGAGAATTGAGATGCCGGTGAACCATCATGCTCACCCACACCGCGAAACATCCGTTAAGTTAACCGTTATCATa     <  6:361310/149‑1 (MQ=255)
                                                                                                                                              cGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGTACGTTTAAATGAGAATTGAGATGCCGGTGAACCATCATGCTCACCCACACCGCGAAACATCCGTTAAGTTAACCGTTATCATa     >  4:161253/1‑147 (MQ=255)
                                                                                                                                              cGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGTACGTTTAAATGAGAATTGAGATGCCGGTGAACCATCATGCTCACCCACACCGCGAAACATCCGTTAAGTTAACCGTTATCATACa   <  5:85499/149‑1 (MQ=255)
                                                                                                                                              cGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTGCAGTATTCATTCGCTGCAGGTACGTTTAAATGAGAATTGAGATGCCGGTGAACCATCATGCTCACCCACACCGCGAAACATCCGTTAAGTTAACCGTTATCATACa   <  2:63547/149‑1 (MQ=255)
                                                                                                                                               gAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGTACGTTTAAATGAGAATTGAGATGCCGGTGAACCATCATGCTCACCCACACCGCGAAACATCCGTTAAGTTAACCGTTATCATACAg  >  3:101682/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
GATATATCCGCGCAGGAATGTTGAAGCAAGATCGGCGGGTGCCTTATCTTCTTCAATGTCGCGTACCGTGGAAACCTTCAGGCAAAGTCGCTCGGCAACGGCCTGCTGACTAAGTCCTAGTTGTTCGCGAGCATTACGCAGGCGAGCGCCGGTAGTAAGTGCTTCATTTTGGTCGTGCGTGGCTTCAGTATTCATTCGCTGCAGGTACGTTTAAATGAGAATTGAGATGCCGGTGAACCATCATGCTCACCCACACCGCGAAACATCCGTTAAGTTAACCGTTATCATACAG  >  NC_000913/2642650‑2642941

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