Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,950,295:1 +C 5.1% coding (228/1773 nt) aspS ← aspartyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,950,2951.C 5.1% 341.8 / 13.1 78coding (228/1773 nt)aspSaspartyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base . (39/35);  new base C (2/2);  total (41/37)
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.88e-01

TCAGACGCGCTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGCTACCAAGATCACGAC  >  NC_000913/1950156‑1950443
                                                                                                                                            |                                                                                                                                                    
tCAGACGCGCTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTg                                                                                                                                             >  4:274029/1‑149 (MQ=255)
tCAGACGCGCTTCTTCGGTGTGGACGAGGGTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTg                                                                                                                                             >  2:1062675/1‑149 (MQ=255)
  aGACGCGCTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGcc                                                                                                                                           >  1:457458/1‑149 (MQ=255)
   gACGCGCTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGccc                                                                                                                                          >  4:876179/1‑149 (MQ=255)
      gcgcTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTg                                                                                                                                       >  4:954986/1‑149 (MQ=255)
      gcgcTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTg                                                                                                                                       >  1:1059782/1‑149 (MQ=255)
       cgcTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGa                                                                                                                                      <  1:256605/149‑1 (MQ=255)
       cgcTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGa                                                                                                                                      <  3:600370/149‑1 (MQ=255)
        gcTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGGGCACCGTACGTGCGCGTGAc                                                                                                                                     <  4:687948/149‑1 (MQ=255)
        gcTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGAc                                                                                                                                     >  1:582435/1‑149 (MQ=255)
        gcTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGAc                                                                                                                                     >  1:735226/1‑149 (MQ=255)
           tcttcGGTGTTGACGTGGTTAGAGTCAAGCGTCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCACCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTg                                                                                                                                  >  1:298241/1‑149 (MQ=255)
           tcttcGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTg                                                                                                                                  >  4:142290/1‑149 (MQ=255)
              tcGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAAt                                                                                                                               <  1:577486/149‑1 (MQ=255)
                ggTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGc                                                                                                                             <  1:471589/149‑1 (MQ=255)
                  tgtTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAg                                                                                                                           <  3:686371/149‑1 (MQ=255)
                    ttGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGaa                                                                                                                         <  3:73559/149‑1 (MQ=255)
                    ttGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGaa                                                                                                                         <  3:274029/149‑1 (MQ=255)
                    ttGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGaa                                                                                                                         <  2:150122/149‑1 (MQ=255)
                    ttGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGaa                                                                                                                         <  2:92998/149‑1 (MQ=255)
                     tGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGAGAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAAc                                                                                                                        <  3:903644/149‑1 (MQ=255)
                      gACGTGGTTAGAGTCAAGCTGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGGGCACCGTACGTGCGCGTGa                                                                                                                                      >  4:1002372/1‑134 (MQ=255)
                      gACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGGGCACCGTACGTGCGCGTGa                                                                                                                                      <  3:1002372/134‑1 (MQ=255)
                        cGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCa                                                                                                                     >  1:496641/1‑149 (MQ=255)
                         gTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCAt                                                                                                                    >  3:313472/1‑149 (MQ=255)
                         gTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCAt                                                                                                                    >  4:365851/1‑149 (MQ=255)
                           ggTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTa                                                                                                                  <  3:157519/149‑1 (MQ=255)
                           ggTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTa                                                                                                                  <  2:1059782/149‑1 (MQ=255)
                            gTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTAc                                                                                                                 >  3:802196/1‑149 (MQ=255)
                                gagTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAg                                                                                                             >  2:459672/1‑149 (MQ=255)
                                 agTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGt                                                                                                            >  3:238491/1‑149 (MQ=255)
                                       gCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCagag                                                                                                      <  3:446396/149‑1 (MQ=255)
                                             gAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGc                                                                                                >  2:517959/1‑149 (MQ=255)
                                               aCATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCtt                                                                                              <  2:457458/149‑1 (MQ=255)
                                                  tCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTaa                                                                                           <  4:355099/149‑1 (MQ=255)
                                                  tCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTaa                                                                                           <  3:661319/149‑1 (MQ=255)
                                                  tCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTaa                                                                                           >  2:227608/1‑149 (MQ=255)
                                                   cTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAAc                                                                                          <  2:411458/149‑1 (MQ=255)
                                                   cTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAAc                                                                                          <  3:257107/149‑1 (MQ=255)
                                                            tgatgaTAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCa                                                                                 <  2:496641/149‑1 (MQ=255)
                                                               tgatAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGa                                                                              >  2:552348/1‑149 (MQ=255)
                                                                 atAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATc                                                                            <  2:282245/149‑1 (MQ=255)
                                                                    gTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCgg                                                                         >  2:178356/1‑149 (MQ=255)
                                                                     tCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGa                                                                        <  3:365851/149‑1 (MQ=255)
                                                                       aGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATc                                                                      <  1:459672/149‑1 (MQ=255)
                                                                        gCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCg                                                                     >  1:583406/1‑149 (MQ=255)
                                                                        gCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCg                                                                     >  1:146919/1‑149 (MQ=255)
                                                                           aGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGaaa                                                                  >  3:370235/1‑149 (MQ=255)
                                                                             gACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGaaaaa                                                                <  2:146919/149‑1 (MQ=255)
                                                                                  cAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCt                                                           >  4:412732/1‑149 (MQ=255)
                                                                                  cAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCt                                                           >  2:1038535/1‑149 (MQ=255)
                                                                                  cAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCt                                                           <  1:517959/149‑1 (MQ=255)
                                                                                      aCTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCAc                                                       >  2:410480/1‑149 (MQ=255)
                                                                                        ttCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGa                                                     >  2:514278/1‑149 (MQ=255)
                                                                                        ttCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGa                                                     >  1:452886/1‑149 (MQ=255)
                                                                                        ttCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGa                                                     >  4:1028979/1‑149 (MQ=255)
                                                                                        ttCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGa                                                     >  2:713740/1‑149 (MQ=255)
                                                                                          cGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATa                                                   <  3:589736/149‑1 (MQ=255)
                                                                                                        gCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAATTCATTACGCAGTTCagag                                                                                                      >  4:452340/1‑84 (MQ=255)
                                                                                                        gCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCagag                                                                                                      <  3:452340/84‑1 (MQ=255)
                                                                                                        gCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTcgcg                                     >  1:862014/1‑149 (MQ=255)
                                                                                                              tCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATc                               >  1:689127/1‑149 (MQ=255)
                                                                                                                 cgGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGAt                            <  2:582435/149‑1 (MQ=255)
                                                                                                                    ttAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGaa                         <  4:313472/149‑1 (MQ=255)
                                                                                                                    ttAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGaa                         <  2:298241/149‑1 (MQ=255)
                                                                                                                    ttAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGaa                         <  3:142290/149‑1 (MQ=255)
                                                                                                                    ttAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGACCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGaa                         >  4:119942/1‑149 (MQ=255)
                                                                                                                      aaTATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGa                       >  3:226174/1‑149 (MQ=255)
                                                                                                                         atTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGa                             <  3:526996/140‑1 (MQ=255)
                                                                                                                         atTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGa                             >  4:526996/1‑140 (MQ=255)
                                                                                                                         atTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCa                    <  3:412732/149‑1 (MQ=255)
                                                                                                                               ttCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGCTAc              <  3:876179/149‑1 (MQ=255)
                                                                                                                                  gtgaCGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGCTCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGCTACCag           >  1:312530/4‑148 (MQ=255)
                                                                                                                                  gTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGCTACaaa           >  1:668375/1‑149 (MQ=255)
                                                                                                                                   tCACGCGCACGTACGGTGCCCGTGACCTGAATGTAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGCTACCAAg          >  4:487432/1‑149 (MQ=255)
                                                                                                                                    cACGGGCACCGTACGTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGCTACCAAGa         >  3:1043384/1‑149 (MQ=255)
                                                                                                                                       gcgcACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGCTACCAAGATCa      <  4:822453/149‑1 (MQ=255)
                                                                                                                                       gcgcACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGCTACCAAGATCa      <  2:328129/149‑1 (MQ=255)
                                                                                                                                           aCGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGCTACCAAGATCacgac  <  4:378903/149‑1 (MQ=255)
                                                                                                                                           aCGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGCTACCAAGATCacgac  >  4:994846/1‑149 (MQ=255)
                                                                                                                                            |                                                                                                                                                    
TCAGACGCGCTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCACGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGCTACCAAGATCACGAC  >  NC_000913/1950156‑1950443

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