Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,180,513:1 +A 5.8% coding (93/1161 nt) ygiC → ATP‑Grasp family ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,180,5131.A 5.8% 523.2 / 27.0 121coding (93/1161 nt)ygiCATP‑Grasp family ATPase
Reads supporting (aligned to +/- strand):  ref base . (52/62);  new base A (2/5);  total (54/67)
Fisher's exact test for biased strand distribution p-value = 4.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

Only 100 of 125 total aligned reads displayed.

ACTATGCAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAA  >  NC_000913/3180365‑3180660
                                                                                                                                                     |                                                                                                                                                   
aCTATGCAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTAc                                                                                                                                                      <  4:577422/149‑1 (MQ=255)
  tATGCAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTg                                                                                                                                                   <  1:777750/149‑1 (MQ=255)
    tGCAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGgt                                                                                                                                                 <  1:1105259/149‑1 (MQ=255)
    tGCAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGgt                                                                                                                                                 >  4:9274/1‑149 (MQ=255)
    tGCAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGgt                                                                                                                                                 <  1:652219/149‑1 (MQ=255)
     gCAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCAGTACGGCTCGCCGTACATGgtg                                                                                                                                                <  1:14602/149‑1 (MQ=255)
      cAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCAGTACGGCTCGCCGTACATGgtgt                                                                                                                                               <  3:511009/149‑1 (MQ=255)
      cAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCAGTACGGCTCGCCGTACATGgtgt                                                                                                                                               <  1:1137693/149‑1 (MQ=255)
       aGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGtgtg                                                                                                                                              >  4:208533/1‑149 (MQ=255)
          gTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGgt                                                                                                                                                 <  3:613141/143‑1 (MQ=255)
          gTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGgt                                                                                                                                                 >  4:613141/1‑143 (MQ=255)
             gTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCAGTACGGCTCGCCGTACATGGTGTGAAAAtt                                                                                                                                        <  1:47695/149‑2 (MQ=38)
              tGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGc                                                                                                                                       <  3:444234/149‑1 (MQ=255)
              tGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGc                                                                                                                                       >  2:327786/1‑149 (MQ=255)
              tGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGc                                                                                                                                       <  3:186860/149‑1 (MQ=255)
              tGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGc                                                                                                                                       <  3:9274/149‑1 (MQ=255)
                 aaCCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTta                                                                                                                                    >  2:488422/1‑149 (MQ=255)
                       tACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACaa                                                                                                                              <  4:718356/149‑1 (MQ=255)
                       tACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCAGTACGGCTCGCCGTAcatggt                                                                                                                                                 >  3:154008/1‑125 (MQ=255)
                       tACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCAGTACGGCTCGCCGTAcatggt                                                                                                                                                 <  4:154008/130‑6 (MQ=255)
                         ccTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGt                                                                                                                            >  2:141290/1‑149 (MQ=255)
                          ctctTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGtt                                                                                                                           <  4:60884/149‑1 (MQ=255)
                           tcttctCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTg                                                                                                                          >  2:535956/1‑149 (MQ=255)
                            cttctCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCTATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGa                                                                                                                         <  2:527745/149‑1 (MQ=255)
                               ctcGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGAccc                                                                                                                      >  2:525638/1‑149 (MQ=255)
                                tcGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGtt                                                                                                                           >  3:968757/1‑143 (MQ=255)
                                     cAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGccc                                                                                                                <  1:652661/149‑1 (MQ=255)
                                        tGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAgg                                                                                                             >  3:373768/1‑149 (MQ=255)
                                           gTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTg                                                                                                          >  4:549185/1‑149 (MQ=255)
                                             ggCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGaa                                                                                                        <  1:525638/149‑1 (MQ=255)
                                                 gATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGc                                                                                                    <  1:565716/149‑1 (MQ=255)
                                                  aTACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCt                                                                                                   <  4:135457/149‑1 (MQ=255)
                                                  aTACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCt                                                                                                   <  2:795848/149‑1 (MQ=255)
                                                  aTACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCt                                                                                                   <  3:264539/149‑1 (MQ=255)
                                                   tACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTg                                                                                                  <  4:569586/149‑1 (MQ=255)
                                                      cGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGgaa                                                                                               <  3:649854/149‑1 (MQ=255)
                                                        aTGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAg                                                                                                              <  3:90014/132‑1 (MQ=255)
                                                        aTGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAg                                                                                                              >  4:90014/1‑132 (MQ=255)
                                                         tGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCTCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGCTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGgaagaa                                                                                            <  3:801147/149‑1 (MQ=255)
                                                         tGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGgaagaa                                                                                            <  3:208533/149‑1 (MQ=255)
                                                           gAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGt                                                                                          >  1:240849/1‑149 (MQ=255)
                                                              agagTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCt                                                                                                   >  2:782866/1‑137 (MQ=255)
                                                                 gTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCAccgc                                                                                    >  4:898470/1‑149 (MQ=255)
                                                                  tCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCAccgcc                                                                                   >  1:559544/1‑149 (MQ=255)
                                                                    aGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCt                                                                                                   <  3:578138/131‑1 (MQ=255)
                                                                      tATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTta                                                                                                                                    <  3:735978/96‑1 (MQ=255)
                                                                      tATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTta                                                                                                                                    >  4:735978/1‑96 (MQ=255)
                                                                       aTTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACt                                                                              >  4:427057/1‑149 (MQ=255)
                                                                        ttACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTg                                                                             >  2:909451/1‑149 (MQ=255)
                                                                        ttACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTg                                                                             >  4:1126427/1‑149 (MQ=255)
                                                                           ccGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCAc                                                                          <  1:535956/149‑1 (MQ=255)
                                                                           ccGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCAc                                                                          <  1:842866/149‑1 (MQ=255)
                                                                                   ccGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACt                                                                              >  1:157319/1‑137 (MQ=255)
                                                                                     ggACTGGCGTGACAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGAtgtgtc                                                                <  1:909451/149‑3 (MQ=255)
                                                                                       aCTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTg                                                              <  3:784805/149‑1 (MQ=255)
                                                                                           gCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAg                                                          >  4:904348/1‑149 (MQ=255)
                                                                                              tGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAgtgg                                                       <  4:373768/149‑1 (MQ=255)
                                                                                               gagaAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAgtggt                                                      >  2:395758/1‑149 (MQ=255)
                                                                                                agaAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAgtggtg                                                     <  3:983733/149‑1 (MQ=255)
                                                                                                  aaaGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGa                                                   >  1:907238/1‑149 (MQ=255)
                                                                                                  aaaGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGa                                                   >  1:52904/1‑149 (MQ=255)
                                                                                                      cccACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGaaaaa                                               <  3:549185/149‑1 (MQ=255)
                                                                                                      cccACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGaaaaa                                               <  1:24896/149‑1 (MQ=255)
                                                                                                       ccACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGc                                                                                                                                       <  3:429223/60‑1 (MQ=255)
                                                                                                       ccACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGc                                                                                                                                       >  4:429223/1‑60 (MQ=255)
                                                                                                       ccACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGAccc                                                                                                                      <  1:194841/77‑1 (MQ=255)
                                                                                                       ccACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGAccc                                                                                                                      >  2:194841/1‑77 (MQ=255)
                                                                                                       ccACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAg                                              >  4:1050847/1‑149 (MQ=255)
                                                                                                         aCGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTg                                            >  2:520388/1‑149 (MQ=255)
                                                                                                          cGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGa                                           <  2:240849/149‑1 (MQ=255)
                                                                                                           gAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGAt                                          <  1:141290/149‑1 (MQ=255)
                                                                                                            aaTACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATc                                         <  1:520388/149‑1 (MQ=255)
                                                                                                                  gTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAGGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGAt                                          <  1:247025/142‑1 (MQ=255)
                                                                                                                  gTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAGGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGAt                                          >  2:247025/1‑142 (MQ=255)
                                                                                                                  gTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGc                                   <  2:907238/149‑1 (MQ=255)
                                                                                                                     tCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGAt                                <  3:427057/149‑1 (MQ=255)
                                                                                                                        aTTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAgtg                                                        <  3:531593/122‑1 (MQ=255)
                                                                                                                        aTTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAgtg                                                        >  4:531593/1‑122 (MQ=255)
                                                                                                                         tcttCACACCAGTACGGCTCGCCGTACATGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGc                            >  1:1052358/3‑149 (MQ=38)
                                                                                                                          tttCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCt                           <  2:559544/149‑1 (MQ=255)
                                                                                                                              acacCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATg                               >  2:466284/1‑141 (MQ=255)
                                                                                                                              acacCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATg                               <  1:466284/141‑1 (MQ=255)
                                                                                                                               cacCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCt                           <  2:638410/144‑1 (MQ=255)
                                                                                                                               cacCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCt                           >  1:638410/1‑144 (MQ=255)
                                                                                                                                acCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGa                                 >  2:605811/1‑137 (MQ=255)
                                                                                                                                acCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCtgat                        <  1:200147/146‑1 (MQ=255)
                                                                                                                                acCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCtgat                        >  2:200147/1‑146 (MQ=255)
                                                                                                                                acCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGAc                     >  1:378356/1‑149 (MQ=255)
                                                                                                                                 ccATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGAt                                                                     >  4:816698/1‑100 (MQ=255)
                                                                                                                                  cATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGAtg                                                                    <  4:418443/100‑1 (MQ=255)
                                                                                                                                  cATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGAtg                                                                    >  3:418443/1‑100 (MQ=255)
                                                                                                                                  cATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATg                               <  3:839466/137‑1 (MQ=255)
                                                                                                                                  cATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATg                               >  4:839466/1‑137 (MQ=255)
                                                                                                                                   cagTACGGCTCGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGAt                                                                     <  4:610013/96‑1 (MQ=255)
                                                                                                                                   cagTACGGCTCGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGAt                                                                     >  3:610013/3‑98 (MQ=255)
                                                                                                                                        cGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTcc             <  3:898470/149‑1 (MQ=255)
                                                                                                                                         ggCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCg            >  3:759237/1‑149 (MQ=255)
                                                                                                                                          gCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGc           >  4:1042252/1‑149 (MQ=255)
                                                                                                                                               ccGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTcc      <  2:52904/149‑1 (MQ=255)
                                                                                                                                                   aCTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCaaaa  <  4:796766/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
ACTATGCAGCGTAGTGCAACCGGTACCTCTTCTCGTTCAATGGGTGGCTGATACCGATGGAAAGAGTCAGTATTACCGAGCGCCCGGACTGGCGTGAGAAAGCCCACGAATACGGTTTCAATTTTCACACCATGTACGGCGAGCCGTACTGGTGTGAAGATGCTTACTACAAGTTGACCCTCGCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAA  >  NC_000913/3180365‑3180660

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: