Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,313,586 T→C 5.1% F79F (TTT→TTC rpiB → allose‑6‑phosphate isomerase/ribose‑5‑phosphate isomerase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,313,5860TC 5.1% 433.5 / 8.3 136F79F (TTT→TTCrpiBallose‑6‑phosphate isomerase/ribose‑5‑phosphate isomerase B
Reads supporting (aligned to +/- strand):  ref base T (64/65);  new base C (3/4);  total (67/69)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.74e-01

Only 100 of 138 total aligned reads displayed.

TGATTGATAAAGGAACCTGGTCGTCAGAGCGTACTGATTATCCACATTACGCCAGTCAAGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACA  >  NC_000913/4313438‑4313732
                                                                                                                                                    |                                                                                                                                                  
tgattgatAAAGGAACCTGGTCGTCAGAGCGTACTGATTATCCACATTACGCCAGTCAAGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTg                                                                                                                                                   <  1:930959/150‑1 (MQ=255)
 gattgatAAAGGAACCTGGTCGTCAGAGCGTACTGATTATCCACATTACGCCAGTCAAGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGc                                                                                                                                                  >  1:581829/1‑150 (MQ=255)
 gattgatAAAGGAACCTGGTCGTCAGAGCGTACTGATTATCCACATTACGCCAGTCAAGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGc                                                                                                                                                  >  1:151359/1‑150 (MQ=255)
   ttgatAAAGGAACCTGGTCGTCAGAGCGTACTGATTATCCACATTACGCCAGTCAAGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCg                                                                                                                                                >  1:1330829/1‑150 (MQ=255)
   ttgatAAAGGAACCTGGTCGTCAGAGCGTACTGATTATCCACATTACGCCAGTCAAGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCAAACTTGTTCGCCg                                                                                                                                                >  2:886781/1‑150 (MQ=255)
    tgatAAAGGAACCTGGTCGTCAGAGCGTACTGATTATCCACATTACGCCAGTCAAGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCgg                                                                                                                                               <  1:1087974/150‑1 (MQ=255)
    tgatAAAGGAACCTGGTCGTCAGAGCGTACTGATTATCCACATTACGCCAGTCAAGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCAAACTTGTTCGCCgc                                                                                                                                               <  1:74118/150‑2 (MQ=255)
     gatAAAGGAACCTGGTCGTCAGAGCGTACTGATTATCCACATTACGCCAGTCAAGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGa                                                                                                                                              >  2:593594/1‑150 (MQ=255)
     gatAAAGGAACCTGGTCGTCAGAGCGTACTGATTATCCACATTACGCCAGTCAAGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGa                                                                                                                                              >  1:780578/1‑150 (MQ=255)
     gatAAAGGAACCTGGTCGTCAGAGCGTACTGATTATCCACATTACGCCAGTCAAGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGa                                                                                                                                              >  1:1141640/1‑150 (MQ=255)
     gatAAAGGAACCTGGTCGTCAGAGCGTACTGATTATCCACATTACGCCAGTCAAGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCAAACTTGTTCGCCgct                                                                                                                                              <  1:303283/150‑3 (MQ=255)
       tAAAGGAACCTGGTCGTCAGAGCGTACTGATTATCCACATTACGCCAGTCAAGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAAt                                                                                                                                            >  1:605410/1‑150 (MQ=255)
        aaaGGAACCTGGTCGTCAGAGCGTACTGATTATCCACATTACGCCAGTCAAGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCAAACTTGTTCGCCgctatc                                                                                                                                           <  1:246366/150‑6 (MQ=255)
         aaGGAACCTGGTCGTCAGAGCGTACTGATTATCCACATTACGCCAGTCAAGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTc                                                                                                                                          >  1:601646/1‑150 (MQ=255)
          aGGAACCTGGTCGTCAGAGCGTACTGATTATCCACATTACGCCAGTCAAGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTcg                                                                                                                                         <  2:315568/150‑1 (MQ=255)
            gAACCTGGTCGTCAGAGCGTACTGATTATCCACATTACGCCAGTCAAGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTcgcg                                                                                                                                       <  1:508778/150‑1 (MQ=255)
             aaCCTGGTCGTCAGAGCGTACTGATTATCCACATTACGCCAGTCAAGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTcgcgc                                                                                                                                      <  2:1184125/150‑1 (MQ=255)
                  ggTCGTCAGAGCGTACTGATTATCCACATTACGCCAGTCAAGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGgtcg                                                                                                                                 >  2:250499/1‑150 (MQ=255)
                  ggTCGTCAGAGCGTACTGATTATCCACATTACGCCAGTCAAGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGgtcg                                                                                                                                 >  1:697693/1‑150 (MQ=255)
                   gtcgtcAGAGCGTACTGATTATCCACATTACGCCAGTCAAGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGgtcgt                                                                                                                                >  1:582569/1‑150 (MQ=255)
                    tcgtcAGAGCGTACTGATTATCCACATTACGCCAGTCAAGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTATCGCGgtcgtc                                                                                                                               <  2:1343183/150‑1 (MQ=255)
                        cAGAGCGTACTGATTATCCACATTACGCCAGTCAAGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTa                                                                                                                           >  1:699043/1‑150 (MQ=255)
                        cAGAGCGTACTGATTATCCACATTACGCCAGTCAAGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTa                                                                                                                           <  1:454402/150‑1 (MQ=255)
                                aCTGATTATCCACATTACGCCAGTCAAGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCt                                                                                                                   <  2:176225/150‑1 (MQ=255)
                                 cTGATTATCCACATTACGCCAGTCAAGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCtt                                                                                                                  >  1:1062/1‑150 (MQ=255)
                                  tGATTATCCACATTACGCCAGTCAAGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTa                                                                                                                 <  2:353369/150‑1 (MQ=255)
                                       aTCCACATTACGCCAGTCAAGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCcg                                                                                                            >  1:859642/1‑150 (MQ=255)
                                       aTCCACATTACGCCAGTCAAGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCcg                                                                                                            >  1:1177551/1‑150 (MQ=255)
                                        tCCACATTACGCCAGTCAAGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCcgc                                                                                                           >  1:1118368/1‑150 (MQ=255)
                                          cacaTTACGCCAGTCAAGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCgcgc                                                                                                         >  1:831896/1‑150 (MQ=255)
                                          cacaTTACGCCAGTCAAGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCgcgc                                                                                                         >  2:722744/1‑150 (MQ=255)
                                          cacaTTACGCCAGTCAAGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCgcgc                                                                                                         >  2:781712/1‑150 (MQ=255)
                                           acaTTACGCCAGTCAAGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCa                                                                                                        <  1:154531/150‑1 (MQ=255)
                                            caTTACGCCAGTCAAGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCaa                                                                                                       >  1:1192637/1‑150 (MQ=255)
                                            caTTACGCCAGTCAAGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCaa                                                                                                       >  1:907565/1‑150 (MQ=255)
                                               tACGCCAGTCAAGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACtt                                                                                                    >  1:192124/1‑150 (MQ=255)
                                                     aGTCAAGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCgg                                                                                              >  1:390653/1‑150 (MQ=255)
                                                         aaGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGgcagc                                                                                          <  2:877374/150‑1 (MQ=255)
                                                          aGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGgcagca                                                                                         >  2:674162/1‑150 (MQ=255)
                                                           gTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCAt                                                                                        >  2:521347/1‑150 (MQ=255)
                                                           gTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGAAAGTTGCGCggaataag                                                                                        <  2:605410/150‑7 (MQ=255)
                                                                   ggCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGacac                                                                                <  2:1184592/150‑1 (MQ=255)
                                                                    gCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACAcc                                                                               <  1:564920/150‑1 (MQ=255)
                                                                    gCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACAcc                                                                               <  2:926488/150‑1 (MQ=255)
                                                                    gCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACAcc                                                                               <  2:660754/150‑1 (MQ=255)
                                                                     cTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGTCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCa                                                                              >  2:827587/1‑150 (MQ=255)
                                                                       gTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAAc                                                                            <  1:323359/150‑1 (MQ=255)
                                                                         tGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGt                                                                          >  1:1408159/1‑150 (MQ=255)
                                                                          gCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTg                                                                         <  2:1071653/150‑1 (MQ=255)
                                                                           cTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGc                                                                        <  2:1411788/150‑1 (MQ=255)
                                                                            tGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCt                                                                       <  2:780578/150‑1 (MQ=255)
                                                                              gCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTTCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTgg                                                                     >  1:1387010/1‑150 (MQ=255)
                                                                               cGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGc                                                                    <  2:1192637/150‑1 (MQ=255)
                                                                                  agagGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCttt                                                                 >  2:235709/1‑150 (MQ=255)
                                                                                   gagGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCtttt                                                                <  2:1142162/150‑1 (MQ=255)
                                                                                   gagGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCtttt                                                                <  2:1305549/150‑1 (MQ=255)
                                                                                     ggTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTgg                                                              <  1:771129/150‑1 (MQ=255)
                                                                                     ggTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTgg                                                              >  2:963469/1‑150 (MQ=255)
                                                                                      gTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGt                                                             >  1:649858/1‑150 (MQ=255)
                                                                                         gATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGATCCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTTGttga                                                          <  1:250499/150‑3 (MQ=255)
                                                                                         gATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCa                                                          <  2:601646/150‑1 (MQ=255)
                                                                                           tGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTTGCTTTTGGTTCACg                                                        <  2:1177551/150‑1 (MQ=255)
                                                                                                    tttGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTgg                                               <  2:1141640/150‑1 (MQ=255)
                                                                                                     ttGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGc                                              >  1:969351/1‑150 (MQ=255)
                                                                                                       gATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCt                                            >  1:85363/1‑150 (MQ=255)
                                                                                                        aTTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTc                                           <  1:674162/150‑1 (MQ=255)
                                                                                                        aTTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTc                                           <  1:673114/150‑1 (MQ=255)
                                                                                                          ttGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGGCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGa                                         <  1:885560/150‑1 (MQ=255)
                                                                                                            gtgGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTTCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAAc                                       <  2:1387010/150‑1 (MQ=255)
                                                                                                              ggTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTg                                     >  2:528203/1‑150 (MQ=255)
                                                                                                              ggTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTg                                     <  2:1064317/150‑1 (MQ=255)
                                                                                                               gTACTGGCGTCGGTATTTCGCTAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTgg                                    <  2:149398/150‑1 (MQ=255)
                                                                                                                tACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGc                                   <  2:192124/150‑1 (MQ=255)
                                                                                                                tACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGc                                   <  1:937103/150‑1 (MQ=255)
                                                                                                                    ggCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCaaaa                               >  1:547554/1‑150 (MQ=255)
                                                                                                                    ggCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCaaaa                               >  2:713411/1‑150 (MQ=255)
                                                                                                                    ggCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAAACTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCaaaa                               >  2:285304/1‑150 (MQ=255)
                                                                                                                     gCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCaaaaa                              <  2:243447/150‑1 (MQ=255)
                                                                                                                         cGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGAt                          <  1:983115/150‑1 (MQ=255)
                                                                                                                          ggTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGAtt                         >  1:465114/1‑150 (MQ=255)
                                                                                                                          ggTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGAGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGAtt                         <  1:101561/150‑1 (MQ=255)
                                                                                                                             aTTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATtgtg                      <  1:711023/150‑1 (MQ=255)
                                                                                                                              tttCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTgg                     <  2:969351/150‑1 (MQ=255)
                                                                                                                               ttCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGa                    >  2:1198178/1‑150 (MQ=255)
                                                                                                                               ttCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGa                    <  2:859642/150‑1 (MQ=255)
                                                                                                                                tCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGAt                   >  2:102204/1‑150 (MQ=255)
                                                                                                                                  aattccggcaAACTTGTTCGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGc                 >  1:1301122/6‑150 (MQ=255)
                                                                                                                                   aTAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCg                >  1:1254867/1‑150 (MQ=255)
                                                                                                                                    ttccggcaAACTTGTTCGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGt               >  1:1251389/4‑150 (MQ=255)
                                                                                                                                        ggcgAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCt           >  2:349400/1‑150 (MQ=255)
                                                                                                                                        ggcgAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCt           >  1:1032670/1‑150 (MQ=255)
                                                                                                                                          cgAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTgg         >  1:398618/1‑150 (MQ=255)
                                                                                                                                          cgAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTgg         <  2:1254867/150‑1 (MQ=255)
                                                                                                                                            aaCAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGgc       >  2:716550/1‑150 (MQ=255)
                                                                                                                                             aCAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGgcg      <  1:346505/150‑1 (MQ=255)
                                                                                                                                               aaGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGGTCACGAGTGGGTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGca    <  2:897059/150‑1 (MQ=255)
                                                                                                                                                aGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGcac   <  1:877628/150‑1 (MQ=255)
                                                                                                                                                 gTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGcaca  >  2:1204308/1‑150 (MQ=255)
                                                                                                                                                  tttGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTgg                     >  1:549401/1‑130 (MQ=255)
                                                                                                                                                  tttGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTgg                     <  2:549401/130‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TGATTGATAAAGGAACCTGGTCGTCAGAGCGTACTGATTATCCACATTACGCCAGTCAAGTCGCACTGGCTGTTGCTGGCGGAGAGGTTGATGGCGGGATTTTGATTTGTGGTACTGGCGTCGGTATTTCGATAGCGGCGAACAAGTTTGCCGGAATTCGCGCGGTCGTCTGTAGCGAACCTTATTCCGCGCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACA  >  NC_000913/4313438‑4313732

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