Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,466,336 C→G 9.2% 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.2% 212.7 / 7.0 71pseudogene (945/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base C (39/25);  new base G (3/4);  total (42/29)
Fisher's exact test for biased strand distribution p-value = 4.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGA  >  NC_000913/1466188‑1466485
                                                                                                                                                    |                                                                                                                                                     
ttGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCAccc                                                                                                                                                      <  2:603726/150‑1 (MQ=255)
  gATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAg                                                                                                                                                    <  2:1088044/150‑1 (MQ=255)
  gATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATAGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAg                                                                                                                                                    >  1:305817/1‑150 (MQ=255)
     atCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGAtt                                                                                                                                                 >  2:1250639/1‑150 (MQ=255)
     atCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGAtt                                                                                                                                                 <  2:630449/150‑1 (MQ=255)
     atCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGAtt                                                                                                                                                 >  2:772757/1‑150 (MQ=255)
      tCACCGTTGACAGCGAGACTGTGGATAACAAGGGCTCCAAGAACGTCACCGATCCGTAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATTACGAATGCGAGATCACCATAACCAACGGTGGCACCGGCACCCAGATTa                                                                                                                                                <  1:262378/150‑1 (MQ=18)
       cACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGAATaa                                                                                                                                               >  2:806715/1‑150 (MQ=255)
        aCCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCAATTGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCCCCGGCACCCAGATTAAt                                                                                                                                              >  2:65029/1‑150 (MQ=255)
         ccGGGGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATg                                                                                                                                             <  1:423351/150‑1 (MQ=255)
            gTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtg                                                                                                                                          <  1:457193/150‑1 (MQ=255)
                      gACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCAc                                                                                                                                <  2:183292/150‑1 (MQ=255)
                          gTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCCCCGGCACCCAGATTAATGGTGATGACGCCACGGCa                                                                                                                            >  2:1131170/1‑150 (MQ=255)
                           tGGATAACAAGGGCACCATGACCGTCACCGATACGGAGTACATGGGTATCCAGATCGACGGTGACAAGGCCATCGTAAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGGGATGCCGCCCCGGCaa                                                                                                                           >  1:1273953/1‑150 (MQ=18)
                             gATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAaac                                                                                                                         <  1:81306/150‑1 (MQ=255)
                              aTAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAaaca                                                                                                                        >  1:1150580/1‑150 (MQ=255)
                               tAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAaacaa                                                                                                                       >  1:733468/1‑150 (MQ=255)
                               tAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGAGAGAGCCCCATCACCACCGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAaacaa                                                                                                                       >  2:947195/1‑150 (MQ=255)
                                 aCAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCCCCGGCACCCAGATTAATGGTGATGACGCCACGGCAaacaaca                                                                                                                     >  2:910791/1‑150 (MQ=255)
                                 aCAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAaacaaca                                                                                                                     <  2:1424880/150‑1 (MQ=255)
                                 aCAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAaacaaca                                                                                                                     >  1:131954/1‑150 (MQ=255)
                                  cAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAaacaacaa                                                                                                                    >  1:1012652/1‑150 (MQ=255)
                                   aaGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGAATAATGGTGATGACGCCACGGCAaacaacaac                                                                                                                   >  2:1517454/1‑150 (MQ=255)
                                     gggCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGCGCACCATAACCCACGGTGGCACCGGGACCCCGATTAATTGTGATGACGCCACGGCAAACAACAACgg                                                                                                                 >  2:1654292/1‑150 (MQ=18)
                                     gggCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACgg                                                                                                                 <  1:1250639/150‑1 (MQ=255)
                                      ggCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGc                                                                                                                <  2:1564701/150‑1 (MQ=255)
                                       gCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCa                                                                                                               <  1:856383/150‑1 (MQ=255)
                                         aCCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCaaa                                                                                                             <  2:1464858/150‑1 (MQ=255)
                                                      ccGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCCCCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACg                                                                                                >  2:1294449/1‑150 (MQ=255)
                                                         aTCCGGAGTCCATGGGTATCCAGATCGACGCTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCTAAACCACCGTTGACGGCa                                                                                             <  1:928794/150‑1 (MQ=255)
                                                           gcGGAGTCCATGGGTATCCAGTTCGACGGTGAAAAGGCTTTAGTAAATAACTAATGCGAGAGCTCCATCACCAACGGTGGCACCGGCACGCAGATTAAAGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAg                                                                                           <  2:1509396/149‑1 (MQ=14)
                                                              gAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATAGTAAATAACGAAGGCGAGAGCACAATCAACAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCAAGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGAt                                                                                        >  1:1231636/1‑150 (MQ=17)
                                                               aGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGAtt                                                                                       <  1:880630/150‑1 (MQ=255)
                                                                gTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTc                                                                                      >  1:1150591/1‑150 (MQ=255)
                                                                   cATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAaacaac                                                                                                                      <  1:364731/115‑1 (MQ=255)
                                                                   cATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAAAGGTGATGACGCCACGGCAaacaac                                                                                                                      >  2:364731/1‑115 (MQ=255)
                                                                       ggTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGt                                                                               >  1:983577/1‑150 (MQ=255)
                                                                        gTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCCCCCAGAGTAATGGTGATGCCGCCCCGGCAAACAACAACGGCCAAACCACCGTGGACGGCAAGGCTTCCACCGGtc                                                                              >  2:1259403/1‑149 (MQ=18)
                                                                            ccAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGa                                                                          >  2:118449/1‑150 (MQ=255)
                                                                            ccAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGa                                                                          >  1:308254/1‑150 (MQ=255)
                                                                             cAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGaa                                                                         >  2:1345767/1‑150 (MQ=255)
                                                                              aGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGaaa                                                                        <  1:218452/150‑1 (MQ=255)
                                                                                aTCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATc                                                                      <  2:1803591/150‑1 (MQ=255)
                                                                                  cGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCCCCGGTACGGAAATCaa                                                                    >  1:867598/1‑150 (MQ=255)
                                                                                  cGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCaa                                                                    >  2:432112/1‑150 (MQ=255)
                                                                                   gACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGCGCACCATCACCAACGGTGGAACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAACCCACCGATGACGGCAAGGATTCCACCGGTACGGAAATCAAc                                                                   >  2:1232688/1‑150 (MQ=255)
                                                                                         gACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGtaat                                                             >  2:1067055/1‑150 (MQ=255)
                                                                                              ggCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACACCAACGGCAAATCCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACgg                                                        >  1:799910/1‑150 (MQ=255)
                                                                                              ggCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACgg                                                        >  1:1433885/1‑150 (MQ=255)
                                                                                                ccATCGTCAATAACGAAGGCGAGAGCTCCATCACCAACGGTGGCACCGGCACGCAGAt                                                                                                                                                  <  2:691055/58‑1 (MQ=14)
                                                                                                ccATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACGCAGAt                                                                                                                                                  >  1:691055/1‑58 (MQ=255)
                                                                                                  aTCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGaaaa                                                    >  2:1785865/1‑150 (MQ=255)
                                                                                                    cGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACGCAGATTAACGGtgat                                                                                                                                        >  1:1002132/1‑64 (MQ=255)
                                                                                                      tCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTa                                                >  1:13164/1‑150 (MQ=255)
                                                                                                         aTAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTc                                             <  2:1033726/150‑1 (MQ=255)
                                                                                                         aTAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTc                                             <  2:1796543/150‑1 (MQ=255)
                                                                                                          tAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAAAGGTGATGACGCCACGGCAAACAACAACGGCAAAACAACAGTTGACGGCAAGTATTCCACCGGTACGTAAATCAACTATAATAACGGAAAAGTTATTCa                                            >  2:336366/1‑150 (MQ=17)
                                                                                                             cGAAGGGGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGAAGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGa                                         >  2:496766/1‑150 (MQ=255)
                                                                                                             cGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGa                                         <  2:1150580/150‑1 (MQ=255)
                                                                                                              gAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGAc                                        >  1:1513739/1‑150 (MQ=255)
                                                                                                                aGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACgg                                      >  1:1769352/1‑150 (MQ=255)
                                                                                                                 ggCGAGAGCACCATCACCAACGGTGGCACCGGCACGCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCTAGGATTCCAcc                                                                                  <  1:199315/105‑1 (MQ=12)
                                                                                                                 ggCGAGAGCACCATCACCAACGGTGGCACCGGCACGCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCAcc                                                                                  >  2:199315/1‑105 (MQ=17)
                                                                                                                 ggCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGc                                     <  2:1150591/150‑1 (MQ=255)
                                                                                                                         cACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGa                             <  1:1371448/150‑1 (MQ=255)
                                                                                                                         cACCATCACCAACGGTGGCACCGGCACCCAGATAAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAGTTATTCAGGACGGCGATCTGGAt                            <  2:305817/150‑1 (MQ=255)
                                                                                                                          aCCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATAACAGGAAATAACGGAAACGTGATTCAGGACGGCGATCTGGAt                            >  2:969698/1‑150 (MQ=11)
                                                                                                                           ccATGACGAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAAGCACCGTTGTCGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATg                           <  1:806715/150‑1 (MQ=255)
                                                                                                                           ccATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATg                           >  1:8948/1‑150 (MQ=255)
                                                                                                                           ccATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATg                           <  1:705416/150‑1 (MQ=255)
                                                                                                                           ccATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATg                           <  1:1517454/150‑1 (MQ=255)
                                                                                                                                aCCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGCTACGGAAATCAACGCTAATAACGGAAAAGATATTCAGGACGGCGATCTGGATGTCAgc                      >  2:943572/1‑150 (MQ=255)
                                                                                                                                         ggcaccggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGt             <  1:1294449/150‑1 (MQ=255)
                                                                                                                                         ggcaccggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGt             >  1:1557070/1‑150 (MQ=255)
                                                                                                                                          gcaccggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTc            >  2:1025423/1‑150 (MQ=255)
                                                                                                                                          gcaccggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTc            >  2:1018095/1‑150 (MQ=255)
                                                                                                                                           caccggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTcc           >  2:1507066/1‑149 (MQ=255)
                                                                                                                                            cccggcaccCAGATTAATGGTGCTGACGCCACGGCAAACAACAACGGCAACACCACCGTTGACGGCAAGGATTCCACCGGTCCGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCAc          <  1:346807/149‑1 (MQ=255)
                                                                                                                                            accggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCAc          <  1:772757/150‑1 (MQ=255)
                                                                                                                                              cggcaccCAGATTATTGGGGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGAACGGAAATCACCGTTAATTACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACgg        <  1:336366/150‑1 (MQ=17)
                                                                                                                                                 caccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTAt     <  2:733468/150‑1 (MQ=255)
                                                                                                                                                   cccAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGcggcggcgg              <  1:1081762/139‑1 (MQ=255)
                                                                                                                                                   cccAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGcggcggcgg              >  2:1081762/1‑139 (MQ=255)
                                                                                                                                                   cccAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTg   >  2:1113399/1‑150 (MQ=255)
                                                                                                                                                    ccAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGa  >  1:1590544/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGA  >  NC_000913/1466188‑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: