Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 576,998 G→T 100% pseudogene (189/309 nt) aaaD → pseudogene; DLP12 prophage; tail fiber assembly protein family;Phage or Prophage Related

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273576,9980GT98.7% 721.0 / ‑0.8 225pseudogene (189/309 nt)aaaDpseudogene; DLP12 prophage; tail fiber assembly protein family;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base T (114/108);  total (116/109)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 225 total aligned reads displayed.

TCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTC  >  CP009273/576853‑577142
                                                                                                                                                 |                                                                                                                                                
 ccGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATg                                                                                                                                               >  1:702011/1‑148 (MQ=255)
      ccGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAg                                                                                                                                          <  2:209645/148‑1 (MQ=255)
        gggCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGCAAGAATCCGCATACCAGGAAGGGCGCTGGGAACCACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCg                                                                                                                                        >  2:454485/1‑148 (MQ=255)
        gggCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCg                                                                                                                                        >  2:1286428/1‑148 (MQ=255)
        gggCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCg                                                                                                                                        >  2:3153954/1‑148 (MQ=255)
        gggCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCg                                                                                                                                        <  1:1900233/148‑1 (MQ=255)
           cTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACt                                                                                                                                     >  2:619363/1‑148 (MQ=255)
           cTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACt                                                                                                                                     >  1:1710229/1‑148 (MQ=255)
              cGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACa                                                                                                                                  <  1:3490538/148‑1 (MQ=255)
                     tGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGa                                                                                                                           >  1:2244634/1‑148 (MQ=255)
                     tGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGa                                                                                                                           >  2:108157/1‑148 (MQ=255)
                      gTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAg                                                                                                                          <  2:2275862/148‑1 (MQ=255)
                      gTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAg                                                                                                                          >  2:1347177/1‑148 (MQ=255)
                       gtCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAgg                                                                                                                         <  2:2026837/147‑1 (MQ=255)
                          aGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGa                                                                                                                      <  1:3542575/148‑1 (MQ=255)
                              cGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAAtgt                                                                                                                  <  1:2314316/148‑1 (MQ=255)
                               ggAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAAtgtg                                                                                                                 >  2:2848320/1‑148 (MQ=255)
                                gAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATgtgg                                                                                                                >  1:2408658/1‑148 (MQ=255)
                                gAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATgtgg                                                                                                                >  1:2408761/1‑148 (MQ=255)
                                   ccACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGa                                                                                                             >  1:1433448/1‑148 (MQ=255)
                                      caGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGt                                                                                                          <  2:1813892/148‑1 (MQ=255)
                                       aGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTc                                                                                                         >  2:2678192/1‑148 (MQ=255)
                                         aCCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTg                                                                                                       >  2:3355367/1‑148 (MQ=255)
                                         aCCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTg                                                                                                       >  1:3078522/1‑148 (MQ=255)
                                          ccgccgTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGc                                                                                                      <  1:3287935/148‑1 (MQ=255)
                                           cgccgTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGcc                                                                                                     >  2:2943330/1‑148 (MQ=255)
                                              cgTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCgt                                                                                                  <  2:1390153/148‑1 (MQ=255)
                                              cgTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCgt                                                                                                  <  2:3165350/148‑1 (MQ=255)
                                              cgTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCgt                                                                                                  <  2:2443279/148‑1 (MQ=255)
                                               gTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCgtg                                                                                                 <  1:369203/148‑1 (MQ=255)
                                               gTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCgtg                                                                                                 <  1:369973/148‑1 (MQ=255)
                                                ttGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCgtgt                                                                                                <  2:3448303/148‑1 (MQ=255)
                                                  gAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCg                                                                                              <  1:1837298/148‑1 (MQ=255)
                                                    aTGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGt                                                                                            <  1:3277326/148‑1 (MQ=255)
                                                     tGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGtt                                                                                           <  2:2908890/148‑1 (MQ=255)
                                                      gggCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTa                                                                                          <  2:5461/148‑1 (MQ=255)
                                                       ggCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTAt                                                                                         <  1:2632138/148‑1 (MQ=255)
                                                           gATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCa                                                                                     >  1:2161731/1‑148 (MQ=255)
                                                            aTGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCaa                                                                                    >  1:1794215/1‑148 (MQ=255)
                                                             tGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCaaa                                                                                   <  1:2597098/148‑1 (MQ=255)
                                                              gCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCaaaa                                                                                  <  2:1800599/148‑1 (MQ=255)
                                                               cTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAAt                                                                                 >  2:3562119/1‑148 (MQ=255)
                                                                tAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAAtg                                                                                <  1:2329178/148‑1 (MQ=255)
                                                                 aaTTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAAtgc                                                                               <  2:1921421/148‑1 (MQ=255)
                                                                   ttACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAAtgctg                                                                             <  1:1414374/148‑1 (MQ=255)
                                                                       tATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTgg                                                                         >  1:3209977/1‑148 (MQ=255)
                                                                       tATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTgg                                                                         <  1:2993372/148‑1 (MQ=255)
                                                                          ctcCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGgtg                                                                      <  1:1353001/148‑1 (MQ=255)
                                                                           tcCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCgtgt                                                                                                <  2:3505238/121‑1 (MQ=255)
                                                                           tcCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGgtgt                                                                     <  2:2805187/148‑1 (MQ=255)
                                                                             ccGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGATGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGttt                                                                   <  2:2681185/148‑1 (MQ=255)
                                                                             ccGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGttt                                                                   <  2:2279681/148‑1 (MQ=255)
                                                                              ctaaaGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTa                                                                  >  1:1331531/3‑148 (MQ=255)
                                                                              cGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTa                                                                  <  1:89258/148‑1 (MQ=255)
                                                                              cGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTa                                                                  <  1:434762/148‑1 (MQ=255)
                                                                                aaaGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATg                                                                >  1:2678190/1‑148 (MQ=255)
                                                                                 aaGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGc                                                               <  1:2986855/148‑1 (MQ=255)
                                                                                  aGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGcc                                                              >  2:2296905/1‑148 (MQ=255)
                                                                                   gAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCt                                                             <  1:625092/148‑1 (MQ=255)
                                                                                     aTCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTAc                                                           <  2:864081/148‑1 (MQ=255)
                                                                                     aTCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTAc                                                           <  2:702011/148‑1 (MQ=255)
                                                                                     aTCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTAc                                                           >  2:223962/1‑148 (MQ=255)
                                                                                      tCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACt                                                          <  1:157165/148‑1 (MQ=255)
                                                                                      tCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACt                                                          <  2:1576476/148‑1 (MQ=255)
                                                                                        cGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACttt                                                        <  2:3075497/148‑1 (MQ=255)
                                                                                         gCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTta                                                       <  2:437989/148‑1 (MQ=255)
                                                                                         gCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTta                                                       >  2:1178500/1‑148 (MQ=255)
                                                                                          cATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTtat                                                      <  1:1691521/148‑1 (MQ=255)
                                                                                          cATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTtat                                                      >  2:558161/1‑148 (MQ=255)
                                                                                           aTACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTtata                                                     >  1:2273760/1‑148 (MQ=255)
                                                                                           aTACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTtata                                                     <  2:333701/148‑1 (MQ=255)
                                                                                               cAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGc                                                 <  2:1264046/148‑1 (MQ=255)
                                                                                                 ggAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGTCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCAt                                               >  2:1970523/1‑148 (MQ=255)
                                                                                                 ggAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCAt                                               >  2:1278020/1‑148 (MQ=255)
                                                                                                 ggAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCAt                                               >  2:1968121/1‑148 (MQ=255)
                                                                                                     gggCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAagc                                           >  1:1531863/1‑148 (MQ=255)
                                                                                                     gggCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAagc                                           <  2:3309653/148‑1 (MQ=255)
                                                                                                      ggCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAagca                                          <  2:2928237/148‑1 (MQ=255)
                                                                                                           tGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCa                                     >  2:801192/1‑148 (MQ=255)
                                                                                                              gAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAaca                                  >  2:2966054/1‑148 (MQ=255)
                                                                                                               aaaCACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAacac                                 >  2:1629753/1‑148 (MQ=255)
                                                                                                                aaCACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACAcc                                <  1:108157/148‑1 (MQ=255)
                                                                                                                 acacTGCCCTGTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACAccc                               <  2:1674832/148‑1 (MQ=255)
                                                                                                                     tGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTAt                           <  2:3275790/148‑1 (MQ=255)
                                                                                                                      gCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTAtt                          >  1:315735/1‑148 (MQ=255)
                                                                                                                       cccTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTc                         >  2:3424235/1‑148 (MQ=255)
                                                                                                                        ccTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTcc                        >  1:2085933/1‑148 (MQ=255)
                                                                                                                         cTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCa                       >  2:1910364/1‑148 (MQ=255)
                                                                                                                          tttCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAg                      >  1:1227822/1‑148 (MQ=255)
                                                                                                                           ttCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGc                     >  1:3141377/1‑148 (MQ=255)
                                                                                                                              aGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTgg                  <  2:2032839/148‑1 (MQ=255)
                                                                                                                                   gCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCg             >  1:3169072/1‑148 (MQ=255)
                                                                                                                                    ccATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGt            >  1:2662957/1‑148 (MQ=255)
                                                                                                                                     catcatGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTg           >  2:3017730/1‑148 (MQ=255)
                                                                                                                                      atcatGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTgg          >  2:1595315/1‑148 (MQ=255)
                                                                                                                                        catGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCt        >  1:2194171/1‑148 (MQ=255)
                                                                                                                                         atGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTg       <  1:2888932/148‑1 (MQ=255)
                                                                                                                                           gAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGtt     <  1:3266446/148‑1 (MQ=255)
                                                                                                                                           gAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGtt     >  1:1055023/1‑148 (MQ=255)
                                                                                                                                              tGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTc  >  1:1147317/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
TCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTC  >  CP009273/576853‑577142

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