Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 566,651 T→C 100% pseudogene (276/312 nt) peaD → pseudogene; DLP12 prophage; replication protein family;Phage or Prophage Related

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913566,6510TC98.5% 8030.9 / 31.4 2809pseudogene (276/312 nt)peaDpseudogene; DLP12 prophage; replication protein family;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base T (20/26);  new base C (1421/1332);  total (1448/1361)
Fisher's exact test for biased strand distribution p-value = 3.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 2809 total aligned reads displayed.

ATTACGGCGCAAGGCTGCCGATGAACTGACCTGTATGACAGCGCGAATTAACCGTGGTGAGACGATACCTGAACCAGTAAAACAACTTCCTGTTATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATTAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATAG  >  NC_000913/566504‑566795
                                                                                                                                                   |                                                                                                                                                
 ttACGCCGTAAGGCCGCAGATGAGCTTGTCCATATGACTGCGAGAATTAACCGTCGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCaa                                                                                                                                                <  5:292074/149‑1 (MQ=255)
                            gtccatgtGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCaa                                                                                                                      <  2:390860/141‑1 (MQ=255)
                            gtccatATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCaaa                                                                                                                                               >  5:185214/6‑123 (MQ=255)
                            gtccatATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAg                                                                                                                     >  2:280690/6‑149 (MQ=255)
                                atATGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAACCAACTTCCTGTCATGGGCGGTAGACCTCTACATCGTGCACAGGCTATGGCGAAGATCGCAGAAATCAAAGCTAAGTTAGGACTGAAAGTAGCAAGTCTa                                                                                                                 >  5:155937/2‑149 (MQ=255)
                                    gACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGa                                                                                                                          >  6:260118/1‑136 (MQ=255)
                                       tgcgaGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACggg                                                                                                          >  7:134720/2‑149 (MQ=255)
                                         cgagAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGc                                                                                                         >  3:334108/4‑148 (MQ=255)
                                         cgagAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGc                                                                                                         <  8:234580/145‑1 (MQ=255)
                                          gagAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCaa                                                                                                       <  7:228763/147‑1 (MQ=255)
                                          gagAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCa                                                                                                        >  4:170995/3‑148 (MQ=255)
                                              aTTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGCTCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAagag                                                                                                   >  4:136673/1‑149 (MQ=255)
                                               ttAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAgg                                                                                                  <  6:383330/149‑1 (MQ=255)
                                                tAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGCCCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGGTCGGACTGAACGGCGCAAGTGAATGCCGGGCAAAGAgg                                                                                                  >  6:198383/1‑148 (MQ=255)
                                                 aaacGTGGTGAGGAGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGc                                                                                                 <  6:99962/145‑1 (MQ=255)
                                                 aaCTGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGc                                                                                                 <  6:105752/148‑1 (MQ=255)
                                                  aCCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGCTCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCaa                                                                                               >  3:52501/1‑149 (MQ=255)
                                                  aCCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCaa                                                                                               <  1:101678/149‑1 (MQ=255)
                                                  aCCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCaa                                                                                               <  6:163842/149‑1 (MQ=255)
                                                   ccGTGGTGAGACGATACCTGAACCAGTAAAACAACTTCCTGTTATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATTAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAat                                                                                              >  3:225500/1‑149 (MQ=255)
                                                        gtgAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAattatt                                                                                          <  7:182957/148‑1 (MQ=255)
                                                        gtgAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTc                                                                                         <  6:285700/149‑1 (MQ=255)
                                                        gggaggCGATCCCGGAACCAGTAAAACAACTTCCTGTCCTGGGCGGTAGACCTCTAAATCGTGCACAGGCGCTGGCGAAGATCGCAGAAAGCCAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTc                                                                                         <  4:95666/147‑1 (MQ=255)
                                                         tgAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCaaa                                                                                                                                               <  8:344712/94‑1 (MQ=255)
                                                             gcGATCCCTGAACCAGTAAAACAACTTCCTGCCATGGGCGGTAGACCTCTAAATCGTGCCCAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTAc                                                                                    <  8:370442/148‑1 (MQ=255)
                                                              cGATACCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACATGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTAcc                                                                                   <  1:76568/149‑1 (MQ=255)
                                                               gATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCt                                                                                  <  7:281025/149‑1 (MQ=255)
                                                               gATCCCCGAACCAGTAAAACAACTTCCTGTCATTGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCt                                                                                  >  5:325919/1‑149 (MQ=255)
                                                                atcccTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACCGAAAGGAGCaa                                                                                                                      <  3:277417/109‑1 (MQ=255)
                                                                 tcccTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTgg                                                                                <  5:266364/147‑1 (MQ=255)
                                                                   ccTGAACCAGTAAAACAACTTCCTGTCATGGGTGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTc                                                                                         <  5:3591/138‑1 (MQ=255)
                                                                    cTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTAc                                                                                    <  1:99132/142‑1 (MQ=255)
                                                                    cTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGa                                                                             >  4:352067/1‑149 (MQ=255)
                                                                     tGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGAc                                                                            >  7:35904/1‑149 (MQ=255)
                                                                       aaCCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGc                                                                          <  1:46730/149‑1 (MQ=255)
                                                                          cAGTAAAACAACTTCCTGTCATGGGCTGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTGCGGACTGCAAGGCGAAAGTGTATGCAGGGCAAAGAGGCAATTATTCATGACCTGTGGACGCAt                                                                        >  4:45562/1‑148 (MQ=255)
                                                                          cAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCTCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCAt                                                                        >  8:119049/1‑148 (MQ=255)
                                                                          cAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCAt                                                                        >  7:300509/1‑148 (MQ=255)
                                                                          cAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGAGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCAt                                                                        >  1:199724/1‑148 (MQ=255)
                                                                           aGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATaa                                                                      <  4:322942/149‑1 (MQ=255)
                                                                            gTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATaa                                                                      >  8:399274/1‑148 (MQ=255)
                                                                            gTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATaa                                                                      >  1:350762/1‑148 (MQ=255)
                                                                             tAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATaata                                                                    <  3:357257/149‑3 (MQ=255)
                                                                             tAAAACAACATCCTGTCATGGGCGGTTGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGCAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATaa                                                                      >  4:203746/1‑147 (MQ=255)
                                                                               aaaCAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAatag                                                                   <  7:154068/148‑1 (MQ=255)
                                                                               aaaCAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAatag                                                                   <  8:260337/148‑1 (MQ=255)
                                                                               aaaCAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGc                                                                  >  5:382953/1‑149 (MQ=255)
                                                                                aaaaacTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGCCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACCGGCAAAGAGGCAATTATTCATTACCGGGGGACGCATAATAGCt                                                                 >  1:235950/4‑149 (MQ=255)
                                                                                 acaacTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGc                                                                                                 <  5:209044/116‑1 (MQ=255)
                                                                                 acaacTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCtt                                                                <  1:374636/149‑1 (MQ=255)
                                                                                  caacTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTc                                                               >  7:375749/1‑149 (MQ=255)
                                                                                   aacTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCt                                                              <  3:141854/149‑1 (MQ=255)
                                                                                   aacTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCt                                                              >  8:207855/1‑149 (MQ=255)
                                                                                    acTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCtg                                                             >  8:147040/1‑149 (MQ=255)
                                                                                     cTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACTCATAATAGCTTCtgt                                                            >  5:236629/1‑149 (MQ=255)
                                                                                     cTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCtgt                                                            >  6:129088/1‑149 (MQ=255)
                                                                                      ttCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCAGAATAGCTTCtgtg                                                           >  5:115637/1‑149 (MQ=255)
                                                                                       tCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCCAAGATCGCAGAAAACAAAGCTAAGTTCGGACTGAAAGCAGCAAGTGTAAGACGGGCAAaga                                                                                                    >  8:310921/1‑107 (MQ=255)
                                                                                         cTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGcgc                                                        >  2:347980/1‑147 (MQ=255)
                                                                                          tGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACACAAAATAGCTTCTGTGcgcc                                                       >  3:111968/1‑146 (MQ=255)
                                                                                           gtcaTGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTg                                                                                 >  4:288125/4‑122 (MQ=255)
                                                                                              aTGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGcgc                                                        <  7:5040/144‑3 (MQ=255)
                                                                                              aTGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGAAGAAATCAAAGCTAAGTGAGGACTGAATGGAGCAAGAGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCAGAATAGCTTCTGTGCGCCGGAc                                                   >  6:226799/1‑149 (MQ=255)
                                                                                                gggCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAgtgt                                                                                                                  <  2:25412/84‑1 (MQ=255)
                                                                                                    ggTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGcc                                             <  2:177412/149‑1 (MQ=255)
                                                                                                      tAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCcgc                                           <  1:288110/149‑1 (MQ=255)
                                                                                                        gACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCcgcgc                                         <  3:308150/149‑2 (MQ=255)
                                                                                                          ccTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCcgcgct                                        >  1:12737/1‑146 (MQ=255)
                                                                                                           ctctAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAaga                                                                                                    >  2:315890/1‑87 (MQ=255)
                                                                                                           ctctAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTaa                                      >  8:285789/1‑149 (MQ=255)
                                                                                                           ctctAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTaa                                      >  4:106464/1‑149 (MQ=255)
                                                                                                             ctAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATaat                                                                     <  7:329054/116‑2 (MQ=255)
                                                                                                             ctAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTaa                                      >  2:201775/1‑147 (MQ=255)
                                                                                                             ctAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTaa                                      >  1:258731/1‑147 (MQ=255)
                                                                                                              tAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAg                                   <  1:149329/149‑1 (MQ=255)
                                                                                                                aaTCGTGCCCAGGCTCTGGCGAAGATCGCAGAAATCAACGCTAAGTTCGGACGGAAAGGAGCAAGTGTAGGACGGGCAAAGAGGCCATTATTCATTACCTGGGGACGCATAATAGCTGCTGTGCGCCGGACGTTGCCGCGCTAACAGgc                                 <  5:85225/149‑1 (MQ=255)
                                                                                                                aaTCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGa                                                                                                                          <  1:318149/60‑1 (MQ=255)
                                                                                                                aaTCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACa                                    <  2:371459/146‑1 (MQ=255)
                                                                                                                aaTCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGgc                                 <  7:208459/149‑1 (MQ=255)
                                                                                                                aaTCGTGCACAGGCTCTGGCGAAGATCGCAGAAAGCAAAGATAAGTTCGGTCTGATAGGAGCAAGTGTATGCCGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGgc                                 <  3:256004/149‑1 (MQ=255)
                                                                                                                    gTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCaa                             <  7:248420/149‑1 (MQ=255)
                                                                                                                       cacaGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGCATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAAca                           <  3:192384/148‑2 (MQ=255)
                                                                                                                        acaGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAAcag                          >  3:167999/1‑146 (MQ=255)
                                                                                                                        acaGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGAAAcag                          >  6:321941/1‑146 (MQ=255)
                                                                                                                          aGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAA                       <  4:228124‑M1/149‑1 (MQ=255)
                                                                                                                           ggCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAc                      <  8:69773‑M1/149‑2 (MQ=255)
                                                                                                                              tctGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccag                   >  4:7501‑M1/1‑145 (MQ=255)
                                                                                                                                tGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAA                       >  6:63068‑M1/1‑143 (MQ=255)
                                                                                                                                  gCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcat                >  3:87974‑M1/1‑141 (MQ=255)
                                                                                                                                  gCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcat                <  2:156529‑M1/148‑8 (MQ=255)
                                                                                                                                    gAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCCACAGTATccagcataaa             >  6:350401‑M1/1‑139 (MQ=255)
                                                                                                                                      aGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaat            >  5:57692‑M1/1‑137 (MQ=255)
                                                                                                                                       gATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTTGGGACGCATAATAGCTTCTGTGCGCCGGACGCTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatca          <  8:173283‑M1/149‑14 (MQ=255)
                                                                                                                                       gATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatca          <  7:85027‑M1/149‑14 (MQ=255)
                                                                                                                                         tCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcagg        <  2:127888‑M1/149‑16 (MQ=255)
                                                                                                                                         tCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcagg        >  2:90757‑M1/1‑134 (MQ=255)
                                                                                                                                           gCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggcc      >  2:224863‑M1/1‑132 (MQ=255)
                                                                                                                                              gAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcg   <  4:255785‑M1/149‑21 (MQ=255)
                                                                                                                                               aaaTCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCaa                                                                                               <  2:60648/56‑1 (MQ=255)
                                                                                                                                               aaaTCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcgg  <  1:176320‑M1/149‑22 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
ATTACGGCGCAAGGCTGCCGATGAACTGACCTGTATGACAGCGCGAATTAACCGTGGTGAGACGATACCTGAACCAGTAAAACAACTTCCTGTTATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATTAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATAG  >  NC_000913/566504‑566795

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