Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,180,505 G→T 12.2% E29* (GAG→TAG) ‡ ygiC → ATP‑Grasp family ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,180,5050GT12.2% 598.0 / 54.8 179E29* (GAG→TAG) ‡ygiCATP‑Grasp family ATPase
Reads supporting (aligned to +/- strand):  ref base G (80/77);  new base T (11/11);  total (91/88)
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.16e-01

Only 100 of 183 total aligned reads displayed.

AACAAAGCACTATGCAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCAT  >  NC_000913/3180357‑3180653
                                                                                                                                                    |                                                                                                                                                     
   aaaGCACTATGCAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGc                                                                                                                                                    >  1:1171004/1‑149 (MQ=255)
           aTGCAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTgg                                                                                                                                           <  2:348626/149‑1 (MQ=255)
              cAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGgtgt                                                                                                                                        <  2:24680/149‑1 (MQ=255)
                 cGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAACACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGaa                                                                                                                                     <  3:1290574/149‑1 (MQ=255)
                        cAACCGGTCCCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCtt                                                                                                                              >  4:1680907/1‑149 (MQ=255)
                        cAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCtt                                                                                                                              >  3:1073571/1‑149 (MQ=255)
                          aCCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTtac                                                                                                                            <  2:1171004/149‑1 (MQ=255)
                           ccGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTtact                                                                                                                           >  1:622696/1‑149 (MQ=255)
                              gTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGc                                                                                                                                >  1:16618/1‑141 (MQ=255)
                              gTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGc                                                                                                                                >  1:177090/1‑141 (MQ=255)
                                aCCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGgtgt                                                                                                                                        >  2:1612989/1‑131 (MQ=255)
                                  ctctTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGaa                                                                                                                                     <  2:1238277/132‑1 (MQ=255)
                                   tcttctCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCt                                                                                                                               >  3:698058/1‑137 (MQ=255)
                                   tcttctCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCt                                                                                                                               <  4:698058/137‑1 (MQ=255)
                                   tcttctCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTtac                                                                                                                            <  3:980106/140‑1 (MQ=255)
                                   tcttctCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTtac                                                                                                                            >  4:980106/1‑140 (MQ=255)
                                   tcttctCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCAGTACGGCTCGCCGTAcatggt                                                                                                                                          <  2:644476/126‑6 (MQ=255)
                                     ttctCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCAGTACGGCTCGCCGTAcatgg                                                                                                                                           <  1:16345/123‑5 (MQ=255)
                                        tcGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGa                                                                                                                                      <  1:550183/125‑1 (MQ=255)
                                        tcGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGa                                                                                                                                      >  2:550183/1‑125 (MQ=255)
                                        tcGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCAGTACGGCTCGCCGTAcatggt                                                                                                                                          >  4:1777624/1‑116 (MQ=39)
                                        tcGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCAGTACGGCTCGCCGTAcatggt                                                                                                                                          <  3:1777624/121‑6 (MQ=39)
                                          gTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCAGTACGGCTCGCCGTACATGgtg                                                                                                                                         <  4:1257401/120‑1 (MQ=38)
                                                  ggTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGtt                                                                                                    >  2:696330/1‑149 (MQ=255)
                                                   gTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTg                                                                                                   <  1:372942/149‑1 (MQ=255)
                                                   gTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTg                                                                                                   <  4:1073571/149‑1 (MQ=255)
                                                        tGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGcc                                                                                                                                                   <  4:1523557/97‑1 (MQ=255)
                                                        tGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGcc                                                                                                                                                   >  3:1523557/1‑97 (MQ=255)
                                                        tGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTta                                                                                                                             <  1:1513356/118‑1 (MQ=255)
                                                        tGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAg                                                                                              >  4:1232975/1‑149 (MQ=255)
                                                        tGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCAGTACGGCTCGCCGTAcatggt                                                                                                                                          <  3:1254912/105‑6 (MQ=38)
                                                                  ggAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCt                                                                                            >  1:1226456/1‑141 (MQ=255)
                                                                    aaaGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTc                                                                                  >  1:917887/1‑149 (MQ=255)
                                                                          tCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCt                                                                                            >  3:438393/1‑133 (MQ=255)
                                                                          tCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCt                                                                                            <  4:438393/133‑1 (MQ=255)
                                                                          tCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCAccgcc                                                                            >  4:560699/1‑149 (MQ=255)
                                                                          tCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCAccgcc                                                                            >  2:726065/1‑149 (MQ=255)
                                                                          tCAGTACTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAg                                                                                                                      <  4:628040/107‑1 (MQ=255)
                                                                           cAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCAccgccg                                                                           >  2:1507036/1‑149 (MQ=255)
                                                                            aGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGa                                                                          >  2:384791/1‑149 (MQ=255)
                                                                             gTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGaa                                                                         <  2:622696/149‑1 (MQ=255)
                                                                                 tACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTg                                                                                           >  2:1763430/1‑127 (MQ=255)
                                                                                   ccGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGcc                                                                                                                                                   <  1:1281227/70‑1 (MQ=255)
                                                                                   ccGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCa                                                                    >  2:1703571/1‑148 (MQ=255)
                                                                                       gcgcCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGa                                                               >  1:290299/1‑149 (MQ=255)
                                                                                        cgcCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGgaagaa                                                                                     >  3:1149867/1‑126 (MQ=255)
                                                                                           ccGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAc                                                                                                                                               <  1:684467/66‑1 (MQ=255)
                                                                                            cGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGa                                                                                                                  <  1:1752659/93‑1 (MQ=255)
                                                                                            cGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGc                                                          <  1:384791/149‑1 (MQ=255)
                                                                                             ggACTGGCATGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGgaagaa                                                                                     <  4:462998/121‑1 (MQ=255)
                                                                                             ggACTGGCATGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGgaagaa                                                                                     >  3:462998/1‑121 (MQ=255)
                                                                                                   gCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCt                                                                                                                               <  3:645744/73‑1 (MQ=255)
                                                                                                   gCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCt                                                                                                                               >  4:645744/1‑73 (MQ=255)
                                                                                                   gCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCAGTACGGCTCGCCGTACATGGTGTGAAAAttg                                                                                                                                >  3:622568/1‑70 (MQ=14)
                                                                                                      tGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAgtgg                                                <  3:560699/149‑1 (MQ=255)
                                                                                                       gagaAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAgtggt                                               <  1:640270/149‑1 (MQ=255)
                                                                                                        agaAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAgtggtg                                              <  2:290299/149‑1 (MQ=255)
                                                                                                        agaAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAgtggtg                                              >  2:1648791/1‑149 (MQ=255)
                                                                                                        agaAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAgtggtg                                              <  1:726065/149‑1 (MQ=255)
                                                                                                         gaAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTTAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAGGGtggtgg                                             >  1:324773/1‑149 (MQ=255)
                                                                                                              cccACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTg                                                                      <  2:1070705/119‑1 (MQ=255)
                                                                                                              cccACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTg                                                                      >  1:1070705/1‑119 (MQ=255)
                                                                                                               ccACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAg                                       >  2:210830/1‑149 (MQ=255)
                                                                                                                cACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGt                                      >  2:479266/1‑149 (MQ=255)
                                                                                                                 aCGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGaaa                                                                                                <  2:117973/90‑1 (MQ=255)
                                                                                                                   gAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGAt                                   <  1:1507036/149‑1 (MQ=255)
                                                                                                                    aaTACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATc                                  <  2:324773/149‑1 (MQ=255)
                                                                                                                     aTACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCg                                 >  2:885664/1‑149 (MQ=255)
                                                                                                                        cGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCa                              >  4:268886/1‑149 (MQ=255)
                                                                                                                           tttCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTc                                                                                                             <  3:1399626/67‑1 (MQ=255)
                                                                                                                           tttCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTc                                                                                                             >  4:1399626/1‑67 (MQ=255)
                                                                                                                           tttCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATc                                  <  4:1513391/142‑1 (MQ=255)
                                                                                                                           tttCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCg                           >  4:867034/1‑149 (MQ=255)
                                                                                                                            ttCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGa                          >  1:1018417/1‑149 (MQ=255)
                                                                                                                             tCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGAt                         <  3:867034/149‑1 (MQ=255)
                                                                                                                               aaTTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAgtggt                                               >  4:1138523/1‑125 (MQ=255)
                                                                                                                               aaTTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAgtggt                                               <  3:1138523/125‑1 (MQ=255)
                                                                                                                               aaTTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGa                       <  2:1372441/149‑1 (MQ=255)
                                                                                                                                 tcttCACACCAGTACGGCTCGCCGTACATGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAgtggt                                               >  3:1517870/3‑123 (MQ=38)
                                                                                                                                   ttCACACCAGTACGGCTCGCCGTACATGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCt                                                        <  3:8168/112‑1 (MQ=38)
                                                                                                                                   ttCACACCAGTACGGCTCGCCGTACATGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCt                                                        >  4:8168/1‑112 (MQ=38)
                                                                                                                                    tCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCtga                  <  1:210830/149‑1 (MQ=255)
                                                                                                                                      acacCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGa                                                               <  1:180099/102‑1 (MQ=255)
                                                                                                                                      acacCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGa                                                               >  2:180099/1‑102 (MQ=255)
                                                                                                                                      acacCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCtgatg                <  3:268886/149‑1 (MQ=255)
                                                                                                                                      acacCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCtgatg                >  3:1598010/1‑149 (MQ=255)
                                                                                                                                      acacCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCtgatg                >  2:1634133/1‑149 (MQ=255)
                                                                                                                                        acCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGcc                                                         <  4:1231048/106‑1 (MQ=255)
                                                                                                                                        acCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGcc                                                         >  3:1231048/1‑106 (MQ=255)
                                                                                                                                        acCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGAt                         >  1:706863/1‑138 (MQ=255)
                                                                                                                                         ccATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGAcc             <  4:116308/149‑1 (MQ=255)
                                                                                                                                          cATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCAc                                                                   <  2:729440/94‑1 (MQ=255)
                                                                                                                                          cATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAg                                       <  1:1586090/122‑1 (MQ=255)
                                                                                                                                          cATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCa            >  2:1562977/1‑149 (MQ=255)
                                                                                                                                            tGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAgtg                                                 <  2:776905/110‑1 (MQ=255)
                                                                                                                                            tGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAgtg                                                 >  1:776905/1‑110 (MQ=255)
                                                                                                                                            tGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAgtg                                                 >  1:444229/1‑110 (MQ=255)
                                                                                                                                             gTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAAt         >  3:887545/1‑149 (MQ=255)
                                                                                                                                                    gAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGc                                                          >  1:810208/1‑93 (MQ=255)
                                                                                                                                                    gAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGc                                                          <  2:810208/93‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AACAAAGCACTATGCAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCAT  >  NC_000913/3180357‑3180653

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