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,9980GT97.2% 824.2 / 11.0 247pseudogene (189/309 nt)aaaDpseudogene; DLP12 prophage; tail fiber assembly protein family;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base G (5/2);  new base T (124/116);  total (129/118)
Fisher's exact test for biased strand distribution p-value = 4.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 247 total aligned reads displayed.

TGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCA  >  CP009273/576851‑577143
                                                                                                                                                   |                                                                                                                                                 
 gTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATg                                                                                                                                                >  2:1162163/1‑150 (MQ=255)
 gTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATg                                                                                                                                                >  1:1734862/1‑150 (MQ=255)
       gcCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGc                                                                                                                                          >  1:2606208/1‑150 (MQ=255)
        ccGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCg                                                                                                                                         <  1:2941856/150‑1 (MQ=255)
         cGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGa                                                                                                                                        >  2:76351/1‑150 (MQ=255)
          gggCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGAc                                                                                                                                       >  1:2450145/1‑150 (MQ=255)
           ggCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACt                                                                                                                                      <  2:1680219/150‑1 (MQ=255)
           ggCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACt                                                                                                                                      <  2:3584654/150‑1 (MQ=255)
           ggCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACt                                                                                                                                      <  1:3997942/150‑1 (MQ=255)
            gCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTa                                                                                                                                     <  2:1545562/150‑1 (MQ=255)
                 gCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATcc                                                                                                                                >  1:4036675/1‑150 (MQ=255)
                 gCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATcc                                                                                                                                <  2:676382/150‑1 (MQ=255)
                  cTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCg                                                                                                                               >  2:2987002/1‑150 (MQ=255)
                   tCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGt                                                                                                                              >  1:1877708/1‑150 (MQ=255)
                       tGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAgg                                                                                                                          >  1:3192399/1‑150 (MQ=255)
                            aGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAAt                                                                                                                     >  1:155126/1‑150 (MQ=255)
                             ggCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAAtg                                                                                                                    >  2:830239/1‑150 (MQ=255)
                                cGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATgtgg                                                                                                                 >  1:3449690/1‑150 (MQ=255)
                                   aGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGa                                                                                                              >  2:3820112/1‑150 (MQ=255)
                                      cacaGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGTAGGGCTCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGt                                                                                                           <  2:1237242/150‑1 (MQ=255)
                                        caGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCt                                                                                                         <  2:1801119/150‑1 (MQ=255)
                                         aGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTg                                                                                                        >  2:655128/1‑150 (MQ=255)
                                         aGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTg                                                                                                        >  2:2266585/1‑150 (MQ=255)
                                          gACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGc                                                                                                       >  1:850344/1‑150 (MQ=255)
                                           aCCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGcc                                                                                                      <  1:2870498/150‑1 (MQ=255)
                                            ccgccgTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGccc                                                                                                     <  2:1075265/150‑1 (MQ=255)
                                             cgccgTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCg                                                                                                    <  2:2321249/150‑1 (MQ=255)
                                             cgccgTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCg                                                                                                    <  2:3996248/150‑1 (MQ=255)
                                              gccgTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCgt                                                                                                   <  1:1421953/150‑1 (MQ=255)
                                                cgTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTCCATCCGTGAGGTGAAGGTGGTGAAGTCTGCCCgtgt                                                                                                 >  2:1147743/1‑150 (MQ=255)
                                                 gTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTc                                                                                                >  2:2529802/1‑150 (MQ=255)
                                                   tGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCgg                                                                                              >  1:3560299/1‑150 (MQ=255)
                                                    gAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGt                                                                                             <  2:2650751/150‑1 (MQ=255)
                                                      aTGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTa                                                                                           <  2:3236967/150‑1 (MQ=255)
                                                        gggCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTAtt                                                                                         >  1:2016378/1‑150 (MQ=255)
                                                              aTGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCaaaa                                                                                   <  2:155129/150‑1 (MQ=255)
                                                                  tAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAAtgct                                                                               <  2:2411725/150‑1 (MQ=255)
                                                                  tAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAAtgct                                                                               >  1:1342262/1‑150 (MQ=255)
                                                                  tAATTACTATCTCCAGATAGAATCCGCATACCAGGTAGGGCGCTGGGTAACACTGCCCTTTCTGCGGGCCATCATGAATGCTATGGGCTGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTACAAAAtgct                                                                               <  1:2514052/150‑1 (MQ=255)
                                                                     ttACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATgctgct                                                                            >  2:3344905/1‑150 (MQ=255)
                                                                     ttACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATgctgct                                                                            >  2:198827/1‑150 (MQ=255)
                                                                      tACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGctgctg                                                                           >  2:509860/1‑150 (MQ=255)
                                                                         tATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGgt                                                                        <  2:3132135/150‑1 (MQ=255)
                                                                         tATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGgt                                                                        <  1:781256/150‑1 (MQ=255)
                                                                          aTCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGgtg                                                                       <  1:2341002/150‑1 (MQ=255)
                                                                           tctcCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGgtgt                                                                      >  1:2574176/1‑150 (MQ=255)
                                                                            ctcCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGtt                                                                     >  1:1415042/1‑150 (MQ=255)
                                                                             tcCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGttt                                                                    <  1:4051011/150‑1 (MQ=255)
                                                                              cccGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTa                                                                   <  1:209014/150‑1 (MQ=255)
                                                                               ccGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTAt                                                                  >  2:3506359/1‑150 (MQ=255)
                                                                               ccGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTAt                                                                  >  1:887318/1‑150 (MQ=255)
                                                                                cGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGGGAGGTGAATGGGGTGAAGTCTGCCCGTGTCGGGTATTCCAAAATGCTGCTGGGTGTTTATg                                                                 >  2:277150/1‑150 (MQ=255)
                                                                                 gAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGc                                                                <  1:58352/150‑1 (MQ=255)
                                                                                 gAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGc                                                                >  2:612472/1‑150 (MQ=255)
                                                                                       aTCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACtt                                                          <  1:760338/150‑1 (MQ=255)
                                                                                        tCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTACAAAATGCTGCTGGGTGTTTATGACTACttt                                                         <  1:3490056/150‑1 (MQ=255)
                                                                                           gCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTtata                                                      <  1:635995/150‑1 (MQ=255)
                                                                                            cATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATag                                                     <  1:1359244/150‑1 (MQ=255)
                                                                                            cATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATag                                                     <  1:611583/150‑1 (MQ=255)
                                                                                             aTACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATaga                                                    <  2:3932202/150‑1 (MQ=255)
                                                                                               aCCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGc                                                  >  2:2546552/1‑150 (MQ=255)
                                                                                               aCCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGc                                                  <  1:1082625/150‑1 (MQ=255)
                                                                                                ccAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCa                                                 >  1:4142490/1‑150 (MQ=255)
                                                                                                 cAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCAt                                                <  1:2917144/150‑1 (MQ=255)
                                                                                                 cAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCAt                                                <  2:2942570/150‑1 (MQ=255)
                                                                                                  aGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATa                                               >  2:2033101/1‑150 (MQ=255)
                                                                                                  aGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATa                                               >  1:1649860/1‑150 (MQ=255)
                                                                                                     aaGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAagc                                            >  1:563806/1‑150 (MQ=255)
                                                                                                     aaGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAagc                                            >  2:2329790/1‑150 (MQ=255)
                                                                                                      aGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAagca                                           <  1:4173248/150‑1 (MQ=255)
                                                                                                      aGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAagca                                           <  2:2967476/150‑1 (MQ=255)
                                                                                                            cTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCaa                                     >  1:1936834/1‑150 (MQ=255)
                                                                                                               ggAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAacac                                  >  1:2801682/1‑150 (MQ=255)
                                                                                                               ggAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAacac                                  <  2:2857346/150‑1 (MQ=255)
                                                                                                               ggAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAacac                                  >  1:1224705/1‑150 (MQ=255)
                                                                                                                 aaaCACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACAccc                                <  1:2948044/150‑1 (MQ=255)
                                                                                                                  aaCACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCt                               >  2:3412703/1‑150 (MQ=255)
                                                                                                                   acacTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCtt                              >  1:1567089/1‑150 (MQ=255)
                                                                                                                      cTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTAtt                           >  2:2788587/1‑150 (MQ=255)
                                                                                                                       tGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTc                          >  2:2615699/1‑150 (MQ=255)
                                                                                                                         cccTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCa                        <  2:2136337/150‑1 (MQ=255)
                                                                                                                            tttCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCt                     >  1:3380935/1‑150 (MQ=255)
                                                                                                                            tttCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCt                     <  1:514722/150‑1 (MQ=255)
                                                                                                                             ttCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTg                    <  2:1290861/150‑1 (MQ=255)
                                                                                                                              tCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTgg                   <  1:2535088/150‑1 (MQ=255)
                                                                                                                               cAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGc                  >  1:3988078/1‑150 (MQ=255)
                                                                                                                                aGCGGGCTTTCTTGAATGCTATGGGCATCGACTACATCCGTGAGGTGAATGTGGTGTAGTCTGCAAGTGTGGGTTATTCCAATATGCTGCTGGGTGTTTATGCCTTCTTTATAGAGCATAATCAGCTCATCACCCTTATTACATCTGGCt                 <  2:205748/150‑1 (MQ=255)
                                                                                                                                aGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCt                 <  1:2208733/150‑1 (MQ=255)
                                                                                                                                  cGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTc               <  1:2273076/150‑1 (MQ=255)
                                                                                                                                    ggCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGt             >  1:613118/1‑150 (MQ=255)
                                                                                                                                     gCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTg            <  2:2953468/150‑1 (MQ=255)
                                                                                                                                     gCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTg            <  1:3353274/150‑1 (MQ=255)
                                                                                                                                       catcatGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGc          >  1:1009996/1‑150 (MQ=255)
                                                                                                                                        atcatGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCt         <  1:2284202/150‑1 (MQ=255)
                                                                                                                                        atcatGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCt         >  1:3843359/1‑150 (MQ=255)
                                                                                                                                          catGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGt       >  2:3539442/1‑150 (MQ=255)
                                                                                                                                          catGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGt       >  2:2819451/1‑150 (MQ=255)
                                                                                                                                           atGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGtt      >  2:2068831/1‑150 (MQ=255)
                                                                                                                                               aTGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCa  <  2:3019667/150‑1 (MQ=255)
                                                                                                                                               aTGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCa  >  1:3027881/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCA  >  CP009273/576851‑577143

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: