Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 926,180 A→G 11.1% intergenic (‑209/+45) serW ← / ← infA tRNA‑Ser/translation initiation factor IF‑1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913926,1800AG11.1% 354.2 / 28.2 117intergenic (‑209/+45)serW/infAtRNA‑Ser/translation initiation factor IF‑1
Reads supporting (aligned to +/- strand):  ref base A (57/47);  new base G (7/6);  total (64/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.28e-02

Only 100 of 117 total aligned reads displayed.

ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAATTTTAATGATAATTATTGTTTGCTCATACTGTAAACAAGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGTCGCCCGTCAGGATGCGGATGTAGTTTTTGCGCATTTTACCGGAGATGTGTGCAGTAACCACGTGACCGTTTTCTAACTCTACGCGGAACATGGTATTAGGCAACGTTTCAAGAACGGTACCTTGCATTTCAATATTGTCTTCTTTGGCC  >  NC_000913/925931‑926441
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                     
            tGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAATTTTAATGATAATTATTGTTTGCTCATACTGTAAACAAGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCgg                                                                                                                                                                                                                                   >  4:69935/1‑274 (MQ=255)
            tGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAATTTTAATGATAATTATTGTTTGCTCATACTGTAAACAAGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCgg                                                                                                                                                                                                                                   >  2:449127/1‑274 (MQ=255)
                                                   gCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAATTTTAATGATAATTATTGTTTGCTCATACTGTAAACAAGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGc                                                                                                                                                                                           <  4:449126/275‑1 (MQ=255)
                                                                                                                      aTTTTAATGATAATTATTGTTTGCTCATACTGTAAACAAGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGTCGCCCGTCAGGATGCGGATGTAGTTTTTgcg                                                                                                                        <  2:70008/275‑1 (MQ=255)
                                                                                                                                         tttGCTCATACTGTAAACAAGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGAGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTAGGGGGTCAGTTCAACAGTCACTTGTACCCCCGTCAGGTGGGGGAGGNAGGTTNNGGGCAgttttaaggtgaggtgg                                                                                                     >  3:9654/1‑262 (MQ=255)
                                                                                                                                           tGCTCATACTGTAAACAAGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGTCGCCCGTCAGGATGCGGATGTAGTTTTTGCGCATTTTACCGGAGAtgtgtg                                                                                                    >  2:147638/1‑274 (MQ=255)
                                                                                                                                           tGCTCATACTGTAAACAAGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGAGACTACGGAAAACAATGCGGCCTTTGCTCAGGTCGGACGGGGTCAGTTCAAAAGTCACTTTTCCGCCCGTCCGGATGCGGCTTTAATTTTTGCGCATTTTACCGgagcttgttg                                                                                                    >  4:80760/1‑267 (MQ=255)
                                                                                                                                               cATACTGTAAACAAGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGtt                                                                                                                                                                                                                                                          >  8:520949/1‑120 (MQ=255)
                                                                                                                                               cATACTGTAAACAAGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGtt                                                                                                                                                                                                                                                          <  7:520949/120‑1 (MQ=255)
                                                                                                                                                    tGTAAACAAGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGttcttt                                                                                                                                                                                                                                                      <  8:216545/119‑1 (MQ=255)
                                                                                                                                                     gTAAACAAGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTaaa                                                                                                                                                                                                                               <  8:575161/141‑1 (MQ=255)
                                                                                                                                                     gTAAACAAGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTaaa                                                                                                                                                                                                                               >  7:575161/1‑141 (MQ=255)
                                                                                                                                                           aaGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGTCGCCCGTCAGGATGCGGATGTAGTTTTTGCGCATTTTACCGGAGATGTGTGCAGTAACCACGTGACCg                                                                                   <  4:147637/275‑1 (MQ=255)
                                                                                                                                                                gtgCAGTATATCTACATCGAGACAAGTTACGGCCTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGCCAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTGTTCGCACGTCAGGCTTCGGATGTAGTTTTTCAGCATTTTACCGTACATGTTTGAAGTACCCCCGCTACCGTTTTc                                                                              >  4:358526/1‑275 (MQ=255)
                                                                                                                                                                     gTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGTCGCCCGTCAGGATGCGGATGTAGTTTTTGCGCATTTTACCGGAGATGTGTGCAGTAACCACGTGACCGTTTTCTAAct                                                                         <  2:358527/275‑1 (MQ=255)
                                                                                                                                                                           cTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTAAGGTCGTACGGGGTCAGGTATACAGTAACTTTGTCGCCCGTCAGGCTGCGGCTGTAGTTTTTGCGGATTTTACCGGAGCTGTGTGCAGTAACCACGGCACACGTGTCTAACTCTatta                                                                   >  4:536778/1‑272 (MQ=255)
                                                                                                                                                                             aCATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGttc                                                                                                                                                                                                                                                         >  8:928450/1‑91 (MQ=255)
                                                                                                                                                                             aCATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGttc                                                                                                                                                                                                                                                         <  7:928450/91‑1 (MQ=255)
                                                                                                                                                                              cATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGTCGCCCGTCAGGATGCGGATGTAGTTTTTGCGCATTTTACCGGAGATGTGTGCAGTAACCACGTGACCGTTTTCTAACTCTACGCGGa                                                                >  2:300143/1‑275 (MQ=255)
                                                                                                                                                                              cATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGTCGCCCGTCAGGATGCGGATGTAGTTTTTGCGCATTTTACCGGAGATGTGTGCAGTAACCACGTGACCGTTTTCTAACTCTACGCGGa                                                                >  2:476118/1‑275 (MQ=255)
                                                                                                                                                                                tCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGccc                                                                                                                                                                                                                                            >  2:276017/1‑101 (MQ=255)
                                                                                                                                                                                tCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCa                                                                                                                                                                                                                                           <  8:915469/102‑1 (MQ=255)
                                                                                                                                                                                tCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCa                                                                                                                                                                                                                                           >  7:915469/1‑102 (MQ=255)
                                                                                                                                                                                tCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTCCGGGGTCAGTTAAACAGTCACTTTGTCGCCCGTCCGGATGCGGATGTAGTTTTTGCGCATTTTACCGGAGTTTTGTGAAGTAACCAGGTGCCCGTTTTCTAAATCTACGCGGaa                                                               >  4:3156/1‑274 (MQ=255)
                                                                                                                                                                                tCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGTCGCCCGTAAGGATGCGGATGTAGTTTTTGCGCATTTTACCGGAGATGTGTGCAGTAACCACGTGACCGTTTTCTAACTCTACGCGGaa                                                               >  2:172205/1‑274 (MQ=255)
                                                                                                                                                                                  gagaCAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGTCGCCCGTCAGGATGCGGATGTAGTTTTTGCGCATTTTACCGGAGATGTGTGCAGTAACCACGTGACCGTTTTCTAACTCTACGCGGAACa                                                             <  2:69936/274‑1 (MQ=255)
                                                                                                                                                                                    gaaaaGTTACGGACTTATACTGCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTTGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGt                                                                                                                                                                           <  2:98198/159‑1 (MQ=255)
                                                                                                                                                                                    gaCAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCa                                                                                                                                                                                                                                        <  4:276016/101‑1 (MQ=255)
                                                                                                                                                                                    gaCAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCa                                                                                                                                                                                                                        >  1:80419/1‑117 (MQ=255)
                                                                                                                                                                                     ccAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATAACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGTCGCCCGTCAGGATGCGGATGTAGTTTTTGCGCATTTTACCGGAGATGTGTGCAGTAACCACGTGACCGTTTTCTAACTCTATGCGGAACa                                                             <  2:536779/270‑1 (MQ=255)
                                                                                                                                                                                      cAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATc                                                                                                                                                                                                                         >  5:992982/1‑114 (MQ=255)
                                                                                                                                                                                      cAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATc                                                                                                                                                                                                                         <  6:992982/114‑1 (MQ=255)
                                                                                                                                                                                        aGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGa                                                                                                                                                                                                                    <  3:80419/117‑1 (MQ=255)
                                                                                                                                                                                                    ataCTTCCAAAGTACTTCATACATATAACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGTCGCCCGTCAGGATGCGGATGTAGTTTTTGCGCATTTTACCGGAGATGTGTGCAGTAACCACGTGACCGTTTTCTAACTCTACGCGGAACATGGTATTAGGCAACGTTTc                                          <  1:9654/275‑1 (MQ=255)
                                                                                                                                                                                                         tCCAAAGTACTTCATACATAACACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGTCGCCCGTCAGGATGCGGATGTAGTTTTTGCGCATTTTACCGGAGATGTGTGCAGTAACCACGTGACCGTTTTCTAACTCTACGCGGAACATGGTATTAGGCAACGTTTCAAGaa                                     <  2:3156/275‑1 (MQ=255)
                                                                                                                                                                                                          ccAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGTCGCCCGTCAGGATGCGGATGTAGTTTTTGCGCATTTTACCGGAGATGTGTGCAGTAACCACGTGACCGTTTTCTAACTCTACGCgg                                                                 >  1:7787/1‑246 (MQ=255)
                                                                                                                                                                                                          ccAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCGCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGTCGCCCGTCAGGATGCGGCTGTAGTTTTTGCTCATTTTACCGGAGCTGGGTGCAGTCACCACGTGACCGTTTTCTAACTCTCAGCGGAACATGGTATTAGGCAACGTTTCAAGAAc                                    >  2:381377/1‑275 (MQ=255)
                                                                                                                                                                                                              aGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGt                                                                                                                                                               <  7:343852/148‑1 (MQ=255)
                                                                                                                                                                                                              aGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGt                                                                                                                                                               >  8:343852/1‑148 (MQ=255)
                                                                                                                                                                                                                  cTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGTCGCCCGTCAGGATGCGGATGTAGTTTTTGCGCATTTTACCGGAGATGTGTGCAGTAACCACGTGACCGTTTTCTAACTCTACGCGGAACATGGTATTAGGCAACGTTTCAAGAACGGTACCtt                            >  2:97797/1‑275 (MQ=255)
                                                                                                                                                                                                                   ttaatacataTAAAAAAATAAAAAGGCCGGTTACACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCGTGCGGTANAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGTCGCCCGTCAGGATGCGGATGTAGTTTTTGCGCATTTTACCGGAGATGTGTGCAGTAACCACGTGACCGTTTTCTAACTCTACGCGGAACATGGTATTAGGCAACGTTTCAAGAACGGTACCTTg                           <  3:127437/272‑1 (MQ=255)
                                                                                                                                                                                                                       tacataTCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGa                                                                                                                                                                                                             <  7:902519/93‑1 (MQ=255)
                                                                                                                                                                                                                       tacataTCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGa                                                                                                                                                                                                             >  8:902519/1‑93 (MQ=255)
                                                                                                                                                                                                                          atatCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGTCGCCCGTCAGGATGCGGATGTAGTTTTTGCGCATTTTACCGGAGATGTGTGCAGTAACCACGTGACCGTTTTCTAACTCTACGCGGAACATGGTATTAGGCAACGTTTCAAGAACGGTACCTTGCATTTCa                    >  4:391950/1‑275 (MQ=255)
                                                                                                                                                                                                                            atCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCa                                                                                                                                                                                                                        >  1:122744/1‑77 (MQ=255)
                                                                                                                                                                                                                            atCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGTCGCCCGTCAGGATGCGGATGTAGTTTTTGCGCATTTTACCGGAGATGTGTGCAGTAACCACGTGACCGTTTTCTAACTCTACGCGGAACATGGTATTAGGCAACGTTTc                                          <  2:424151/251‑1 (MQ=255)
                                                                                                                                                                                                                            atCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGTCGCCCGTCAGGATGCGGATGTAGTTTTTGCGCATTTTACCGGAGATGTGTGCAGTAACCACGTGACCGTTTTCTAACTCTACGCGGAACATGGTATTAGGCAACGTTTCAAGAACGTTAc                               >  1:127437/1‑262 (MQ=255)
                                                                                                                                                                                                                            atCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGTCGCCCGTCAGGATGCGGATGTAGTTTTTGCGCATTTTACCGGAGATGTGTGCAGTAACCACGTGACCGTTTTCTAACTCTACGCGGAACATGGTATTAGGCAACGTTTCAAGAACGGTACCt                             >  4:205775/1‑264 (MQ=255)
                                                                                                                                                                                                                            atCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGTCGCCCGTCAGGATGCGGATGTAGTTTTTGCGCATTTTACCGGAGATGTGTGCAGTAACCACGTGACCGTTTTCTAACTCTACGCGGAACATGGTATTAGGCAACGTTTCAAGAACGGTACCTTGCATTTCaa                   >  2:364549/1‑274 (MQ=255)
                                                                                                                                                                                                                              cacaAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAAc                                                                                                                                                                  >  7:327175/1‑129 (MQ=255)
                                                                                                                                                                                                                                caAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCgg                                                                                                                                                                                                                                   <  8:159990/62‑1 (MQ=255)
                                                                                                                                                                                                                                caAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCgg                                                                                                                                                                                                                                   >  7:159990/1‑62 (MQ=255)
                                                                                                                                                                                                                                caAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCgg                                                                                                                                                                                                                                   >  7:106917/1‑62 (MQ=255)
                                                                                                                                                                                                                                caAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCgg                                                                                                                                                                                                                                   <  8:106917/62‑1 (MQ=255)
                                                                                                                                                                                                                                caAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCATCGa                                                                                                                                                                                                                    <  3:122744/77‑1 (MQ=255)
                                                                                                                                                                                                                                caAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCttt                                                                                                                                                                                             <  6:270672/100‑1 (MQ=255)
                                                                                                                                                                                                                                caAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCttt                                                                                                                                                                                             >  5:270672/1‑100 (MQ=255)
                                                                                                                                                                                                                                caAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGTCGCCCGTCAGGATGCGGATGTAGTTTTTGCGCATTTTACCGGAGATGTGTGCAGTAACCACGTGACCGTTTTCTAACTCTACGCGGAACATGGTATTAGGCAACGTTTCAAGAACGGTACCTTGCa                         <  2:205776/264‑1 (MQ=255)
                                                                                                                                                                                                                                 aaaaTAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGTCGCCCGTCAgg                                                                                                                                            >  5:775260/1‑148 (MQ=255)
                                                                                                                                                                                                                                  aaaTAAAAAGGCCGGTTAAACCGACCTTTTACTCGt                                                                                                                                                                                                                                                           >  8:740667/1‑36 (MQ=255)
                                                                                                                                                                                                                                  aaaTAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACggg                                                                                                                                                                             <  8:353003/114‑1 (MQ=255)
                                                                                                                                                                                                                                   aaTAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGTCGCCCGTCAGGATGCGGATGTAGTTTTTGCGCATTTTACCGGAGATGTGTGCAGTAACCACGTGACCGTTTTCTAACTCTACGCgg                                                                 >  2:507921/1‑221 (MQ=255)
                                                                                                                                                                                                                                    aTAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGTCGCCCGTCAGGATg                                                                                                                                         <  6:775260/148‑1 (MQ=255)
                                                                                                                                                                                                                                    aTAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGTCGCCCGTCAGGATg                                                                                                                                         >  7:528917/1‑148 (MQ=255)
                                                                                                                                                                                                                                     tAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATg                                                                                                                                                                                                     >  7:1095631/1‑87 (MQ=255)
                                                                                                                                                                                                                                     tAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGt                                                                                                                                                       >  8:531771/1‑133 (MQ=255)
                                                                                                                                                                                                                                     tAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGt                                                                                                                                                       <  7:531771/133‑1 (MQ=255)
                                                                                                                                                                                                                                      taaaaGGTCGGTTTAACCGGCCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCATTTCAACAGTCACTTTg                                                                                                                                                        >  5:82041/2‑131 (MQ=255)
                                                                                                                                                                                                                                      aaaaaGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACgg                                                                                                                                                                              >  8:1090976/1‑109 (MQ=255)
                                                                                                                                                                                                                                      aaaaaGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACgg                                                                                                                                                                              <  7:1090976/109‑1 (MQ=255)
                                                                                                                                                                                                                                       aaaaGGTCGGTTTAACCGGCCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGc                                                                                                                                                                                                                                     >  5:1105657/1‑53 (MQ=25)
                                                                                                                                                                                                                                       aaaaGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGTCGCCCGTCAGGATGCGGATGTAGTTTTTGCGCATTTTACCGGAGATGTGTGCAGTAACCACGTGACCGTTTTCTAACTCTACGCGGAACa                                                             <  4:507920/221‑1 (MQ=255)
                                                                                                                                                                                                                                        aaaGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGTCGCCCGTCAGGATGCGGATGTAGTTTTTGCGCATTTTAc                                                                                                                >  4:526637/1‑169 (MQ=255)
                                                                                                                                                                                                                                         aaGGTCGGTTTAACCGGCCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGc                                                                                                                                                                                                 <  5:1046275/87‑1 (MQ=39)
                                                                                                                                                                                                                                         aaGGTCGGTTTAACCGGCCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGc                                                                                                                                                                                                 <  5:744247/87‑1 (MQ=39)
                                                                                                                                                                                                                                         aaGGTCGGTTTAACCGGCCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGc                                                                                                                                                                                                 >  6:1046275/1‑87 (MQ=39)
                                                                                                                                                                                                                                         aaGGTCGGTTTAACCGGCCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGc                                                                                                                                                                                                 >  6:744247/1‑87 (MQ=39)
                                                                                                                                                                                                                                          aggtcGGTTTAACCGGCCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAAc                                                                                                                                                                  <  6:190810/113‑1 (MQ=255)
                                                                                                                                                                                                                                          aGGTCGGTTTAACCGGCCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAg                                                                                                                                                                                                           >  7:291480/1‑76 (MQ=38)
                                                                                                                                                                                                                                          aGGTCGGTTTAACCGGCCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAAc                                                                                                                                                                  >  5:190810/1‑117 (MQ=255)
                                                                                                                                                                                                                                          aGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGTCGCCCGTCAGGATGCGGATGTAGTTTTTGCGCATTTTACCGGAGATGTGTGCAGTAACCACGTGACCGTTTTCTAACTCTACGCGGAACATGGTATTAGGCAACGTTTCAAGAACGGTACCTTGCATTTCAATATTGTCTTCTTTgg    <  2:434308/275‑1 (MQ=255)
                                                                                                                                                                                                                                           ggtcGGTTTAACCGGCCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCa                                                                                                                                                                    >  7:190118/4‑114 (MQ=255)
                                                                                                                                                                                                                                           ggtcGGTTTAACCGGCCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCa                                                                                                                                                                    <  8:190118/111‑1 (MQ=255)
                                                                                                                                                                                                                                            gCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACaa                                                                                                                                                                                                                           >  6:58648/1‑58 (MQ=255)
                                                                                                                                                                                                                                            gCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACaa                                                                                                                                                                                                                           <  5:58648/58‑1 (MQ=255)
                                                                                                                                                                                                                                            gCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGTCGCCCGTCAGGATGCGGATGTAGTTTTTGCGCATTTTACCGga                                                                                                            <  2:526638/169‑1 (MQ=255)
                                                                                                                                                                                                                                            gCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGTCGCCCGTCAGGATGCGGATGTAGTTTTTGCGCATTTTACCGGAGATGTGTGCAGt                                                                                                <  1:96380/181‑1 (MQ=255)
                                                                                                                                                                                                                                            gCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGTCGCCCGTCAGGATGCGGATGTAGTTTTTGCGCATTTTACCGGAGATGTGTGCAGTAACCACGTGACCGTTTTCTAACTCTACGCGGAACATGGTATTAGGCAACGTTTCAAGAACGGTACCTTGCATTTCAATATTGTCTTCTTTGGcc  <  4:464660/275‑1 (MQ=255)
                                                                                                                                                                                                                                                 ttAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTa                                                                                                                                                                                                                 <  8:619580/63‑1 (MQ=255)
                                                                                                                                                                                                                                                 ttAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTa                                                                                                                                                                                                                 >  7:619580/1‑63 (MQ=255)
                                                                                                                                                                                                                                                 ttAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACgg                                                                                                                                                                              >  7:997656/1‑98 (MQ=255)
                                                                                                                                                                                                                                                 ttAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGTCGCCCGTc                                                                                                                                               >  7:308967/1‑129 (MQ=255)
                                                                                                                                                                                                                                                 ttAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGTCGCCCGTc                                                                                                                                               <  8:308967/129‑1 (MQ=255)
                                                                                                                                                                                                                                                    aaCCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCtt                                                                                                                                                                                              <  6:690712/79‑1 (MQ=255)
                                                                                                                                                                                                                                                      ccGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGTCGCCCGTCAGGATGCGGa                                                                                                                                     >  7:727845/1‑134 (MQ=255)
                                                                                                                                                                                                                                                        gACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGc                                                                                                                                                                                           >  8:323586/1‑78 (MQ=255)
                                                                                                                                                                                                                                                        gACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGc                                                                                                                                                                                           <  7:323586/78‑1 (MQ=255)
                                                                                                                                                                                                                                                        gACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGt                                                                                                                                                                           >  8:292203/1‑94 (MQ=255)
                                                                                                                                                                                                                                                        gACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGt                                                                                                                                                                           <  7:292203/94‑1 (MQ=255)
                                                                                                                                                                                                                                                        gACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGTCGCCCGTCAGGATGCGGATGTAGTTTTTGCGCAt                                                                                                                     >  8:964010/1‑148 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                     
ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAATTTTAATGATAATTATTGTTTGCTCATACTGTAAACAAGTTGTGCAGTATATCTACATCGAGACAAGTTACGGACTTATACTTCCAAAGTACTTCATACATATCACAAAATAAAAAGGCCGGTTAAACCGACCTTTTACTCGTTCTTTCTCTTCGCCCATCAGGCGGTAAAACAATCAGCGACTACGGAAGACAATGCGGCCTTTGCTCAGGTCGTACGGGGTCAGTTCAACAGTCACTTTGTCGCCCGTCAGGATGCGGATGTAGTTTTTGCGCATTTTACCGGAGATGTGTGCAGTAACCACGTGACCGTTTTCTAACTCTACGCGGAACATGGTATTAGGCAACGTTTCAAGAACGGTACCTTGCATTTCAATATTGTCTTCTTTGGCC  >  NC_000913/925931‑926441

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