Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,624,195 C→G 5.9% D384E (GAC→GAG)  ppk → polyphosphate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,624,1950CG 5.9% 287.7 / 7.1 85D384E (GAC→GAG) ppkpolyphosphate kinase
Reads supporting (aligned to +/- strand):  ref base C (45/35);  new base G (3/2);  total (48/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.80e-01

CATTCGCGAACGCGATGTGTTGCTCTATTATCCTTATCACACCTTTGAGCATGTGCTGGAACTGCTGCGTCAGGCTTCGTTCGACCCGAGCGTACTGGCGATTAAAATTAACATTTACCGCGTGGCGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTCGACGAAGAAGCCAACATTCACTGGGCGAAGCGCCTGACCGAAGCAGGCGTGCACGTTATCTTCTCTGCGCCGGGGCTGAAAAT  >  NC_000913/2624048‑2624344
                                                                                                                                                   |                                                                                                                                                     
cATTCGCGAACGCGATGTGTTGCTCTATTATCCTTATCACACCTTTGAGCATGTGCTGGAACTGCTGCGTCAGGCTTCGTTCGACCCGAGCGTACTGGCGATTAAAATTAACATTTACCGCGTGGCGAAAGATTCACGCATCATCGACTc                                                                                                                                                     >  1:1373093/1‑150 (MQ=255)
 aTTCGCGGACGCGATGTGTTGCTCTATTATCCTTATCACACCTTTGAGCATGTGCTGGAACTGCTGCGTCAGGCTTCGTTCGACCCGAGCGTACTGGCGATTAAAATTAACATTTACCGCGTGGCGAAAGATTCACGCATCATCGACTCg                                                                                                                                                    >  2:675308/1‑150 (MQ=255)
   tCGCGAACGCGATGTGTTGCTCTATTATCCTTATCACACCTTTGAGCATGTGCTGGAACTGCTGCGTCAGGCTTCGTTCGACCCGAGCGTACTGGCGATTAAAATTAACATTTACCGCGTGGCGAAAGATTCACGCATCATCGACTCgat                                                                                                                                                  <  2:1373093/150‑1 (MQ=255)
    cgcgAACGGGATGTGTTGCTCTATTATCCTTATCACACCTTTGAGCATGTGCTGGAACTGCTGCGTCAGGCTTCGTTCGACCCGAGCGTACTGGCGATTAAAATTAACATTTACCGCGTGGCGAAAGATTCACGCATCATCGACTCgatg                                                                                                                                                 <  2:1034572/150‑1 (MQ=255)
    cgcgAACGCGATGTGTTGCTCTATTATCCTTATCACACCTTTGAGCATGTGCTGGAACTGCTGCGTCAGGCTTCGTTCGACCCGAGCGTACTGGCGATTAAAATTAACATTTACCGCGTGGCGAAAGATTCACGCATCATCGACTCgatg                                                                                                                                                 >  2:719910/1‑150 (MQ=255)
      cgAACGCGATGTGTTGCTCTATTATCCTTATCACACCTTTGAGCATGTGCTGGAACTGCTGCGTCAGGCTTCGTTCGACCCGAGCGTACTGGCGATTAAAATTAACATTTACCGCGTGGCGAAAGATTCACGCATCATCGACTCgatgat                                                                                                                                               <  1:679426/150‑1 (MQ=255)
         aCGCGATGTGTTGCTCTATTATCCTTATCACACCTTTGAGCATGTGCTGGAACTGCTGCGTCAGGCTTCGTTCGACCCGAGCGTACTGGCGATTAAAATTAACATTTACCGCGTGGCGAAAGATTCACGCATCATCGAGTCgatgatgcg                                                                                                                                            <  1:408701/150‑4 (MQ=255)
         aCGCGATGTGTTGCTCTATTATCCTTATCACACCTTTGAGCATGTGCTGGAACTGCTGCGTCAGGCTTCGTTCGACCCGAGCGTACTGGCGATTAAAATTAACATTTACCGCGTGGCGAAAGATTCACGCATCATCGAGTCgatgatgcg                                                                                                                                            <  2:13477/150‑4 (MQ=255)
            cgATGTGTTGCTCTATTATCCTTATCACACCTTTGAGCATGTGCTGGAACTGCTGCGTCAGGCTTCGTTCGACCCGAGCGTACTGGCGATTAAAATTAACATTTACCGCGTGGCGAAAGATTCACGCATCATCGACTCGATGATCCAcgc                                                                                                                                         <  1:332291/150‑1 (MQ=255)
            cgATGTGTTGCTCTATTATCCTTATCACACCTTTGAGCATGTGCTGGAACTGCTGCGTCAGGCTTCGTTCGACCCGAGCGTACTGGCGATTAAAATTAACATTTACCGCGTGGCGAAAGATTCACGCATCATCGACTCGATGATCCAcgc                                                                                                                                         >  1:908649/1‑150 (MQ=255)
             gATGTGTTGCTCTATTATCCTTATCACACCTTTGAGCATGTGCTGGAACTGCTGCGTCAGGCTTCGTTCGACCCGAGCGTACTGGCGATTAAAATTAACATTTACCGCGTGGCGAAAGATTCACGCATCATCGACTCGATGATCCAcgcc                                                                                                                                        <  1:719910/150‑1 (MQ=255)
               tgtgTTGCTCTATTATCCTTATCACACCTTTGAGCATGTGCTGGAACTGCTGCGTCAGGCTTCGTTCGACCCGAGCGTACTGGCGATTAAAATTAACATTTACCGCGTGGCGAAAGATTCACGCATCATCGACTCGATGATCCAcgccgc                                                                                                                                      <  2:1115783/150‑1 (MQ=255)
                 tgtTGCTCTATTATCCTTATCACACCTTTGAGCATGTGCTGGAACTGCTGCGTCAGGCTTCGTTCGACCCGAGCGTACTGGCGATTAAAATTAACATTTACCGCGTGGCGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGcac                                                                                                                                    >  1:1061390/1‑150 (MQ=255)
                      ctctATTATCCTTATCACACCTTTGAGCATGTGCTGGAACTGCTGCGTCAGGCTTCGTTCGACCCGAGCGTACTGGCGATTAAAATTAACATTTACCGCGTGGCGAAAGATTCACGCATCATCGCCTCGATGATCCACGCCGCACATAAc                                                                                                                               >  2:296519/1‑150 (MQ=255)
                               ccTTATCACACCTTTGAGCATGTGCTGGAACTGCTGCGTCAGGCTTCGTTCGCCCCGAGCGTACTGGCGATTAAAATTAACATTTACCGCGTGGCGAAAGATTCACGCATCATCGCCTCGATGCTCCACGCCGCACATAACGGTAAGaaa                                                                                                                      >  1:1547264/1‑150 (MQ=255)
                               ccTTATCACACCTTTGAGCATGTGCTGGAACTGCTGCGTCAGGCTTCGTTCGACCCGAGCGTACTGGCGATTAAAATTAACATTTACCGCGTGGCGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGaaa                                                                                                                      >  2:898924/1‑150 (MQ=255)
                                   aTCACACCTTTGAGCATGTGCTTGAACTGCTGCGTCAGGCTTCGTTCGACCCGAGCGTACTGGCGATTAAAATTAACATTTACCGCGTGGCGATAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCa                                                                                                                  <  2:645817/150‑1 (MQ=255)
                                    tCACACCTTTGAGCATGTGCTGGAACTGCTGCGTCAGGCTTCGTTCGACCCGAGCGTACTGGCGATTAAAATTAACATTTACCGCGTGGCGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCAc                                                                                                                 <  2:1485196/150‑1 (MQ=255)
                                     cacaCCTTTGAGCATGTGCTGGAACTGCTGCGTCAGGCTTCGTTCGACCCGAGCGTACTGGCGATTAAAATTAACATTTACCGCGTGGCGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCAcc                                                                                                                >  2:12406/1‑150 (MQ=255)
                                          cTTTGAGCATGTGCTGGAACTGCTGCGTCAGGCTTCGTTCGATCCGAGCGTACTGGCGATTAAAATTAACATTTACCGCGTGGCGAAAGATTCACGCATCATCGACTCGATGATCCCCGCCGCACATAACGGTAAGAAAGTCACCgtggg                                                                                                           >  2:438816/1‑149 (MQ=255)
                                          cTTTGAGCATGTGCTGGAACTGCTGCGTCAGGCTTCGTTCGACCCGAGCGTACTGGCGATTAAAATTAACATTTACCGCGTGGCGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCgtggt                                                                                                           <  1:586930/150‑1 (MQ=255)
                                          cTTTGAGCATGTGCTGGAACTGCTGCGTCAGGCTTCGTTCGACCCGAGCGTACTGGCGATTAAAATTAACATTTACCGCGTGGCGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCatggt                                                                                                           >  1:748934/1‑150 (MQ=255)
                                            ttGAGCATGTGCTGGAACTGCTGCGTCAGGCTTCGTTCGACCCGAGCGTACTGGCGATTAAAATTAACATTTACCGCGTGGCGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGtggtgg                                                                                                         <  2:1132092/150‑1 (MQ=255)
                                            ttGAGCATGTGCTGGAACTGCTGCGTCAGGCTTCGTTCGACCCGAGCGTACTGGCGATTAAAATTAACATTTACCGCGTGGCGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGtggtgg                                                                                                         <  1:330310/150‑1 (MQ=255)
                                            ttGAGCATGTGCTGGAACTGCTGCGTCAGGCTTCGTTCGACCCGAGCGTACTGGCGATTAAAATTAACATTTACCGCGTGGCGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGtggtgg                                                                                                         >  2:1616161/1‑150 (MQ=255)
                                             tGAGCATGTGCTGGAACTGCTGCGTCAGGCTTCGTTCGACCCGAGCGTACTGGCGATTAAAATTAACATTTACCGCGTGGCGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTggtggt                                                                                                        >  2:1044871/1‑150 (MQ=255)
                                               aGCATGTGCTGGAACTGCTGCGTCAGGCTTCGTTCGACCCGAGCGTACTGGCGATTAAAATTAACATTTACCGCGTGGCGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTg                                                                                                      <  2:1302605/150‑1 (MQ=255)
                                                 cATGTGCTGGAACTGCTGCGTCAGGCTTCGTTCGACCCGAGCGTACTGGCGATTAAAATTAACATTTACCGCGTGGCGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAg                                                                                                    <  1:1570318/150‑1 (MQ=255)
                                                  aTGTGCTGGAACTGCTGCGTCAGGCTTCGTTCGACCCGAGCGTACTGGCGATTAAAATTAACATTTACCGCGTGGCGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGt                                                                                                   <  1:954708/150‑1 (MQ=255)
                                                  aTGTGCTGGAACTGCTGCGTCAGGCTTCGTTCGACCCGAGCGTACTGGCGATTAAAATTAACATTTACCGCGTGGCGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGt                                                                                                   <  1:1332203/150‑1 (MQ=255)
                                                     tgCTGGAACTGCTGCGTCAGGCTTCGTTCGACCCGAGCGTACTGGCGATTAAAATTAACATTTACCGCGTGGCGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTAc                                                                                                <  1:1485854/150‑1 (MQ=255)
                                                      gCTGGAACTGCTGCGTCAGGCTTCGTTCGACCCGAGCGTACTGGCGATTAAAATTAACATTTACCGCGTGGCGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACa                                                                                               >  2:990384/1‑150 (MQ=255)
                                                      gCTGGAACTGCTGCGTCAGGCTTCGTTCGACCCGAGCGTACTGGCGATTAAAATTAACATTTACCGCGTGGCGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACa                                                                                               <  1:1235256/150‑1 (MQ=255)
                                                               gctgcGTCAGGCTTCGTTCGACCCGAGCGTACTGGCGATTAAAATTAACATTTACCGCGTGGCGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGttt                                                                                      >  2:194149/1‑150 (MQ=255)
                                                                ctgcGTCAGGCTTCGTTCGACCCGAGCGTACTGGCGATTAAAATTAACATTTACCGCGTGGCGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTc                                                                                     <  2:1061390/150‑1 (MQ=255)
                                                                 tgcGTCAGGCTTCGTTCGACCCGAGCGTACTGGCGATTAAAATTAACATTTACCGCGTGGCGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTcg                                                                                    <  1:34059/150‑1 (MQ=255)
                                                                  gcGTCAGGCTTCGTTCGACCCGAGCGTACTGGCGATTAAAATTAACTTTTCCCGGGGGGCGAAGGATTCCCGCAACTGCGCCGGGTGGATCCACGCCGCCCATAACGGTACGAAAGTCCCGGGGGGGGTTGATTTACAGGCGCGTTGcgc                                                                                   >  2:373197/1‑149 (MQ=255)
                                                                    gTCAGGCTTCGTTCGACCCGAGCGTACTGGCGATTAAAATTAACATTTACCGCGTGGCGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTcgacg                                                                                 >  2:1222477/1‑150 (MQ=255)
                                                                     tCAGGCTTCGTTCGACCCGAGCGTACTGGCGATTAAAATTAACATTTACCGCGTGGCGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTcgacga                                                                                <  1:489077/150‑1 (MQ=255)
                                                                     tCAGGCTTCGTTCGACCCGAGCGTACTGGCGATTAAAATTAACATTTACCGCGTGGCGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTcgacga                                                                                >  1:1438153/1‑150 (MQ=255)
                                                                      cAGGCTTCGTTCGACCCGAGCGTACTGGCGATTAAAATTAACATTTACCGCGTGGCGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTCGACgaa                                                                               >  1:1484633/1‑150 (MQ=255)
                                                                      cAGGCTTCGTTCGACCCGAGCGTACTGGCGATTAAAATTAACATTTACCGCGTGGCGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTCGACgaa                                                                               >  2:742420/1‑150 (MQ=255)
                                                                       aGGCTTCGTTCGACCCGAGCGTACTGGCGATTAAAATTAACATTTACCGCGTGGCGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTCGACgaag                                                                              <  1:1044871/150‑1 (MQ=255)
                                                                           ttcgttcgACCCGAGCGTACTGGCGATTAAAATTAACATTTACCGCGTGGCGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTCGACGAAGAAGc                                                                          <  1:1279167/150‑1 (MQ=255)
                                                                             cgttcgACCCGAGCGTACTGGCGATTAAAATTAACATTTACCGCGTGGCGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTCGACGAAGAAGCCa                                                                        <  1:675308/150‑1 (MQ=255)
                                                                             cgttcgACCCGAGCGTACTGGCGATTAAAATTAACATTTACCGCGTGGCGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTCGACGAAGAAGCCa                                                                        >  1:1534442/1‑150 (MQ=255)
                                                                               ttcgACCCGAGCGTACTGGCGATTAAAATTAACATTTACCGCGTGGCGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTCGACGAAGAAGCCAAc                                                                      <  1:898924/150‑1 (MQ=255)
                                                                               ttcgACCCGAGCGTACTGGCGATTAAAATTAACATTTACCGCGTGGCGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTCGACGAAGAAGCCAAc                                                                      >  2:598558/1‑150 (MQ=255)
                                                                                      cGAGCGTACTGGCGATTAAAATTAACATTTACCGCGTGGCGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTCGACGAAGAAGCCAACATTCACt                                                               >  2:591066/1‑150 (MQ=255)
                                                                                      cGAGCGTACTGGCGATTAAAATTAACATTTACCGCGTGGCGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTCGACGAAGAAGCCAACATTCACt                                                               <  1:12406/150‑1 (MQ=255)
                                                                                              cTGGCGATTAAAATTAACATTTACCGCGTGGCGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGTACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTCGACGAAGAAGCCAACATTCACTGGGCGAAg                                                       >  1:1467466/1‑150 (MQ=255)
                                                                                              cTGGCGATTAAAATTAACATTTACCGCGTGGCGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTCGACGAAGAAGCCAACATTCACTGGGCGAAg                                                       >  1:496847/1‑150 (MQ=255)
                                                                                              cTGGCGATTAAAATTAACATTTACCGCGTGGCGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTCGACGAAGAAGCCAACATTCACTGGGCGAAg                                                       >  1:696611/1‑150 (MQ=255)
                                                                                              cTGGCGATTAAAATTAACATTTACCGCGTGGCGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTCGACGAAGAAGCCAACATTCACTGGGCGAAg                                                       >  1:959557/1‑150 (MQ=255)
                                                                                                ggCGATTAAAATTAACATTTACCGCGTGGCGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCTTTTCGACGAAGAAGCCAACATTCACTGGGCGAAgcg                                                     <  2:1534442/150‑1 (MQ=255)
                                                                                                ggCGATTAAAATTAACATTTACCGCGTGGCGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTcga                                                                                   <  1:1064282/120‑1 (MQ=255)
                                                                                                ggCGATTAAAATTAACATTTACCGCGTGGCGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTcga                                                                                   >  2:1064282/1‑120 (MQ=255)
                                                                                                ggCGATAAAAATTAACATTTACCGAGTCGCCAAAGATTCACGCATCATCGACTCGATGATACACGCCCCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTCGACGAAGAAGCCAACATTCACTGGGCGAAgcg                                                     >  2:450833/1‑150 (MQ=255)
                                                                                                 gCGATTAAAATTAACATTTACCGCGTGGCGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTCGACGAAGAAGCCAACATTCACTGGGCGAAgcgc                                                    >  1:1276366/1‑150 (MQ=255)
                                                                                                  cGATTAAAATTAACATTTACCGCGTGGCGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGTACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTCGACGAAGAAGCCAACATTCACTGGGCGAAGCGcc                                                   <  2:1467466/150‑1 (MQ=255)
                                                                                                  cGATTAAAATTAACATTTACCGCGTGGCGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTCGACGAAGAAGCCAACATTCACTGGGCGAAGCGcc                                                   <  2:959557/150‑1 (MQ=255)
                                                                                                   gATTAAAATTAACATTTACCGCGTGGCGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTCGACGAAGAAGCCAACATTCACTGGGCGAAGCGCCt                                                  <  2:1484633/150‑1 (MQ=255)
                                                                                                   gATTAAAATTAACATTTACCGCGTGGAGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTCGACGAAGAAGCCAACATTCACTGGGCGAAGCGCCt                                                  >  2:210674/1‑150 (MQ=255)
                                                                                                     ttAAAATTAACATTTACCGCGTGGCGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTCGACGAAGAAGCCAACATTCACTGGGCGAAGCGCCTGa                                                >  2:502507/1‑150 (MQ=255)
                                                                                                       aaaaTTAACATTTACCGCGTGGCGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTCGACGAAGAAGCCAACATTCACTGGGCGAAGCGCCTGAcc                                              >  1:150763/1‑150 (MQ=255)
                                                                                                         aaTTAACATTTACCGCGTGGCGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTCGACGAAGAAGCCAACATTCACTGGGCGAAGCGCCTGACCGa                                            >  2:1052487/1‑150 (MQ=255)
                                                                                                          aTTAACATTTACCGCGTGGCGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTCGACGAAGAAGCCAACATTCACTGGGCGAAGCGCCTGACCGaa                                           >  1:80364/1‑150 (MQ=255)
                                                                                                          aTTAACATTTACCGCGTGGCGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCAACGTGGTGGTTGAGTTACAGGCGCGTTTCGACGAAGAAGCCAACATTCACTGGGCGAAGCGCCTGACCGaa                                           >  2:1082254/1‑150 (MQ=255)
                                                                                                                aTTTACCGCGTGGCGAAAGATTCACGCTTCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTCGACGAAGAAGCCAACATTCACTGGGCGAAGCGCCTGACCGAAGCAGGc                                     <  2:1547264/150‑1 (MQ=255)
                                                                                                                   tACCGCGTGGCGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTCGACGAAGAAGCCAACATTCACTGGGCGAAGCGCCTGACCGAAGCAGGCGTg                                  >  2:834938/1‑150 (MQ=255)
                                                                                                                     ccGCGTGGCGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTCGACGAAGAAGCCAACATTCACTGGGCGAAGCGCCTGACCGAAGCAGGCGTGCa                                <  1:750805/150‑1 (MQ=255)
                                                                                                                       gcgTGGCGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTCGACGAAGAAGCCAACATTCACTGGGCGAAGCGCCTGACCGAAGCAGGCGTGCACg                              <  2:1572016/150‑1 (MQ=255)
                                                                                                                           ggCGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTCGACGAAGAAGCCAACATTCACTGGGCGAAGCGCCTGACCGAAGCAGGCGTGCACGTTAt                          <  1:742420/150‑1 (MQ=255)
                                                                                                                           ggCGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTCGACGAAGAAGCCAACATTCACTGGGCGAAGCGCCTGACCGAAGCAGGCGTGCACGTTAt                          >  1:271926/1‑150 (MQ=255)
                                                                                                                                  gATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTCGACGAAGAAGCCAACATTCACTGGGCGAAGCGCCTGACCGAAGCAGGCGTGCACGTTATCTTctct                   <  1:297357/150‑1 (MQ=255)
                                                                                                                                    ttCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTCGACGAAGAAGCCAACATTCACTGGGCGAAGCGCCTGACCGAAGCAGGCGTGCACGTTATCTTCTCTgc                 >  2:1635956/1‑150 (MQ=255)
                                                                                                                                     gcgtggatcatcGAGTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTCGACGAAGAAGCCAACATTCACTGGGCGAAGCGCCTGACCGAAGCAGGCGTGCACGTTATCTTCTCTgcg                >  1:1451995/7‑150 (MQ=255)
                                                                                                                                      cACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCCCCGTGGTGGTTGAGTTACAGGCGCGTTTCGACGAAGAAGCCAACATTCACTGGGCGAAGCGCCTGACCGAAGCAGGCGTGCACGTTATCTTCTCTgcgc               >  2:67578/1‑150 (MQ=255)
                                                                                                                                      cACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTCGACGAAGAAGCCAACATTCACTGGGCGAAGCGCCTGACCGAAGCAGGCGTGCACGTTATCTTCTCTgcgc               >  1:1152860/1‑150 (MQ=255)
                                                                                                                                      cACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTCGACGAAGAAGCCAACATTCACTGGGCGAAGCGCCTGACCGAAGCAGGCGTGCACGTTATCTTCTCTgcgc               >  1:1620625/1‑150 (MQ=255)
                                                                                                                                       gtggatcatcGAGTCGATGATGCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTCGACGAAGAAGCCAACATTCACTGGGCGAAGCGCCTGACCGAAGCAGGCGTGCACGTTATCTTCTCTGCGcc              >  1:281022/5‑150 (MQ=255)
                                                                                                                                        tggatcatcGAGTAGATGATCAAAGCAGCACATAAAGGTAAGAAAGTCACCCTGGTTGTTTAGATAAAGGCGCGTTTCGACGAAGAAGCCAAAATTCACTGGGCGAAGCGCCTGACCGAAGCACGCGTGCACGTTATCTTCTCTGCGCCg             >  2:312068/4‑150 (MQ=255)
                                                                                                                                         gCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTCGACGAAGAAGCCAACATTCACTGGGCGAAGCGCCTGACCGAAGCAGGCGTGCACGTTATCTTCTCTGCGCCgg            >  1:9040/1‑150 (MQ=255)
                                                                                                                                              atcGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTCGACGAAGAAGCCAACATTCACTGGGCGAAGCGCCTGACCGAAGCAGGCGTGCACGTTATCTTCTCTGCGCCGGGGCTg       >  2:286578/1‑150 (MQ=255)
                                                                                                                                                 gACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTCGACGAAGAAGCCAACATTCACTGGGCGAAGCGCCTGACCGAAGCAGGCGTGCACGTTATCTTCTCTGCGCCGGGGCTGaaa    >  1:17358/1‑150 (MQ=255)
                                                                                                                                                   cTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTCGACGAAGAAGCCAACATTCACTGGGCGAAGCGCCTGACCGAAGCAGGCGTGCACGTTATCTTCTCTGCGCCGGGGCTGAAAAt  <  1:591066/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CATTCGCGAACGCGATGTGTTGCTCTATTATCCTTATCACACCTTTGAGCATGTGCTGGAACTGCTGCGTCAGGCTTCGTTCGACCCGAGCGTACTGGCGATTAAAATTAACATTTACCGCGTGGCGAAAGATTCACGCATCATCGACTCGATGATCCACGCCGCACATAACGGTAAGAAAGTCACCGTGGTGGTTGAGTTACAGGCGCGTTTCGACGAAGAAGCCAACATTCACTGGGCGAAGCGCCTGACCGAAGCAGGCGTGCACGTTATCTTCTCTGCGCCGGGGCTGAAAAT  >  NC_000913/2624048‑2624344

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