Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,286,917 C→T 5.7% pseudogene (1998/2525 nt) yejO ← pseudogene, autotransporter outer membrane homology;putative transport; Not classified; putative ATP‑binding component of a transport system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,286,9170CT 5.7% 220.6 / 2.4 70pseudogene (1998/2525 nt)yejOpseudogene, autotransporter outer membrane homology;putative transport; Not classified; putative ATP‑binding component of a transport system
Reads supporting (aligned to +/- strand):  ref base C (27/39);  new base T (2/2);  total (29/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

AGATAGGGTCTAACACTCCACAATCCGTCAACCCAACGGAACCCGCTCTCAACATGAGCACCCGCGCCGTTACTATTGTAATCGCCAAACGCTGTTGCCCCATTACTCATCTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATGAGAGTAACCAAAGATCAAGCCGCGAATTGTACTGC  >  NC_000913/2286769‑2287060
                                                                                                                                                    |                                                                                                                                                
aGATAGGGGCTAACACTCCACAATCCGTCAACCCAACGGAACCCGCTCTCAACATGAGCACCCGCGCCGTTACTATTGTAATCGCCAAACGCTGTTGCCCCATTACTCATCTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTcacc                                                                                                                                                 <  1:1119153/150‑1 (MQ=255)
 gATAGGGTCTAACACTCCACAATCCGTCAACCCAACGGAACCCGCTCTCAACATGAGCACCCGCGCCGTTACTATTGTAATCGCCAAACGCTGTTGCCCCATTACTCATCTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTcacca                                                                                                                                                <  1:132530/150‑1 (MQ=255)
        tCTAACACTCCACAATCCGTCAACCCAACGGAACCCGCTCTCAACATGAGCACCCGCGCCGTTACTATTGTAATCGCCAAACGCTGTTGCCCCATTACTCATCTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATc                                                                                                                                         >  1:375903/1‑150 (MQ=255)
        tCTAACACTCCACAATCCGTCAACCCAACGGAACCCGCTCTCAACATGAGCACCCGCGCCGTTACTATTGTAATCGCCAAACGCTGTTGCCCCATTACTCATCTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATc                                                                                                                                         >  1:369777/1‑150 (MQ=255)
        tCTAACACTCCACAATCCGTCAACCCAACGGAACCCGCTCTCAACATGAGCACCCGCGCCGTTACTATTGTAATCGCCAAACGCTGTTGCCCCATTACTCATCTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATc                                                                                                                                         <  2:1487491/150‑1 (MQ=255)
        gcTAACACTCCACAATCCGTCAACCCAACGGAACCCGCTCTCAACATGAGCACCCGCGCCGTTACTATTGTAATCGCCAAACGCTGTTGCCCCATTACTCATCTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATc                                                                                                                                         <  2:1303136/149‑1 (MQ=255)
         cTAACACTCCACAATCCGTCAACCCAACGGAACCCGCTCTCAACATGAGCACCCGCGCCGTTACTATTGTAATCGCCAAACGCTGTTGCCCCATTACTCATCTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCa                                                                                                                                        <  1:660692/150‑1 (MQ=255)
          tAACACTCCACAATCCGTCAACCCAACGGAACCCGCTCTCAACATGAGCACCCGCGCCGTTACTATTGTAATCGCCAAACGCTGTTGCCCCATTACTCATCTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCaa                                                                                                                                       <  2:945466/150‑1 (MQ=255)
           aaCACTCCACAATCCGTCAACCCAACGGAACCCGCTCTCAACATGAGCACCCGCGCCGTTACTATTGTAATCGCCAAACGCTGTTGCCCCATTACTCATCTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAAc                                                                                                                                      >  1:430950/1‑150 (MQ=255)
           aaCACTCCACAATCCGTCAACCCAACGGAACCCGCTCTCAACATGAGCACCCGCGCCGTTACTATTGTAATCGCCAAACGCTGTTGCCCCATTACTCATCTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAAc                                                                                                                                      >  2:445841/1‑150 (MQ=255)
             cacTCCACAATCCGTCAACCCAACGGAACCCGCTCTCAACATGAGCACCCGCGCCGTTACTATTGTAATCGCCAAACGCTGTTGCCCCATTACTCATCTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCCCCCCCTCACCAt                                                                                                                                    >  1:1224453/1‑150 (MQ=255)
               cTCCACAATCCGTCAACCCAACGGAACCCGCTCTCAACATGAGCACCCGCGCCGTTACTATTGTAATCGCCAAACGCTGTTGCCCCATTACTCATCTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAg                                                                                                                                  <  1:722680/150‑1 (MQ=255)
                tCCACAATCCGTCAACCCAACGGAACCCGCTCTCAACATGAGCACCCGCGCCGTTACTATTGTAATCGCCAAACGCTGTTGCCCCATTACTCATCTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAgg                                                                                                                                 <  2:1589382/150‑1 (MQ=255)
                  cacaATCCGTCAACCCAACGGAACCCGCTCTCAACATGAGCACCCGCGCCGTTACTATTGTAATCGCCAAACGCTGTTGCCCCATTACTCATCTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCATCACCCCATCAACATAGGCa                                                                                                                               >  1:299418/1‑150 (MQ=255)
                  cacaATCCGTCAACCCAACGGAACCCGCTCTCAACATGAGCACCCGCGCCGTTACTATTGTAATCGCCAAACGCTGTTGCCCCATTACTCATCTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCATCACCCCATCAACATAGGCa                                                                                                                               >  1:346780/1‑150 (MQ=255)
                    caATCCGTCAACCCAACGGAACCCGCTCTCAACATGAGCACCCGCGCCGTTACTATTGTAATCGCCAAACGCTGTTGCCCCATTACTCATCTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCCACATAGGCAcc                                                                                                                             >  2:289729/1‑150 (MQ=255)
                    caATCCGTCAACCCAACGGAACCCGCTCTCAACATGAGCACCCGCGCCGTTACTATTGTAATCGCCAAACGCTGTTGCCCCATTACTCATCTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCAcc                                                                                                                             >  1:1345213/1‑150 (MQ=255)
                    caATCCGTAAACCCAAAGGAACCCGCTCTCAACATGTGCACCCGAGCCGTTACTATTGTAATAGCCAATCGCTGGTGCCCCATTACACATCTTGCCATTGATGGTGTTGGCAAAACGGTCTACTTTCACCACCCCATCTACATAGGCAca                                                                                                                             <  2:1093866/150‑2 (MQ=255)
                     aaTCCGTCAACCCAACGGAACCCGCTCTCAACATGAGCACCCGCGCCGTTACTATTGTAATCGCCAAACGCTGTTGCCCCATTACTCATCTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCg                                                                                                                            <  1:999220/150‑1 (MQ=255)
                          gTCAACCCAACGGAACCCGCTCTCAACATGAGCACCCGCGCCGTTACTATTGTAATCGCCAAACGCTGTTGCCCCATTACTCATCTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTg                                                                                                                       <  2:985053/150‑1 (MQ=255)
                          gTCAACCCAACGGAACCCGCTCTCAACATGAGCACCCGCGCCGTTACTATTGTAATCGCCAAACGCTGTTGCCCCATTACTCATCTTGCCACGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTg                                                                                                                       <  2:1727829/150‑1 (MQ=255)
                           tCAACCCAACGGAACCCGCTCTCAACATGAGCACCCGCGCCGTTACTATTGTAATCGCCAAACGCTGTTGCCCCATTACTCATCTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGa                                                                                                                      <  2:505804/150‑1 (MQ=255)
                            cAACCCAACGGAACCCGCTCTCAACATGAGCACCCGCGCCGTTACTATTGTAATCGCCAAACGCTGTTGCCCCATTACTCATCTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGAt                                                                                                                     <  1:445841/150‑1 (MQ=255)
                             aaCCCAACGGAACCCGCTCTCAACATGAGCACCCGCGCCGTTACTATTGTAATCGCCAAACGCTGTTGCCCCATTACTCATCTTGCCATGGATGGTGTTGGCAAAACGGGCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATg                                                                                                                    <  2:376831/150‑1 (MQ=255)
                              aCCCAACGGAACCCGCTCTCAACATGAGCACCCGCGCCTTTACTATTGTAATCGCCAAACGCTGTTGCCCCATTACTCATCTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGc                                                                                                                   <  2:1251663/150‑1 (MQ=255)
                               cccAACGGAACCCGCTCTCAACATGAGCACCCGCGCCGTTACTATTGTAATCGCCAAACGCTGTTGCCCCATTACTCATCTTGCCATGGATGGTGTTGGCAAAACGGTCACCTTTCACCACCCCATCAACATAGGCACAGTTCTGATGCt                                                                                                                  >  1:24772/1‑150 (MQ=255)
                                             ctctCAACATGAGCACCCGCGCCGTTACTATTGTAATCGCCAAACGCTGTTGCCCCATTACTCATCTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCATCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATaa                                                                                                    <  2:346780/150‑1 (MQ=255)
                                             ctctCAACATGAGCACCCGCGCCGTTACTATTGTAATCGCCAAACGCTGTTGCCCCATTACTCATCTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCATCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATaa                                                                                                    <  2:299418/150‑1 (MQ=255)
                                                tcAACATGAGCACCCGCGCCGTTACTATTGTAATCGCCAAACGCTGTTGCCCCATTACTCATCTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGcc                                                                                                 <  1:1386436/150‑1 (MQ=255)
                                                 cAACATGAGCACCCGCGCCGTTACTATTGTAATCGCCAAACGCTGTTGCCCCATTACTCATCTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACAACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGccc                                                                                                >  2:1056904/1‑150 (MQ=255)
                                                    cATGAGCACCCGCGCCGTTACTATTGTAATCGCCAAACGCTGTTGCCCCATTACTCATCTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCCa                                                                                            >  1:1023694/1‑150 (MQ=255)
                                                           aCCCGCGCCGTTACTATTGTAATCGCCAAACGCTGTTGCCCCATTACTCATCTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCCAGGGtata                                                                                     <  1:1481197/150‑1 (MQ=255)
                                                             ccGCGCCGTTACTATTGTAATCGCCAAACGCTGTTGCCCCATTACTCATCTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCCAGGGTATAGc                                                                                   <  1:1636437/150‑1 (MQ=255)
                                                             ccGCGCCGTTACTATTGTAATCGCCAAACGCTGTTGCCCCATTACTCATCTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCCAGGGTATAGc                                                                                   <  1:1558458/150‑1 (MQ=255)
                                                             ccGCGCCGTTACTATTGTAATCGCCAAACGCTGTTGCCCCATTACTCATCATGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTCATGCTCCCAACCCGCAAAAGCCCCCAGGGTATAGc                                                                                   >  2:1039347/1‑150 (MQ=255)
                                                                 gcCGTTACTATTGTAATCTCCAAACGCTGTTGCCCCATTACTCATCTTTCCATGGATGGTGTTGGCAAAACGGGCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACAGGCATAAGCCCCCAGGGTATAGCTATc                                                                               <  1:1468976/150‑1 (MQ=255)
                                                                    gTTACTATTGTAATCGCCAAACGCTGTTGCCCCATTACTCATCTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGCTGCTCCCAACCGGCATAAGCCCCCAGGGTATAGCTATCGat                                                                            >  2:1019847/1‑150 (MQ=255)
                                                                       aCTATTGTAATCGCCAAACGCTGTTGCCCCATTACTCATCTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCCAGGGTATAGCTATCGATAtt                                                                         <  1:1240680/150‑1 (MQ=255)
                                                                        cTATTGTAATCGCCAAACGCTGTTGCCCCATTACTCATCTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCCAGGGTATAGCTATCGATATTa                                                                        >  2:1725500/1‑150 (MQ=255)
                                                                           ttGTAATCTCCAAACGCTGTTGCGCCATTACTCATCTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCAACCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCCAGGGTATAGCTATCGATATTACCt                                                                     <  2:634047/150‑1 (MQ=255)
                                                                            tGTAATCGCCAAACGCTGTTGCCCCATTACTCATCTTGCCATGGATGGTGTGGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCCAGGGTATAGCTATCGATATTACCtt                                                                    <  1:335166/150‑1 (MQ=255)
                                                                                aTCGCCAAACGCTGTTGCCCCATTACTCATCTTGCCTTGGTTTGTGTTTGCATAACGGTCAACTTTCACCACCACATCAACATAGGCTAAGTTCTGATGCTCCCATCCGGCATAAGCCCGCTGGGTATAGCTATCGATATTACCTTTgcc                                                                <  2:701534/150‑1 (MQ=255)
                                                                                  cGCCAAACGCTGTTGCCCCATTACTCATCTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCCAGGGTATAGCTATCGATATTACCTTTgccgc                                                              >  1:1576745/1‑150 (MQ=255)
                                                                                   gCCAAACGCTGTTGCCCCATTACTCATCTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGTCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCCAGGGTATAGCTATCGATATTACCTTTgccgcc                                                             >  1:1309555/1‑150 (MQ=255)
                                                                                     cAAACGCTGTTGCCCCATTACTCATCTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCCAGGGTATAGCTATCGATATTACCTTTGCcgccgc                                                           >  2:1402750/1‑150 (MQ=255)
                                                                                        aCGCTGTTGCCCCATTACTCATCTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGAt                                                        <  2:333441/150‑1 (MQ=255)
                                                                                          gCTGTTGCCCCATTACTCATCTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCCACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCa                                                      <  1:289729/150‑1 (MQ=255)
                                                                                            tGTTGCCCCATTACTCATCTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCCAGGGTATAGCTAACGATATTACCTTTGCCGCCGCGATCaaa                                                    <  1:856652/150‑1 (MQ=255)
                                                                                                  cccATTACTCATCTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAat                                              <  2:474995/150‑1 (MQ=255)
                                                                                                   ccATTACTCATCTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCCCCGTTCTGATGCTCCCACCCGGCATACGCCCCCAGGGTATAGCTATCGATATTGCCTTTGCCGCCGCGATCCAAACCAatt                                             >  2:592000/1‑149 (MQ=255)
                                                                                                         cTCATCTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCTGCATAAGCCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGaa                                       <  1:1056904/150‑1 (MQ=255)
                                                                                                         cTCATCTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGaa                                       <  2:1310493/150‑1 (MQ=255)
                                                                                                           cATCTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATg                                     <  2:533928/150‑1 (MQ=255)
                                                                                                            aTCTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATgag                                   >  2:392793/1‑150 (MQ=255)
                                                                                                            aTCTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATgag                                   >  2:392788/1‑150 (MQ=255)
                                                                                                               ttGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATGagag                                 <  2:1224453/150‑1 (MQ=255)
                                                                                                                tGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCACAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATGagag                                 <  2:1712438/150‑1 (MQ=255)
                                                                                                                      ggatggTGTTGGCAAAACGGTCAACTTTCACCACCCAATCAACATAGGCAACGTTCTGATGGTCCCAACCGGCATACGCCCCCAGGGTTTAGCTATCGATATTACCTTTGCCGCCGCGAGCAAAACCAATATCAGAATGAGAGTAACCaa                          >  2:313786/1‑150 (MQ=255)
                                                                                                                       gatggTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATGAGAGTAACCaaa                         >  2:1686858/1‑150 (MQ=255)
                                                                                                                        atggTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATGAGAGTAACCAAAg                        <  2:375903/150‑1 (MQ=255)
                                                                                                                        atggTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATGAGAGTAACCAAAg                        <  2:369777/150‑1 (MQ=255)
                                                                                                                                ggCAAAACGGTCAACTTTCACCACCCCATCAACATAGTAACCGTTCTGATGCTCCCAACCGGCATAAGCCCCCAGGGTAAAGCTATCGATATTACCTATGCCGCCGCGATCAAAACCAATAACAGAATGCGAGTAACCAAAGATCAAGcc                >  2:1075379/1‑150 (MQ=255)
                                                                                                                                 gCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATGAGAGTAACCAAAGATCAAGCcg               <  1:492434/150‑1 (MQ=255)
                                                                                                                                   aaaaCGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATGAGAGTAACCAAAGATCAAGCcgcg             <  2:1309555/150‑1 (MQ=255)
                                                                                                                                        ggTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATGAGAGTAACCAAAGATCAAGCCGCGAATTg        <  1:807666/150‑1 (MQ=255)
                                                                                                                                         gTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATGAGAGTAACCAAAGATCAAGCCGCGAATTGt       >  1:455921/1‑150 (MQ=255)
                                                                                                                                         gTCAACTTTCACCACCCCATCAACATAGGCACCGATCTGATGCTCCCAACCGGCATAAGCCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATGAGAGTAACCAAAGATCAAGCCGCGAATTGt       >  1:594967/1‑150 (MQ=255)
                                                                                                                                           cAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAACCCAATATCAGAATGAGAGTAACCAAAGATCAAGCCGCGAAGTGTAc     >  1:1378907/1‑150 (MQ=255)
                                                                                                                                           cAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATGAGAGTAACCAAAGATCAAGCCGCGAATTGTAc     <  1:1686858/150‑1 (MQ=255)
                                                                                                                                              cTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCCAGGGTTTAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATGAGAGAAACCAAAGATCAAGCCGCGAATTGTACTGc  >  2:852373/1‑150 (MQ=255)
                                                                                                                                              cTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCCAGGGTATAGCTATCTATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATGAGAGTAACCAAAGATCAAGCCGCGAATTGTACTGc  >  1:575013/1‑150 (MQ=255)
                                                                                                                                              cTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATGAGAGTAACCAAAGATCAAGCCGCGAATTGTACTGc  <  2:1576575/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
AGATAGGGTCTAACACTCCACAATCCGTCAACCCAACGGAACCCGCTCTCAACATGAGCACCCGCGCCGTTACTATTGTAATCGCCAAACGCTGTTGCCCCATTACTCATCTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATGAGAGTAACCAAAGATCAAGCCGCGAATTGTACTGC  >  NC_000913/2286769‑2287060

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