Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,204,408 G→A 5.5% pseudogene (5/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,4080GA 5.5% 289.7 / 4.0 90pseudogene (5/1134 nt)aaaEpseudogene, phage terminase protein A family, e14 prophage;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base G (37/48);  new base A (3/2);  total (40/50)
Fisher's exact test for biased strand distribution p-value = 6.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.28e-01

TGATTTCAGTAATGACGAAGGTGTGTGCTGGCCATCAATTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACGCCAGCCACTGATGTGGGCCATTACCACCGCCGGGTACAACATTGAGGGGCCGTGC  >  NC_000913/1204259‑1204556
                                                                                                                                                     |                                                                                                                                                    
tGATTTCAGTAATGACGAAGGTGTGTGCTGGCCATCAATTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATgg                                                                                                                                                      >  2:1034000/1‑150 (MQ=255)
 gATTTCAGTAATGACGAAGGTGTGTGCTGGCCATCAATTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGt                                                                                                                                                     >  1:1621413/1‑150 (MQ=255)
    ttCAGTACTGACGAAGGTGTGTGCTGGCCATCAATTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTtca                                                                                                                                                  <  2:1488341/150‑1 (MQ=255)
    ttCAGTAATGACGAAGGTGTGTGCTGGCCATCAATTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTtca                                                                                                                                                  >  2:1674177/1‑150 (MQ=255)
     tCAGTAATGACGAAGGTGTGTGCTGGCCATCAATTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTtcat                                                                                                                                                 <  1:1387197/150‑1 (MQ=255)
     tCAGTAATGACGAAGGTGTGTGCTGGCCATCAATTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTtcat                                                                                                                                                 >  1:872314/1‑150 (MQ=255)
        gTAATGACGAAGGTGTGTGCTGGCCATCAATTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCAc                                                                                                                                              >  1:453492/1‑150 (MQ=255)
        gTAATGACGAAGGTGTGTGCTGGCCATCAATTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCAc                                                                                                                                              <  1:71698/150‑1 (MQ=255)
        gTAATGACGAAGGTGTGTGCTGGCCATCAATTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCAc                                                                                                                                              >  2:215542/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)
                                                                                                                                                     |                                                                                                                                                    
TGATTTCAGTAATGACGAAGGTGTGTGCTGGCCATCAATTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACGCCAGCCACTGATGTGGGCCATTACCACCGCCGGGTACAACATTGAGGGGCCGTGC  >  NC_000913/1204259‑1204556

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: