Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 732,660 G→A intergenic (+114/‑426) tyrS → / → PA14_RS03450 tyrosine‑‑tRNA ligase/16S ribosomal RNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_008463732,6600GA93.3% 294.5 / 8.4 90intergenic (+114/‑426)tyrS/PA14_RS03450tyrosine‑‑tRNA ligase/16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base G (2/4);  new base A (55/28);  total (58/32)
Fisher's exact test for biased strand distribution p-value = 1.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCTTTCGCGCGCATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGGGTGGTGAAGCGCTCCGGCAGGTTGTTCGCCGCAGCGGTTCGGTCCCGAAAGGGGCTCGAAACGAAGCTTCATCGAGGTGCTTGACAGCGAATTTGAACGCTGTAGAATGCGCCTCCCGCTGATCGGAAGATAGTTTGAAGGTCG  >  NC_008463/732513‑732804
                                                                                                                                                   |                                                                                                                                                
ggCTTTCGCGCGCATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGagt                                                                                                                                                <  2:795993/150‑1 (MQ=255)
 gCTTTCGCGCGCATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGagtg                                                                                                                                               <  1:9134173/150‑1 (MQ=255)
 gCTTTCGCGCGCATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGagtg                                                                                                                                               >  2:5358700/1‑150 (MQ=255)
  cTTTCGCGCGCATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTNACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGagtgg                                                                                                                                              >  1:5798579/1‑150 (MQ=255)
  cTTTCGCGCGCATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGagtgg                                                                                                                                              >  1:4321777/1‑150 (MQ=255)
  cTTTCGCGCGCATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGagtgg                                                                                                                                              >  1:1328285/1‑150 (MQ=255)
    ttCGCGCGCATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCTCCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAgtggtg                                                                                                                                            <  2:7744926/150‑1 (MQ=255)
    ttCGCGCGCATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAgtggtg                                                                                                                                            >  1:9608345/1‑150 (MQ=255)
     tCGCGCGCATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGa                                                                                                                                           <  1:5118103/150‑1 (MQ=255)
       gcgcgcATCACGCTGAAGGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAg                                                                                                                                         <  2:3033800/150‑1 (MQ=255)
        cgcgcATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAgc                                                                                                                                        <  1:5428930/150‑1 (MQ=255)
        cgcgcATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAgc                                                                                                                                        >  2:2762612/1‑150 (MQ=255)
        cgcgcATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAgc                                                                                                                                        >  2:8640574/1‑150 (MQ=255)
        cgcgcATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAgc                                                                                                                                        >  1:7836868/1‑150 (MQ=255)
         gcgcATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAgcg                                                                                                                                       >  1:5690976/1‑150 (MQ=255)
         gcgcATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAgcg                                                                                                                                       >  2:6568290/1‑150 (MQ=255)
         gcgcATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAgcg                                                                                                                                       >  1:32183/1‑150 (MQ=255)
          cgcATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAgcgc                                                                                                                                      >  1:8166805/1‑150 (MQ=255)
           gcATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCt                                                                                                                                     >  2:9500709/1‑150 (MQ=255)
           gcATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCt                                                                                                                                     >  1:3405713/1‑150 (MQ=255)
            cATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTc                                                                                                                                    <  1:401881/150‑1 (MQ=255)
            cATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTc                                                                                                                                    >  2:6139485/1‑150 (MQ=255)
            cATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTc                                                                                                                                    <  2:4154266/150‑1 (MQ=255)
             aTCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTcc                                                                                                                                   >  1:795695/1‑150 (MQ=255)
             aTCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTcc                                                                                                                                   >  1:2635933/1‑150 (MQ=255)
             aTCACGCTGAAAGCTGAATAAGTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAAGGCGCACCACTCCCGGCGACGACGAAGTGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTcc                                                                                                                                   <  1:7197864/150‑1 (MQ=255)
              tCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAATCGCTCCg                                                                                                                                  >  1:6407372/1‑150 (MQ=255)
              tCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCg                                                                                                                                  >  2:8615527/1‑150 (MQ=255)
               cACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCgg                                                                                                                                 <  1:5070527/150‑1 (MQ=255)
               cACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCgg                                                                                                                                 <  1:5060260/150‑1 (MQ=255)
                aCGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGc                                                                                                                                <  2:8352416/150‑1 (MQ=255)
                aCGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGc                                                                                                                                >  2:1026289/1‑150 (MQ=255)
                 cGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTTCGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCa                                                                                                                               <  1:1205162/150‑1 (MQ=255)
                 cGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGct                                                                                                                               >  1:7779290/1‑149 (MQ=255)
                 cGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCa                                                                                                                               >  1:3336867/1‑150 (MQ=255)
                 cGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCa                                                                                                                               >  1:3327088/1‑150 (MQ=255)
                   cTGACAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAgg                                                                                                                             >  2:9754382/1‑150 (MQ=255)
                    tGACAGCTGACTAATTCGTTCATATTAGTGTTGACGCCAGTTTGTTTCTCCCTCTAAGGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCACGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGgt                                                                                                                            <  2:6421182/150‑1 (MQ=11)
                    tGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGgt                                                                                                                            >  1:6772707/1‑150 (MQ=255)
                    tGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGgt                                                                                                                            >  2:4597793/1‑150 (MQ=255)
                     gAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGgtt                                                                                                                           <  1:9809279/150‑1 (MQ=255)
                     gAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGgtt                                                                                                                           >  2:1147164/1‑150 (MQ=255)
                     gAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGgtt                                                                                                                           >  1:7889317/1‑150 (MQ=255)
                      aaaGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGgttg                                                                                                                          <  1:10829752/150‑1 (MQ=255)
                       aaGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCTGgttgt                                                                                                                         >  2:2689830/1‑150 (MQ=255)
                       aaGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGgttgt                                                                                                                         <  2:8404382/150‑1 (MQ=38)
                       aaGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGgttgt                                                                                                                         >  1:5704801/1‑150 (MQ=38)
                       aaGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGgttgt                                                                                                                         >  1:5227442/1‑150 (MQ=38)
                       aaGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGgttgt                                                                                                                         >  2:7713527/1‑150 (MQ=38)
                       aaGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGgttgt                                                                                                                         >  2:2689865/1‑150 (MQ=38)
                        aGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATACTGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACAGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGgttgtt                                                                                                                        <  1:4134386/150‑1 (MQ=255)
                        aGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGgttgtt                                                                                                                        >  1:696778/1‑150 (MQ=38)
                        aGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGgttgtt                                                                                                                        >  1:695550/1‑150 (MQ=38)
                        aGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGgttgtt                                                                                                                        <  2:1466963/150‑1 (MQ=38)
                        aGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTAGAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGgttgtt                                                                                                                        >  2:6510009/1‑150 (MQ=255)
                         gCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTc                                                                                                                       >  2:2435516/1‑150 (MQ=37)
                         gCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTc                                                                                                                       >  2:8410157/1‑150 (MQ=37)
                         gCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTc                                                                                                                       >  2:8412605/1‑150 (MQ=37)
                         gCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTc                                                                                                                       >  2:3841788/1‑150 (MQ=37)
                         gCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTc                                                                                                                       <  2:9076557/150‑1 (MQ=37)
                         gCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTc                                                                                                                       >  2:7918050/1‑150 (MQ=37)
                         gCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTc                                                                                                                       <  1:7084365/150‑1 (MQ=37)
                          cTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGATGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGCGTGGTGAAGCGCTCCGGCAGGGTGTTcg                                                                                                                      >  2:8177446/1‑150 (MQ=255)
                          cTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTcg                                                                                                                      >  2:2675201/1‑150 (MQ=25)
                          cTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTcg                                                                                                                      >  1:4736795/1‑150 (MQ=25)
                           tGAATAATTCTTTCAAATCAGTGTTGACGCCAGTTTCTATCTCCCGATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTcgc                                                                                                                     <  2:186201/150‑1 (MQ=255)
                           tGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTcgc                                                                                                                     >  2:8009443/1‑150 (MQ=25)
                           tGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTcgc                                                                                                                     <  1:3993458/150‑1 (MQ=25)
                           tGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTcgc                                                                                                                     <  2:7797107/150‑1 (MQ=25)
                           tGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTcgc                                                                                                                     <  1:10869239/150‑1 (MQ=25)
                           tGAATAATTCTTTCAAATAAGTGTAGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTcgc                                                                                                                     <  1:10864573/150‑1 (MQ=25)
                            gAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGACGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTcgcc                                                                                                                    >  2:8206384/1‑150 (MQ=21)
                            gAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTcgcc                                                                                                                    >  1:2632811/1‑150 (MQ=21)
                            gAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTcgcc                                                                                                                    >  2:8744846/1‑150 (MQ=21)
                            gAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTcgcc                                                                                                                    >  2:8772114/1‑150 (MQ=21)
                            gAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTcgcc                                                                                                                    >  2:7616458/1‑150 (MQ=21)
                            gAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTcgcc                                                                                                                    >  2:10066024/1‑150 (MQ=21)
                             aataatTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTcgccg                                                                                                                   >  2:1300658/1‑150 (MQ=18)
                              ataatTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTcgccgc                                                                                                                  >  1:3564889/1‑150 (MQ=17)
                              ataatTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTcgccgc                                                                                                                  <  2:8793579/150‑1 (MQ=17)
                              ataatTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTcgccgc                                                                                                                  >  1:9800571/1‑150 (MQ=17)
                              ataatTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTcgccgc                                                                                                                  <  2:6803628/150‑1 (MQ=17)
                              ataatTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTcgccgc                                                                                                                  >  1:4949692/1‑150 (MQ=17)
                               taatTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTCGCCGCa                                                                                                                 >  2:2636080/1‑150 (MQ=17)
                                aatTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTCGCCGCAg                                                                                                                <  1:4056604/150‑1 (MQ=14)
                                aatTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTCGCCGCAg                                                                                                                >  1:5060772/1‑150 (MQ=14)
                                aatTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTCGCCGCAg                                                                                                                >  1:5692560/1‑150 (MQ=14)
                                aatTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTCGCCGCAg                                                                                                                >  1:5692619/1‑150 (MQ=14)
                                aatTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTCGCCGCAg                                                                                                                <  1:9869064/150‑1 (MQ=14)
                                aatTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGAGTGGTGAAGCGCTCCGGCAGGTTGTTCGCCGCAg                                                                                                                >  1:10827272/1‑150 (MQ=14)
                                           atcaGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGGGTGGTGAAGCGCTCCGGCAGGTTGTTCGCCGCAGCGGTTCGGTcc                                                                                                     <  2:7218955/147‑1 (MQ=1)
                                                           gtgtCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGGGTGGTGAAGCGCTCCGGCAGGTTGTTCGCCGCAGCGGTTCGGTCCCGAAAGGGGCTCGaaa                                                                                     <  1:3781381/147‑1 (MQ=1)
                                                                    tcacTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGGGTGGTGAAGCGCTCCGGCAGGTTGTTCGCCGCAGCGGTTCGGTCCCGAAAGGGGCTCGAAACGAAGCTTc                                                                            >  2:2087226/4‑150 (MQ=1)
                                                                    tcacTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGGGTGGTGAAGCGCTCCGGCAGGTTGTTCGCCGCAGCGGTTCGGTCCCGAAAGGGGCTCGAAACGAAGCTTc                                                                            >  2:2087216/4‑150 (MQ=1)
                                                                      ccgataATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGGGTGGTGAAGCGCTCCGGCAGGTTGTTCGCCGCAGCGGTTCGGTCCCGAAAGGGGCTCGAAACGAAGCTTCAt                                                                          <  2:6063863/147‑1 (MQ=1)
                                                                                                           aactGCTTGAAAATCAAGTATTTACGGATTCTCTCCCAGGGGTGGTGAAGCGCTCCGGCAGGTTGTTCGCCGCAGCGGTTCGGTCCCGAAAGGGGCTCGAAACGAAGCTTCATCGAGGTGCTTGACAGCGATTTTGAACGCTGTAGAATg                                     <  2:9575133/147‑1 (MQ=0)
                                                                                                                                              cccggggCGGTGAAGCGCTCCGGCAGGTAGTTCGCCGCCGCGGTTCGGTCCCGCAAGGGGCTCGAAACGCAGCTTCATCGAGGTGCTTGCCAGCGAATTTGAACGCTGTAGAATGCGCCTCCCGCTGACCGGACGATAGTTTGAAGGTCg  <  2:4702310/147‑1 (MQ=1)
                                                                                                                                                   |                                                                                                                                                
GGCTTTCGCGCGCATCACGCTGAAAGCTGAATAATTCTTTCAAATAAGTGTTGACGCCAGTTTCTATCTCCCTATAATGCGCACCACTCCCGGCGACGACGAAGAGTAAATGCTTGAAAATCAAGGATTTACGGATTCTCTCCCAGGGGTGGTGAAGCGCTCCGGCAGGTTGTTCGCCGCAGCGGTTCGGTCCCGAAAGGGGCTCGAAACGAAGCTTCATCGAGGTGCTTGACAGCGAATTTGAACGCTGTAGAATGCGCCTCCCGCTGATCGGAAGATAGTTTGAAGGTCG  >  NC_008463/732513‑732804

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