Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,204,409 T→A 7.4% 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.4% 215.1 / 4.8 90pseudogene (6/1134 nt)aaaEpseudogene, phage terminase protein A family, e14 prophage;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base T (37/46);  new base A (3/4);  total (40/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATTTCAGTAATGACGAAGGTGTGTGCTGGCCATCAATTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACGCCAGCCACTGATGTGGGCCATTACCACCGCCGGGTACAACATTGAGGGGCCGTGCTA  >  NC_000913/1204260‑1204558
                                                                                                                                                     |                                                                                                                                                     
gATTTCAGTAATGACGAAGGTGTGTGCTGGCCATCAATTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGt                                                                                                                                                       >  1:1621413/1‑150 (MQ=255)
   ttCAGTACTGACGAAGGTGTGTGCTGGCCATCAATTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTtca                                                                                                                                                    <  2:1488341/150‑1 (MQ=255)
   ttCAGTAATGACGAAGGTGTGTGCTGGCCATCAATTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTtca                                                                                                                                                    >  2:1674177/1‑150 (MQ=255)
    tCAGTAATGACGAAGGTGTGTGCTGGCCATCAATTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTtcat                                                                                                                                                   >  1:872314/1‑150 (MQ=255)
    tCAGTAATGACGAAGGTGTGTGCTGGCCATCAATTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTtcat                                                                                                                                                   <  1:1387197/150‑1 (MQ=255)
       gTAATGACGAAGGTGTGTGCTGGCCATCAATTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCAc                                                                                                                                                >  2:215542/1‑150 (MQ=255)
       gTAATGACGAAGGTGTGTGCTGGCCATCAATTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCAc                                                                                                                                                <  1:71698/150‑1 (MQ=255)
       gTAATGACGAAGGTGTGTGCTGGCCATCAATTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCAc                                                                                                                                                >  1:453492/1‑150 (MQ=255)
         aaTGACGAAGGTGTGTGCTGGCCATCAATTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCAccc                                                                                                                                              >  2:248274/1‑150 (MQ=255)
          aTGACGAAGGTGTGTGCTGGCCATCAATTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAATGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCAcccc                                                                                                                                             <  1:1026950/150‑1 (MQ=255)
           tGACGAAGGTGTGTGCTGGCCATCAATTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGCAGGCTGGTTAACGCGTAAGGCGCGTCGCCCGGGTGATGGTTCATCCCCCCa                                                                                                                                            >  2:678552/1‑150 (MQ=255)
           tGACGAAGGTGTGTGCTGGCCATCAATTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCa                                                                                                                                            >  2:541212/1‑150 (MQ=255)
            gACGAAGGTGTGTGCTGGCCATCAATTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGCACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCAc                                                                                                                                           >  2:95389/1‑150 (MQ=255)
            gACGAAGGTGTGTGCTGGCCATCAATTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCAc                                                                                                                                           >  2:640667/1‑150 (MQ=255)
            gACGAAGGTGTGTGCTGGCCATCAATTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGAACCATCAccctgg                                                                                                                                           <  1:1488499/150‑4 (MQ=255)
            gACGAAGGTGGGTGCTGGCCATCAATTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGAACCATCAccctgg                                                                                                                                           <  1:517437/150‑4 (MQ=255)
              cGAAGGTGAGTGCTCGCCTTCAATTGAAACCATTCCCCGTCAGATTGGCCCGGGGATCAGTACCGTCAGAACCGCAATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCTTCGCCAGGGTGATGGATCATCACCCCACtg                                                                                                                                         <  1:1209953/150‑1 (MQ=255)
                  ggTGTGTGCTGGCCATCAATTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGcc                                                                                                                                     >  2:784867/1‑150 (MQ=255)
                      tgtgCTGGCCATCAATTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCgtgg                                                                                                                                 >  1:858591/1‑150 (MQ=255)
                      tgggCTGGCCATCACTTGACACCATTGCCCGGCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAATCAGAAGGCTGGTTAACGCTTAAGGCGAGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCgtgg                                                                                                                                 <  1:1209450/147‑1 (MQ=255)
                       gtgCTGGCCATCAATTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCgtggt                                                                                                                                >  1:507160/1‑150 (MQ=255)
                         gCTGGCCATCAATTGATACCATTGCCCGGCAGATTGGCGCGGGGATGAGTACCGGCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGtggtgg                                                                                                                              <  2:1620768/150‑1 (MQ=255)
                          cTGGCCATCAATTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCTCGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGa                                                                                                                             <  2:1491117/150‑1 (MQ=255)
                            ggCCATCAATTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATg                                                                                                                           <  1:1694497/150‑1 (MQ=255)
                            ggCCATCAATTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCGGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATg                                                                                                                           <  2:14523/150‑1 (MQ=255)
                             gCCATCAATTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGa                                                                                                                          <  2:1242030/150‑1 (MQ=255)
                               cATCAATTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAGGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAat                                                                                                                        >  1:991959/1‑150 (MQ=255)
                                aTCAATTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAata                                                                                                                       <  1:1674177/150‑1 (MQ=255)
                                aTCAATTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGAAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAata                                                                                                                       >  2:1053778/1‑150 (MQ=255)
                                  cAATTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATc                                                                                                                     >  2:1175358/1‑150 (MQ=255)
                                      tGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGa                                                                                                                 >  1:946919/1‑150 (MQ=255)
                                       gAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAg                                                                                                                <  2:1621413/150‑1 (MQ=255)
                                        aaaCCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGc                                                                                                               >  1:1560606/1‑150 (MQ=255)
                                        aaaCCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGc                                                                                                               <  2:311734/150‑1 (MQ=255)
                                        aaaCCATAGCCCGTCAGATTGGCGGGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAATGCGCGTCGCCAGGGTGATGGATCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGc                                                                                                               <  2:1056938/150‑1 (MQ=255)
                                             aTTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGcc                                                                                                          <  1:1034000/150‑1 (MQ=255)
                                                gCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGCAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCcaca                                                                                                       <  2:1136845/150‑1 (MQ=255)
                                                gCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCcaca                                                                                                       <  1:1106625/150‑1 (MQ=255)
                                                  ccGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGa                                                                                                     >  2:163977/1‑150 (MQ=255)
                                                   cGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCCGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCTCTAGCACGCCACAGAt                                                                                                    >  1:1453142/1‑150 (MQ=255)
                                                        caTTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCt                                                                                               <  1:1300195/149‑1 (MQ=255)
                                                        caTTGGCGCGGGGAAGAGTACCGTCAGAACGGCTATCGCACGGCTGGACGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGGGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCt                                                                                               <  2:215894/149‑1 (MQ=255)
                                                         aTTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCtt                                                                                              >  2:18972/1‑150 (MQ=255)
                                                         aTTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCtt                                                                                              <  1:1343742/150‑1 (MQ=255)
                                                          ttGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCttt                                                                                             <  2:1017625/150‑1 (MQ=255)
                                                                cgGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCTCAGATGCGCTTTAcacca                                                                                       <  2:858591/150‑1 (MQ=255)
                                                                 ggggATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTAcaccac                                                                                      <  2:1468816/150‑1 (MQ=255)
                                                                 ggggATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTAcaccac                                                                                      <  1:54159/150‑1 (MQ=255)
                                                                  gggATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACg                                                                                     <  2:176721/150‑1 (MQ=255)
                                                                  gggATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACg                                                                                     <  2:983246/150‑1 (MQ=255)
                                                                   ggATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGGGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGa                                                                                    <  1:1276401/150‑1 (MQ=255)
                                                                           aCCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTCCACCACGATGCTTAcc                                                                            >  1:126443/1‑150 (MQ=255)
                                                                           aCCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTAcc                                                                            <  1:901431/150‑1 (MQ=255)
                                                                             cGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCgg                                                                          <  1:686246/150‑1 (MQ=255)
                                                                              gTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCAGGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCggg                                                                         <  1:248274/150‑1 (MQ=255)
                                                                              gTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCggg                                                                         <  1:979935/150‑1 (MQ=255)
                                                                              gTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCggg                                                                         >  2:1395946/1‑150 (MQ=255)
                                                                                  gAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATgg                                                                     >  2:496523/1‑150 (MQ=255)
                                                                                       gCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGgcg                                                                <  2:901601/150‑1 (MQ=255)
                                                                                        cTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGgcgc                                                               >  1:1121006/1‑150 (MQ=255)
                                                                                          aTCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGa                                                             >  2:604007/1‑150 (MQ=255)
                                                                                              cACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACgcc                                                         <  2:679753/150‑1 (MQ=255)
                                                                                                      gAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACGCCAGCCACTg                                                 >  2:627167/1‑150 (MQ=255)
                                                                                                         gCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACGCCAGCCACTGAtg                                              <  1:640667/150‑1 (MQ=255)
                                                                                                         gCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACAGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACGCCAGCCACTGAtg                                              <  1:95389/150‑1 (MQ=255)
                                                                                                             caGGCTGGGTTACGGGTTATGCGCGTCGGCATGGTGATGGTTTATTCCCCCACTGTGGCGTGGGGGATGAATATCACGACCACGCCACAGTTGCGGTTTACACCACGATGCTTACCGGTATGGGGGCGCGACGCCAGCCACTGATGTggg                                          <  1:1125320/149‑1 (MQ=255)
                                                                                                               ggCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACGCCAGCCACTGATGTGGGcc                                        <  2:872314/150‑1 (MQ=255)
                                                                                                               ggCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACAGGGATGGGGGCGCGACGCCAGCCACTGATGTGGGcc                                        >  2:252151/1‑150 (MQ=255)
                                                                                                                 cTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACGCCAGCCACTGATGTGGGCCAt                                      >  1:632545/1‑150 (MQ=255)
                                                                                                                  tGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACGCCAGCCACTGATGTGGGCCAtt                                     <  2:1560606/150‑1 (MQ=255)
                                                                                                                   ggTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACTCCAGCCACTGATGTGGGCCATTa                                    <  1:1244559/150‑1 (MQ=255)
                                                                                                                   ggTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACGCCAGCCACTGATGTGGGCCATTa                                    <  2:244390/150‑1 (MQ=255)
                                                                                                                      tAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACGCCAGCCACTGATGTGGGCCATTacca                                 >  1:439580/1‑150 (MQ=255)
                                                                                                                           cgTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACGCCAGCCACTGATGTGGGCCATTACCAccgcc                            >  1:1216867/1‑150 (MQ=255)
                                                                                                                           cgTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACGCCAGCCACTGATGTGGGCCATTACCAccgcc                            >  1:670807/1‑150 (MQ=255)
                                                                                                                            gTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACGCCAGCCACTGATGTGGGCCATTACCAccgccg                           <  1:739280/150‑1 (MQ=255)
                                                                                                                              aaGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACGCCAGCCACTGATGTGGGCCATTACCACCGCCggg                         <  1:163977/150‑1 (MQ=255)
                                                                                                                              aaGGCGCGTCGCAATGGTGATGGTTCATCACCCCACTGTGCCGTGGTGTATGAATATCACGAGAACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACGCCAGCCACTGATGTGGGCCATTACCACCGCCggg                         <  1:376323/150‑1 (MQ=255)
                                                                                                                                 gcgcGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACGCCAGCCACTGATGTGGGCCATTACCACCGCCGGGTac                      <  2:453492/150‑1 (MQ=255)
                                                                                                                                    cgTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACGCCAGCCACTGATGTGGGCCATTACCACCGCCGGGTacaac                   >  1:957377/1‑150 (MQ=255)
                                                                                                                                      acagtggggTGATGAACCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACGCCCGCCACTGATGTGGGCCATTACCACCGCCGGGTACAACAt                 >  1:73655/7‑150 (MQ=255)
                                                                                                                                         ccAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACGCCAGCCACTGATGTGGGCCATTACCACCGCCGGGTACAACATTGa              >  2:1626466/1‑150 (MQ=255)
                                                                                                                                         ccAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACGCCAGCCACTGATGTGGGCCATTACCACCGCCGGGTACAACATTGa              >  1:1339784/1‑150 (MQ=255)
                                                                                                                                          cAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACGCCAGCCACTGATGTGGGCCATTACCACCGCCGGGTACAACATTGAg             >  1:1210881/1‑150 (MQ=255)
                                                                                                                                           ggggTGATGAACCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACGCCAGCCACTGATGTGGGCCATTACCACCGCCGGGTACAACATTGAgg            >  2:530915/2‑150 (MQ=255)
                                                                                                                                           ggggTGATGAACCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACGCCAGCCACTGATGTGGGCCATTACCACCGCCGGGTACAACATTGAgg            >  1:756040/2‑150 (MQ=255)
                                                                                                                                             ggTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACGCCAGCCACTGATGTGGGCCATTACCACCGCCGGGTACAACATTGAgggg          <  2:946919/150‑1 (MQ=255)
                                                                                                                                             ggTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGAAGCTTACCGGGATGGGGGCGCGACGCCAGCCACTGATGTGGGCCATTACCACCGCCGGGTACAACATTGAgggg          <  2:1602456/150‑1 (MQ=255)
                                                                                                                                                  tGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACGCCAGCCACTGATGTGGGCCATTACCACCGCCGGGTACAACATTGAGGGGCCGTg     >  1:320603/1‑150 (MQ=255)
                                                                                                                                                  tGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACGCCAGCCACTGATGTGGGCCATTACCACCGCCGGGTACAACATTGAGGGGCCGTg     <  2:439580/150‑1 (MQ=255)
                                                                                                                                                   ggTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACGCCAGCCACTGATGTGGGCCATTACCACCGCCGGGTACAACATTGAGGGGCCGTGc    >  1:950425/1‑150 (MQ=255)
                                                                                                                                                     ttCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACGCCAGCCACTGATGTGGGCCATTACCACCGCCGGGTACAACATTGAGGGGCCGTGCTa  >  1:1404191/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GATTTCAGTAATGACGAAGGTGTGTGCTGGCCATCAATTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACGCCAGCCACTGATGTGGGCCATTACCACCGCCGGGTACAACATTGAGGGGCCGTGCTA  >  NC_000913/1204260‑1204558

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