Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,204,409 T→A 7.6% pseudogene (6/1134 nt) aaaE → pseudogene, phage terminase protein A family, e14 prophage;Phage or Prophage Related

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,204,4090TA 7.6% 302.6 / 12.5 116pseudogene (6/1134 nt)aaaEpseudogene, phage terminase protein A family, e14 prophage;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base T (53/54);  new base A (7/2);  total (60/56)
Fisher's exact test for biased strand distribution p-value = 1.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.50e-01

Only 100 of 117 total aligned reads displayed.

TTCAGTAATGACGAAGGTGTGTGCTGGCCATCAATTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACGCCAGCCACTGATGTGGGCCATTACCACCGCCGGGTACAACATTGAGGGGCCGTGCTA  >  NC_000913/1204263‑1204558
                                                                                                                                                  |                                                                                                                                                     
ttCAGTAATGACGAAGGTGTGTGCTGGCCATCACTTGAAACCATTCCCCGGCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCGATCGCACGGCGGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTtca                                                                                                                                                    <  1:153128/150‑1 (MQ=255)
ttCAGTAATGACGAAGGTGTGTGCTGGCCATCAATTGAAACCATTGCCCGGCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTtca                                                                                                                                                    <  2:90275/150‑1 (MQ=255)
 tCAGTAATGACGAAGGTGTGTGCTGGCCATCAATAGAAACCATTGCCCGTCCGATTGGCGCTGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTtcat                                                                                                                                                   <  2:108416/150‑1 (MQ=255)
  cAGTAATGACGAAGGTGTGTGCTGGCCATCAATTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTtcatc                                                                                                                                                  >  1:629442/1‑150 (MQ=255)
    gTAATGACGAAGGTGTGTGCTGGCCATCAATTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCAc                                                                                                                                                >  1:1071661/1‑150 (MQ=255)
        tGACGAAGGTGTGTGCTGGCCATCAATTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCa                                                                                                                                            >  1:1524030/1‑150 (MQ=255)
        tGACGAAGGTGTGTGCTGGCCATCAATTGAAACCATTGACCGTCAGATTGGCGCGGGGATGAGTACAGTCAGAACGGATATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCa                                                                                                                                            >  1:1519834/1‑150 (MQ=255)
        tGAAGAAGGTGTGTGCTGGCCATCATTTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGAACCATCAccctg                                                                                                                                            <  2:887094/150‑3 (MQ=255)
           cGAAGGTGTGTGGTGGGCCTCAATTGAAACCATTGCCCTTCAGATTGGCGCGGGGATGAGTACCGTCTGAACGTCTATAGAACTGCAGGAGGCAGCCGGCTGGTTACGGCGTAAGGCGCGTCGCCAGGGTGAGGGTTCATCGCCCCAcgt                                                                                                                                         >  2:300303/1‑148 (MQ=255)
              aGGTGGGTGCGGGCCATCAAGTGTAACCATTGCCCGGCTGATTGGCGCGGGGATGAGTACCGTCAGATCGGCTATCGCTCGGCTGGAAGCAGAAGGCTGGATAACGCGTAAGGCGCGGGGCCAGGGTGATGGTTCTTCACCCCACTGTGc                                                                                                                                      <  1:331929/150‑1 (MQ=255)
                   tgtgCTGGCCATCAATTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCgtgg                                                                                                                                 <  2:1426140/150‑1 (MQ=255)
                   tgtgCTGGCCATCAATTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCgtgg                                                                                                                                 <  2:1812935/150‑1 (MQ=255)
                   tgtgCTGGCCATCAATTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCgtgg                                                                                                                                 <  1:425482/150‑1 (MQ=255)
                     tgCTGGCCATCAATTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTTATTGTTCATCACCCCACTGTGCCgtggtg                                                                                                                               <  1:1227984/150‑1 (MQ=255)
                      gCTGGCCATCAATTGAAACCATTGCCCGTCAGTTTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGtggtgg                                                                                                                              <  1:659634/150‑1 (MQ=255)
                       cTGGCCATCAATTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGa                                                                                                                             <  1:468602/150‑1 (MQ=255)
                       cTGGCCATCAATTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCATAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGa                                                                                                                             <  2:1426685/150‑1 (MQ=255)
                        tGGCCATCAATTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGAt                                                                                                                            <  1:1025663/150‑1 (MQ=255)
                            cATCAATTGAAACCATTGCCCGTCAGCTTGGCGCGGGGATGAGTCCCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAat                                                                                                                        >  2:177164/1‑150 (MQ=255)
                            cATCAATTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAat                                                                                                                        <  2:554873/150‑1 (MQ=255)
                            cATCAATTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGGAAGGCGCGTCGCCAGGGGGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAat                                                                                                                        <  2:706644/150‑1 (MQ=255)
                             aTCAATTCACACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAata                                                                                                                       <  2:1500510/150‑1 (MQ=255)
                                aaTTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCa                                                                                                                    >  1:127790/1‑150 (MQ=255)
                                aaTTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCa                                                                                                                    >  2:1725421/1‑150 (MQ=255)
                                 aTTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCCCTGTGACGTGGTGGATGAATATCAc                                                                                                                   >  1:471938/1‑150 (MQ=255)
                                 aTTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCAc                                                                                                                   >  1:688676/1‑150 (MQ=255)
                                 aTTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCAc                                                                                                                   <  2:1334720/150‑1 (MQ=255)
                                  ttGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAAAGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACg                                                                                                                  <  2:1688921/150‑1 (MQ=255)
                                    gAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAg                                                                                                                <  2:1737864/150‑1 (MQ=255)
                                      aaCCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGCTGGATGAATATCACGAGCa                                                                                                              >  1:1192600/1‑150 (MQ=255)
                                       aCCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCAc                                                                                                             <  1:796740/150‑1 (MQ=255)
                                         cATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGc                                                                                                           <  1:1009415/150‑1 (MQ=255)
                                          aTTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGcc                                                                                                          >  1:1347439/1‑150 (MQ=255)
                                          aTTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGcc                                                                                                          <  1:1145536/150‑1 (MQ=255)
                                             gCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCcaca                                                                                                       <  1:404749/150‑1 (MQ=255)
                                               ccGTCTGATTGGCGCGGGGATGAGTAACGTCAGAACGGCTATCGCACGGCTGTAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGa                                                                                                     <  2:810175/150‑1 (MQ=255)
                                               ccGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCCCCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGa                                                                                                     >  2:424671/1‑150 (MQ=255)
                                               ccGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGa                                                                                                     <  2:450571/150‑1 (MQ=255)
                                                cGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGAt                                                                                                    >  1:1662103/1‑150 (MQ=255)
                                                 gTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATg                                                                                                   >  2:479287/1‑150 (MQ=255)
                                                  tCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATgc                                                                                                  <  2:860760/150‑1 (MQ=255)
                                                     gATTGGCGCGGGGATGAGTACCGTCAGACCGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGGGGATGAATATCACGCGCACGCCACAGATGCGCt                                                                                               >  1:1362119/1‑150 (MQ=255)
                                                     gATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTTACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGCGCACGCCACAGCTGCGCt                                                                                               >  1:1354698/1‑150 (MQ=255)
                                                     gATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCt                                                                                               <  1:1545146/150‑1 (MQ=255)
                                                     gATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCt                                                                                               <  1:1039630/150‑1 (MQ=255)
                                                     gATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCt                                                                                               >  2:1653400/1‑150 (MQ=255)
                                                     gATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCt                                                                                               <  2:208271/150‑1 (MQ=255)
                                                     gATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCt                                                                                               >  2:1335635/1‑150 (MQ=255)
                                                        tGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTa                                                                                            <  2:901367/150‑1 (MQ=255)
                                                          gcgcGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCTCTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTaca                                                                                          >  2:118371/1‑150 (MQ=255)
                                                          gcgcGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTaca                                                                                          >  1:753166/1‑150 (MQ=255)
                                                              ggggAGGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGGGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTAcaccac                                                                                      <  2:1512758/150‑1 (MQ=255)
                                                               gggATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACg                                                                                     >  1:1433650/1‑150 (MQ=255)
                                                                   tGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGc                                                                                 <  2:379736/150‑1 (MQ=255)
                                                                   tGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGc                                                                                 <  2:1449776/150‑1 (MQ=255)
                                                                      gTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTa                                                                              >  2:1555502/1‑150 (MQ=255)
                                                                         ccGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTTCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCg                                                                           >  2:218840/1‑150 (MQ=255)
                                                                           gTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCggg                                                                         >  1:1024366/1‑150 (MQ=255)
                                                                               gAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGCTGCTTACCGGGATgg                                                                     >  2:1227890/1‑150 (MQ=255)
                                                                                aaCGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATggg                                                                    >  2:1677264/1‑150 (MQ=255)
                                                                                   gggtATCGCACGGCTGGAAGCAGACGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGgc                                                                 <  2:78384/147‑1 (MQ=255)
                                                                                   ggCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCCGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACCCCACGATGCTTACCGGGATGGGGgc                                                                 >  1:625496/1‑150 (MQ=255)
                                                                                   ggCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGgc                                                                 >  2:1511908/1‑150 (MQ=255)
                                                                                      tATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGcgcg                                                              >  1:596443/1‑150 (MQ=255)
                                                                                        tCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGcgcgcc                                                            >  1:1177224/1‑148 (MQ=255)
                                                                                        tCGCACGGAGGGCAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGATCATCACCCCACTGTGCCGGGGTGGAGGAATATCCCGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGCTGGGGGCGCGAc                                                            <  2:76313/150‑1 (MQ=255)
                                                                                                cTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCCGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATACTTACCGGGATGGGGGCGCGACgccagccc                                                    >  1:864518/1‑149 (MQ=255)
                                                                                                  ggAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACGCCAGCCACt                                                  >  1:1490309/1‑150 (MQ=255)
                                                                                                         gAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACGCCAGCCACTGATGTgg                                           >  2:1360215/1‑150 (MQ=255)
                                                                                                            ggCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACGCCAGCCACTGATGTGGGcc                                        >  2:366571/1‑150 (MQ=255)
                                                                                                             gCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACAGGGATGGGGGCGCGACGCCAGCCACTGATGTGGGCCa                                       >  1:1386419/1‑150 (MQ=255)
                                                                                                              cTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACGCCAGCCACTGATGTGGGCCAt                                      >  2:561045/1‑150 (MQ=255)
                                                                                                               tGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACGCCAGCCACTGATGTGGGCCAtt                                     >  1:183974/1‑150 (MQ=255)
                                                                                                                ggctAACGCGCAAGGCGCGTCGCCAGGGGGATGGTTCATCACCCCACTGAGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACGCCAGCCACTGATGTGGGCCATTa                                    <  1:1743712/147‑1 (MQ=255)
                                                                                                                 gTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACGCCAGCCACTGATGTGGGCCATTac                                   <  1:215081/150‑1 (MQ=255)
                                                                                                                    aaCGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACCGATGCGCTTTACACCAGGATGCTTACCGGGATGGGGGCGCGACGCCAGCCACTGATGTGGGCGATTaccac                                >  1:506259/1‑150 (MQ=255)
                                                                                                                    aaCGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACGCCAGCCACTGATGTGGGCCATTaccac                                <  2:635304/150‑1 (MQ=255)
                                                                                                                      cgcgTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACCCCCCGATGCTTACCGGGATGGGGGCGCGCCGCCAGCCACTGATGTGGGCCATTaccaccc                              >  1:1114397/1‑149 (MQ=255)
                                                                                                                          tAAGGCGCGTCGTCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACGCCAGCCACTGATGTGGGCCATTACCCCCGCCgg                          >  2:1395964/1‑150 (MQ=255)
                                                                                                                          tAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGCTGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACGCCAGCCACTGATGTGGGCCATTACCACCGCCgg                          >  1:1824117/1‑150 (MQ=255)
                                                                                                                          tAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACGCCAGCCACTGATGTGGGCCATTACCACCGCCgg                          >  1:36246/1‑150 (MQ=255)
                                                                                                                           aaGGCGCGTCGCCAGGGGGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACGCCAGCCACTGATGTGGGCCATTACCACCGCCggg                         <  1:987785/150‑1 (MQ=255)
                                                                                                                            aGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACGCCAGCCACTGATGTGGGCCATTACCACCGCCGGGt                        <  2:564217/150‑1 (MQ=255)
                                                                                                                              gcgcGTCGCCCGGGGGATGGTTGATCACCCCACTGGGCCGGGGTGGAAGAATATCACGAGCACGCCAAAGAAGCGCGTTACACCACGATGCTTACCGGGATGGGGGCGCGACGCCAGCCACTGATGTGGGCCATTACCACCGCCGGGTac                      <  2:915160/150‑1 (MQ=255)
                                                                                                                                 cgTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACGCCAGCCACTGATGTGGGCCATTACCACCGCCGGGTacaac                   >  2:1670109/1‑150 (MQ=255)
                                                                                                                                  gTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACGCCAGCCACTGATGTGGGCCATTACCACCGCCGGGTacaaca                  <  2:819999/150‑1 (MQ=255)
                                                                                                                                    cagtggggTGATGAACCATCACCCCACTGTGCCGTGGTGGATGAATATGACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGAGGGGGGCGCGACGCCCGCCACTGATGTGGGCCATTACCCCCGCCGGGTACAACAtt                >  1:927298/6‑150 (MQ=255)
                                                                                                                                    cagtggggTGATGAACCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGAAGGGGGCGCGACGCCAGCCACTGATGTGGGCCATTACCACCGCCGGGTACAACAtt                >  2:949130/6‑150 (MQ=255)
                                                                                                                                     gCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACGCCAGCCACTGATGTGGGCCATTACCACCGCCGGGTACAACATTg               >  1:1049085/1‑150 (MQ=255)
                                                                                                                                      gtggggTGATGAACCATCACCCTACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACGCCAGCCACTGATGTGGGCCATTACCACCGCCGGGTACAACATTgc              >  1:546586/4‑149 (MQ=255)
                                                                                                                                      gtggggTGATGAACCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACGCCAGCCACTGATGTGGGCCATTACCACCGCCGGGTACAACATTGa              >  2:558684/4‑150 (MQ=255)
                                                                                                                                      ccAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACGCCAGCCACTGATGTGGGCCATTACCACCGCCGGGTACAACATTGa              <  1:1300601/150‑1 (MQ=255)
                                                                                                                                        ggggTGATGAACCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGCCGCCAGCCACTGATGTGGGCCATTACCACCGCCGGGTACAACATTGAgg            >  1:878461/2‑150 (MQ=255)
                                                                                                                                        ggggTGATGAACCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACGCCAGCCACTGATGTGGGCCATTACCACCGCCGGGTACAACATTGAgg            >  1:1734017/2‑150 (MQ=255)
                                                                                                                                        aGGGTGAGGGTTCATCACCCCCCTGTGCCGTGGTGGATGAAAATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACGCCAGCCACTGATGTGGGCCATTACCACCGCCGGGTACAACATTGAgg            <  2:1192600/150‑1 (MQ=255)
                                                                                                                                          ggTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACGCCAGCCACTGATGTGGGCCATTACCACCGCCGGGTACAACATTGAgggg          >  2:1672611/1‑150 (MQ=255)
                                                                                                                                              aTGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACGCCAGCCACTGATGTGGGCCATTACCACCGCCGGGTACAACATTGAGGGGCCGt      <  2:1524030/150‑1 (MQ=255)
                                                                                                                                               tGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACGCCAGCCACTGATGTGGGCCATTACCACCGCCGGGTACAACATTGAGGGGCCGTg     <  2:1490309/150‑1 (MQ=255)
                                                                                                                                                 gTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACGCCAGCCACTGATGTGGGCCATTACCACCGCCGGGTACAACATTGAGGGGCCGTGCt   <  2:45609/150‑1 (MQ=255)
                                                                                                                                                  ttCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCCCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACGCCAGCCACTGATGTGGGCCATTACCACCGCCGGGTACAACATTGAGGGGCCGTGCTa  >  1:955536/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
TTCAGTAATGACGAAGGTGTGTGCTGGCCATCAATTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACGCCAGCCACTGATGTGGGCCATTACCACCGCCGGGTACAACATTGAGGGGCCGTGCTA  >  NC_000913/1204263‑1204558

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