Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,050,225 A→T 7.7% L181M (TTG→ATG)  yihA ← cell division GTP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,050,2250AT 7.7% 188.0 / 3.3 70L181M (TTG→ATG) yihAcell division GTP‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (30/34);  new base T (4/2);  total (34/36)
Fisher's exact test for biased strand distribution p-value = 4.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGAACGGCGATATTGCTGTCTACCGCCCAC  >  NC_000913/4050077‑4050373
                                                                                                                                                    |                                                                                                                                                    
tATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCAATAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCaa                                                                                                                                                      <  1:65982/149‑1 (MQ=255)
     aaTTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACgaag                                                                                                                                                 <  1:789991/149‑1 (MQ=255)
          aCGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGAATCTCGCTAAACCAGGTATCCCGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAAc                                                                                                                                            <  2:686791/149‑1 (MQ=255)
           cGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCGTTGAAGaaacta                                                                                                                                           <  1:17752/149‑4 (MQ=255)
           cGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGATTCTTCAACGAAGAAACCg                                                                                                                                           >  1:305217/1‑149 (MQ=255)
            tggCCGGATACGCCACATCCTGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTATTCAACGAAGAAAACGt                                                                                                                                          <  1:630162/148‑1 (MQ=255)
             ggCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGGTTCTTCAACGAAGAAAACGtt                                                                                                                                         <  1:396911/149‑1 (MQ=255)
                  taTACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTTCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCAGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAAc                                                                                                                                    <  1:576875/148‑1 (MQ=255)
                   aTACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAAcc                                                                                                                                   >  1:731996/1‑149 (MQ=255)
                     aCGCCACATCCGGCACAAGCAATAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTg                                                                                                                                 <  1:1005057/149‑1 (MQ=255)
                      cGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGc                                                                                                                                >  1:235322/1‑149 (MQ=255)
                       gCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCCCGCCTTGTTTCGTCAACGAAGAAAACGGTTCAACCTGcc                                                                                                                               >  2:58880/1‑148 (MQ=255)
                       gCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGc                                                                                                                                >  1:64206/1‑148 (MQ=255)
                       gCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGc                                                                                                                                <  2:64206/148‑1 (MQ=255)
                        ccACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGcac                                                                                                                              >  1:276774/1‑149 (MQ=255)
                         cacaTCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTACCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGcaca                                                                                                                             >  1:353814/1‑149 (MQ=255)
                            aTCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCa                                                                                                                          <  2:565196/149‑1 (MQ=255)
                            aTCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCa                                                                                                                          <  2:916344/149‑1 (MQ=255)
                               cGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTAGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACAGCACCg                                                                                                                       >  2:376306/1‑149 (MQ=255)
                                    caAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTaaa                                                                                                                  <  2:9386/149‑1 (MQ=255)
                                     aaGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGGTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAAc                                                                                                                 <  1:554970/149‑1 (MQ=255)
                                        cATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGcc                                                                                                              <  1:362722/149‑1 (MQ=255)
                                             aGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTAc                                                                                                         <  1:725882/149‑1 (MQ=255)
                                             aGGAAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAAAGTTTCAAACTGCAAATCACCGATAAAAGCCAGAAc                                                                                                         >  2:617326/1‑149 (MQ=255)
                                                cAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGc                                                                                                      <  1:719866/149‑1 (MQ=255)
                                                 aaGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGAACATCACCGTTAAACGCCAGTACAGCt                                                                                                     >  1:714846/1‑149 (MQ=255)
                                                     aaaTTATTCGACGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTGGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCAc                                                                                                 <  2:569480/149‑1 (MQ=255)
                                                         tattCGCCGTCCTGCGTTTCTTCTACAGTCTGCATCTCGCTAAAACAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCAc                                                                                             <  1:580416/149‑1 (MQ=255)
                                                         tattCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAACCCAGGGATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGCAGAAAACGTTTCAACCTGCACATCCCCGGTAAGCGCCAGTACAGCTTCACGccc                                                                                             >  2:195876/1‑147 (MQ=255)
                                                              gCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGATTCACGCACCatat                                                                                        >  1:812125/1‑149 (MQ=255)
                                                               ccGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGACAATCACAGTTAAACGCCAGTACAGCTTCACGCACCATAtt                                                                                       >  1:218198/1‑149 (MQ=255)
                                                               ccGTCCTGCGTTTCTTCTACAGGCTGCATCTCGATAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACTTTTCAACCTGCACATAACCGTTAATCTCCAGTACAGCTTCACGCACCATAtt                                                                                       >  1:315816/1‑149 (MQ=255)
                                                               cagTCCTGAGTTTATTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATAtt                                                                                       <  2:700908/147‑1 (MQ=255)
                                                                     tGCGCTTTTTCTACAGGGTGCATCGCGCTAATCAATGTATCCACTTGCAGCCGCAGCTTGTCAACGACTTGTTTCTTCAACTAAGAAAAAGTTTCACCCTGCACATCACCCTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTg                                                                                 <  2:714348/149‑1 (MQ=255)
                                                                      gCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTgc                                                                                >  2:384448/1‑149 (MQ=255)
                                                                      gCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTgc                                                                                <  2:945267/149‑1 (MQ=255)
                                                                        gTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTgcgg                                                                              >  2:261892/1‑148 (MQ=255)
                                                                        gTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGAATCCAGTATCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAACAAAACGTTTCACCCTGCACAACACCGTAAAACGCCAGTACACCTTCACGCACCATATTCAATTgcgc                                                                              >  2:423887/1‑149 (MQ=255)
                                                                         tttCTTCTTCAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCt                                                                             >  2:384474/1‑149 (MQ=255)
                                                                                 aCAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTgc                                                                     >  1:292460/1‑149 (MQ=255)
                                                                                 aCAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTgc                                                                     >  1:373651/1‑149 (MQ=255)
                                                                                  cAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAACCGTTTCAACCTGCACATCACCGGTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTgcg                                                                    >  1:24182/1‑149 (MQ=255)
                                                                                      cTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGcgccgc                                                                <  2:731996/149‑1 (MQ=255)
                                                                                        gCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCTCTTTACGTGCGCCGCtt                                                              <  2:235322/149‑1 (MQ=255)
                                                                                        gCATCTCGCTAAACCAGGTATCCAGTTTCTGCCACAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCtt                                                              <  1:639363/149‑1 (MQ=255)
                                                                                         cATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTAc                                                                         <  2:699732/137‑1 (MQ=255)
                                                                                         cATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTAc                                                                         >  1:699732/1‑137 (MQ=255)
                                                                                         cATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTg                                                             <  2:453288/149‑1 (MQ=255)
                                                                                          aTCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTTCAGCTTCACGCACCATATTCAATTGCTCTTTACGTGCGCCGCTTGc                                                            <  2:522256/149‑1 (MQ=255)
                                                                                             tcGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCATCGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAg                                                         >  2:396913/1‑149 (MQ=255)
                                                                                             tcGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAg                                                         <  1:462989/149‑1 (MQ=255)
                                                                                             tcGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTAAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTTCCAg                                                         >  1:393891/1‑149 (MQ=255)
                                                                                              cGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTGCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGt                                                        <  1:58880/149‑1 (MQ=255)
                                                                                               gCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACCGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGccgggt                                                       >  2:757773/1‑145 (MQ=255)
                                                                                                   aaCCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCAGATTCAATTGCGCTTTACGTGCGCCGGTTGCCAGTTTGTc                                                   >  1:86960/1‑149 (MQ=255)
                                                                                                          tATCCATTTTCTGCCGCATCGTGTCCACGCCTTGTTTCTTCATCGAAGAAAACGGTTCAACCTGCACATCCCCGTTAAACGGCAGTACAGCTTCACGCACCATATTCAATAGCGCTTTACGTGAGCCGCTTGCCAGTTTGTCCGCTTTg                                            <  1:217629/149‑1 (MQ=255)
                                                                                                          tATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTg                                            <  1:764762/149‑1 (MQ=255)
                                                                                                          tATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTg                                            <  2:126971/149‑1 (MQ=255)
                                                                                                           aTCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTgg                                           <  2:262544/149‑1 (MQ=255)
                                                                                                               aGTTTCTGACGCAGCATGTCCTCGCATTTTTTCTTCAACGAAGAAAACGTTTCATCCTGCACATCACCGTTAAACTCCAGTACAGCATCAAGCACCATATTAAATTGCTCTTTACCTCCGCCGCTTGCCAGTTTGTCCGCTTTGGTcag                                       <  2:375220/149‑1 (MQ=255)
                                                                                                                     tGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCAc                                 <  1:68761/149‑1 (MQ=255)
                                                                                                                      gCCGCAGCTTGTCCACGCCTTGTTACTTCAACGAAGAAAACGATTCAACCTGCACATCACCGTAAAACTCCAGTACAGCTTCACGCCCCATATTCAATTCCGCTTTACGTGCGCCTCTTGCCAGTTTGTCAGCTTAGGTCAGCaccacc                                >  1:1011126/1‑149 (MQ=255)
                                                                                                                             cTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTAAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGAACgg                         >  1:483181/1‑149 (MQ=255)
                                                                                                                                   cACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCTCCATATTCAATTGCGCTTTACGTGCGCCGCTGGCCAGTTTGTCCGCTTTGGTCCGCAGCCCCAGAACGGCGatat                   >  2:67039/1‑149 (MQ=255)
                                                                                                                                        cTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGAACGGCGATATTGCTg              <  1:595027/149‑1 (MQ=255)
                                                                                                                                         ttGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCTCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGAACGGCGATATTGCTGt             <  1:171166/149‑1 (MQ=255)
                                                                                                                                         cgttttCTTCGTTGAAGAAAACGTTTCACCCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCAGCAGAACGGCGCTATTGCTGt             >  1:942417/4‑149 (MQ=255)
                                                                                                                                          gttttCTTCGTTGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGAACGGCGATATTGCTGTc            >  1:952709/3‑149 (MQ=255)
                                                                                                                                          gttttCTTCGTTGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGAACGGCGATATTGCTGTc            >  1:156790/3‑149 (MQ=255)
                                                                                                                                               cttcAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGAACGGCGATATTGCTGTCTACCg       >  1:177561/1‑149 (MQ=255)
                                                                                                                                                ttcAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGAACCATAATCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGACAGCATAACCATAACTGCGATATTGCTGTCTCCCGc      >  2:176043/1‑149 (MQ=255)
                                                                                                                                                    aCGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGAACGGCGATATTGCTGTCTACCGCCCAc  >  2:60447/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TATTGAATTTACGGGCCGGATACGCCACATCCGGCACAAGCATTAAGGCAAGAAAATTATTCGCCGTCCTGCGTTTCTTCTACAGGCTGCATCTCGCTAAACCAGGTATCCAGTTTCTGCCGCAGCTTGTCCACGCCTTGTTTCTTCAACGAAGAAAACGTTTCAACCTGCACATCACCGTTAAACGCCAGTACAGCTTCACGCACCATATTCAATTGCGCTTTACGTGCGCCGCTTGCCAGTTTGTCCGCTTTGGTCAGCAGCACCAGAACGGCGATATTGCTGTCTACCGCCCAC  >  NC_000913/4050077‑4050373

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