Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,466,336 C→G 9.8% pseudogene (945/2513 nt) ydbA → pseudogene, autotransporter homolog; interrupted by IS2 and IS30

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,466,3360CG 9.8% 193.2 / 3.0 73pseudogene (945/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base C (35/30);  new base G (4/4);  total (39/34)
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

ATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGA  >  NC_000913/1466187‑1466485
                                                                                                                                                     |                                                                                                                                                     
aTTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTggcaccggcacc                                                                                                                                                       <  1:811211/150‑1 (MQ=255)
   gATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAg                                                                                                                                                    <  2:1347436/150‑1 (MQ=255)
   gATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAg                                                                                                                                                    <  2:266875/150‑1 (MQ=255)
    aaatCACCGGGGTCATCGCGACGGTGGCAAACCACGGCACCAAGACCGTCACCGATCCCGAGACCAAGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGa                                                                                                                                                   <  1:450949/148‑1 (MQ=17)
      atCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGACGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGAtt                                                                                                                                                 >  2:502016/1‑150 (MQ=255)
        cACCGGTGGCAGCGCGACGGTGGATACCAAGGGCACCATGACCGTCACCGATCCGGAGTCCAGGGGTATCCAGCGCGCCGGTGCCAAGGCCATCGTCAATAACGAAGGCGAGGGCACCATCACCAACGGTGGCACCGGCACCCAGATaaa                                                                                                                                               >  1:134221/1‑150 (MQ=17)
         aCCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAAt                                                                                                                                              <  1:1813510/150‑1 (MQ=255)
         aCCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAAt                                                                                                                                              <  2:831554/150‑1 (MQ=255)
          cctgttaCAGCGGGACGGTGGATAACACCGGCAGCATGACCGTTACCGATCCGGAGTCCATGGGTAGCCATATCGACGGTGACAAGGCAATCGTCAAAAACGAAGGAGAGAGCACGATCACCAACGGTGGCACCGGCACCCAGATTAAtt                                                                                                                                             <  2:201657/144‑2 (MQ=17)
          ccGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATg                                                                                                                                             <  1:1382534/150‑1 (MQ=255)
          ccGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATg                                                                                                                                             <  2:1362462/150‑1 (MQ=255)
            ggTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGt                                                                                                                                           <  1:1989276/150‑1 (MQ=255)
            ggTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGGCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGt                                                                                                                                           <  2:35835/150‑1 (MQ=255)
             gTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtg                                                                                                                                          >  1:1353984/1‑150 (MQ=255)
             gTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtg                                                                                                                                          <  2:549574/150‑1 (MQ=255)
             gTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGGCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAAGGATGGCTAGAGCACCATCACCTACGGTGGCACCGGCACCCAGATTAATGGtg                                                                                                                                          <  2:829183/150‑1 (MQ=18)
              ggACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtga                                                                                                                                         <  1:383349/149‑1 (MQ=255)
               gACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACTGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtgat                                                                                                                                        <  2:483444/150‑1 (MQ=255)
                aCAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCTGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtgatg                                                                                                                                       <  1:1652510/150‑1 (MQ=255)
                 cAGCGCAACGGTGGATAACAAGGACACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCCACGGTGGCACCGGCCCCCAGAGTTAGGGtgctgc                                                                                                                                      >  2:343927/1‑149 (MQ=14)
                   gcgcGACGGTGGATAACAAGGGCACCATCACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACg                                                                                                                                    <  1:1004295/150‑1 (MQ=255)
                      cgACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCCCCCAGAGTAATGGTGATGACGcac                                                                                                                                 >  1:682633/1‑148 (MQ=255)
                           gTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCa                                                                                                                            >  1:1143517/1‑150 (MQ=255)
                             ggATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGAGACGGCCATCGTCAATAACGAAGGCGAGAGAACCATCACCCTCGGTGGCCGCGGGACCGAGCTTAAGGGTGTGGACGCGCCGGCaaa                                                                                                                          >  2:128199/1‑150 (MQ=17)
                               aTAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATAGTAAAAAACGAAGGCGAGAGCCCCATCACCAACGGTGGCACCGGCACGCAGATTAAAGGTGATGACGCCACGGCAaaca                                                                                                                        >  2:760501/1‑150 (MQ=14)
                               aTAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCAAGTTCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGTGAGCACAATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAaaca                                                                                                                        <  2:453712/150‑1 (MQ=255)
                                  aCAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTGTACATATCGAAGGTGACAAGGCCATCGTCAATAACTAAGTATACATCCCCATCACCCACGTTGGCACCGGCACCCAGATTTATTGTGATGCATCCCCGGAAaacaaaa                                                                                                                     >  2:350202/1‑148 (MQ=12)
                                  aCAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAaacaaca                                                                                                                     >  1:1813736/1‑150 (MQ=255)
                                  aCAAGGGCACCATGACCGTCACCGATCCGGAGACCATGGGTATCCAGATCGACGGTGACAAGGCCATAGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAAAGGTGATGACGCCACCGCAaacaaca                                                                                                                     >  1:1402549/1‑150 (MQ=255)
                                    aaGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAaacaacaac                                                                                                                   <  2:554088/150‑1 (MQ=255)
                                      gggCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCCTCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACgg                                                                                                                 >  2:394138/1‑150 (MQ=255)
                                      gggCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATAGACGGGGACAAGGCCCTAGTCAATAACGAAGGCGAGAGCGCCATCACCAACGGTGGCCCCGGCACCCAGATTAAGGGTGATGACGCCACGGCAAACAACAACgg                                                                                                                 >  2:21099/1‑150 (MQ=14)
                                                       ccGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATTGTCAATAACGAAGGCGAGAGCTCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGtgatg                                                                                                                                       <  2:462846/111‑1 (MQ=17)
                                                       ccGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACg                                                                                                <  2:988765/150‑1 (MQ=255)
                                                       ccGATCCGGAGTCCATGGGTATCCAGATAGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCTCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGtgatg                                                                                                                                       >  1:462846/1‑111 (MQ=17)
                                                         gagccGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGGCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGc                                                                                              <  2:834832/147‑1 (MQ=255)
                                                         gATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCCCCCAGAATAATGGTGATGACGCCCCGGCAAACAACAACGGCAAACCCACCCTTGACGGc                                                                                              >  2:717698/1‑150 (MQ=18)
                                                         gATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGAATAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGc                                                                                              >  1:486314/1‑150 (MQ=255)
                                                           tCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCTCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAaaca                                                                                                                        <  1:848854/122‑1 (MQ=255)
                                                           tCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAaaca                                                                                                                        >  2:848854/1‑122 (MQ=255)
                                                               gAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGAt                                                                                        <  2:1752874/150‑1 (MQ=255)
                                                                   ccATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCCCCATCACCAACGGTGGCACCGGCACCCAGATTATTGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTccc                                                                                    >  1:94762/1‑149 (MQ=255)
                                                                    cATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACAATCACCAACGGTGGCGCCGGCACCCAGAGTAATGGTGATGACGCCCCGGCAAACAACAAATGCAAAACCACCGTTGACGGCAAGGATTGCAc                                                                                   >  2:1460788/1‑150 (MQ=18)
                                                                      tGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTCGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCg                                                                                 >  2:26619/1‑150 (MQ=255)
                                                                      tGGGGATCCAGATGGACGGTGACTAGGCCTTCGTCTCTCAGGCAGCATAGAGGACGATCTACAACCTTGGCACCGTGAACCAGATTACTGGTGATTAGGCCACTGTAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCg                                                                                 <  2:620400/150‑1 (MQ=9)
                                                                         gTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTa                                                                              >  1:1509029/1‑150 (MQ=255)
                                                                          tATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTAc                                                                             >  1:1577011/1‑150 (MQ=255)
                                                                            tCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCTCCATCACCAACGGTGGCACCGGCACGCAGATTAAGGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGAt                                                                                        <  2:1172729/137‑1 (MQ=255)
                                                                            tCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACGCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGAt                                                                                        <  2:1909751/137‑1 (MQ=255)
                                                                            tCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGAt                                                                                        >  1:1909751/1‑137 (MQ=255)
                                                                            tCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAAGGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGAt                                                                                        >  1:1172729/1‑137 (MQ=255)
                                                                             ccAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGCCGCCCCGGCCAACACCCCCGGCAACACCACCGTTGCCGGCAAGGCTTCCACCGGTCCGGa                                                                          >  2:1283126/1‑150 (MQ=18)
                                                                             ccAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCAACAACGGATGCACCTGCACGAAGTTAAATGTTGATTACGCCACGGCAAACAACAAAGGCATAAACACAGATGAAAGCAAGGATTCCACCGGTACggg                                                                          >  2:2008178/1‑149 (MQ=17)
                                                                              cAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGaa                                                                         >  2:991594/1‑150 (MQ=255)
                                                                                 ttCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGATAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATc                                                                      <  2:1806987/149‑1 (MQ=255)
                                                                                    gtcGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAATACCACCGTTGACGGCAAGGATTCCACCGTTACGGAAATCAAc                                                                   <  1:91807/148‑1 (MQ=255)
                                                                                    gACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCACCGGTGGCACCGGCACCCAGTTTAATGGTGATGACGCCCCGGCAACCAACACCGGCAAAACCCCCGTTGACGGCAAGGCTAGCACCGGTACGGAAATctcc                                                                   >  1:371347/1‑147 (MQ=25)
                                                                                    gACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAAc                                                                   >  1:1130641/1‑150 (MQ=255)
                                                                                    gACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAAc                                                                   >  1:1424815/1‑150 (MQ=255)
                                                                                          gACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCCCCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGACATCAACGGtaat                                                             >  1:350490/1‑150 (MQ=255)
                                                                                            cAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTacggaaataaacggaaataa                                                           >  2:1953118/1‑150 (MQ=255)
                                                                                            cAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGtaataa                                                           >  1:896408/1‑150 (MQ=255)
                                                                                             aaGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTTATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCGACCGGTACGGAAATCAACGGTAATAAc                                                          <  2:70109/150‑1 (MQ=255)
                                                                                                gCCATCGTCAATAACGAAGGCGAGAGCACCATCACCACCGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGa                                                       <  2:958979/150‑1 (MQ=255)
                                                                                                        cAATAACGAAGGCGAGAGCACCATCACCCACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGCCAAGGATTCCACCGGTACGGAAATCAACGGAAATAACGGAAAAGatat                                               >  1:1835641/1‑150 (MQ=255)
                                                                                                                  ggCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGAAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGc                                     >  2:1795793/1‑150 (MQ=255)
                                                                                                                         gCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTgg                              >  2:417815/1‑150 (MQ=255)
                                                                                                                         ccACCATCACCACCGCTGGCACCGGCACCCAGATTAATCGTGATGACGCCACGGCAAACAACAACCGCAAAACCACCGTTGACGGCAAGGATTCCAGCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTgg                              <  1:1984376/149‑1 (MQ=17)
                                                                                                                          cACCATCACCAACGGTGGCACCGGCACCCAGATTAATTGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGa                             <  2:2009955/150‑1 (MQ=255)
                                                                                                                            ccATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATg                           <  1:394138/150‑1 (MQ=255)
                                                                                                                                    aaCGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAaac                                                                                                                         <  2:517130/48‑1 (MQ=255)
                                                                                                                                    aaCGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAaac                                                                                                                         >  1:517130/1‑48 (MQ=255)
                                                                                                                                      cGGTGGAACCGGAAACCATAATAAAGAAGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGCCGGCAAGGATTCCCCCGGTACGGAAATTAAAGGGAATAACGGAAAAGTGATTCAGGACGGCGATCTGGATGTCAgcggcgg                 >  2:290897/1‑150 (MQ=11)
                                                                                                                                         tGGCACCGGCACCCAGATTAATGGTGATGACGCCACTGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCTAACGAAATAACGGAAATGTTATTCAGGACGGCGATATTGATGTCAGaggcggcgg              >  2:980659/1‑150 (MQ=14)
                                                                                                                                          ggcaccggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGt             <  1:502016/150‑1 (MQ=255)
                                                                                                                                          ggcaccggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGt             >  2:1142895/1‑150 (MQ=255)
                                                                                                                                          ggcaccggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGt             <  1:343927/150‑1 (MQ=255)
                                                                                                                                            caccggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCAcc                                                                                  >  1:113324/1‑79 (MQ=39)
                                                                                                                                            caccggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCAcc                                                                                  <  2:113324/79‑1 (MQ=39)
                                                                                                                                            caccggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTcc           >  2:1025194/1‑149 (MQ=255)
                                                                                                                                            caccggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTAGACGGCAAGGATGCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCa           <  1:128199/150‑1 (MQ=255)
                                                                                                                                              ccggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAg                                                                                           >  2:154345/1‑68 (MQ=40)
                                                                                                                                              ccggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAg                                                                                           <  1:154345/68‑1 (MQ=40)
                                                                                                                                              ccggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAATGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACg         <  1:608726/150‑1 (MQ=255)
                                                                                                                                              ccggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTTCGGATATAAACGGTAAGTACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACg         >  1:1823972/1‑150 (MQ=17)
                                                                                                                                                  caccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTAt     >  2:1010584/1‑150 (MQ=255)
                                                                                                                                                    cccAGATTAATGGTGATGACGCCACTGCAAACATCAACTGCAAAACCACCGGTGACGGCAATGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAACGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTg   <  2:134221/150‑1 (MQ=17)
                                                                                                                                                    cccAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCCCCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGAGGGCGATCTGGATGTCCGCGGCGGCGGTCACGGTATTg   >  2:348913/1‑150 (MQ=255)
                                                                                                                                                    cccAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTg                                                                                                   >  1:1693748/1‑54 (MQ=40)
                                                                                                                                                    cccAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTg                                                                                                   <  2:1693748/54‑1 (MQ=40)
                                                                                                                                                    cccAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAAGGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGTATGTCAGCGGCGGCGGTCACGGGATTg   >  2:1209625/1‑150 (MQ=255)
                                                                                                                                                    cccAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGAAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTg   >  2:298046/1‑150 (MQ=255)
                                                                                                                                                    cccAGATTAATGGTGATGAAGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCAACGGTACGGAAAACAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTg   >  1:977967/1‑150 (MQ=255)
                                                                                                                                                     ccAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGa  <  1:26619/150‑1 (MQ=255)
                                                                                                                                                     ccAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGa  >  1:846037/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGA  >  NC_000913/1466187‑1466485

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: