Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,223,930 T→G 6.1% intergenic (+23/+349) ycgI → / ← minE pseudogene/cell division topological specificity factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,223,9300TG 6.1% 236.8 / 5.6 82intergenic (+23/+349)ycgI/minEpseudogene/cell division topological specificity factor
Reads supporting (aligned to +/- strand):  ref base T (35/42);  new base G (2/3);  total (37/45)
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.13e-01

CCAGGCAGCTATTAAAGCCTCATTCAGCAGTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACA  >  NC_000913/1223784‑1224077
                                                                                                                                                  |                                                                                                                                                   
ccAGGCAGCTATTAAAGCCTCATTCAGCAGTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCgg                                                                                                                                                  >  1:1212582/1‑150 (MQ=255)
 cAGGCAGCTATTAAAGCCTCATTCAGCAGTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGgc                                                                                                                                                 >  1:98496/1‑150 (MQ=255)
 cAGGCAGCTATTAAAGCCTCATTCAGCAGTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAAAGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGgc                                                                                                                                                 >  2:860617/1‑150 (MQ=255)
    gCAGCTATTAAAGCCTCATTCAGCAGTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCa                                                                                                                                              >  1:993352/1‑150 (MQ=255)
      aGCTATTAAAGCCTCATTCAGCAGTACGTTTAGCGGACATCTCGGGGTGGGGGATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGt                                                                                                                                            <  1:736538/150‑1 (MQ=255)
       gCTATTAAAGCCTCATTCAGCAGTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCCGAAGCGCAGtt                                                                                                                                           <  2:1152924/150‑1 (MQ=255)
         tATTAAAGCCTCATTCAGCAGTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGtttt                                                                                                                                         <  2:406011/150‑1 (MQ=255)
          aTTAAAGCCTCATTCAGCAGTACGTTTAGCGGCCATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGCCCATCAACGATTAAACTGCGCTTCGGCGCAGttttg                                                                                                                                        >  2:1002474/1‑149 (MQ=255)
           ttAAAGCCTCATTCAGCAGTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAACCTGCGCTTCGGCGCAGTGTTCg                                                                                                                                       >  2:469138/1‑150 (MQ=255)
           ttAAAGCCTCATTCAGCAGTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCCGAAGCGCAGtttaat                                                                                                                                       <  1:899494/150‑4 (MQ=17)
           ttAAAGCCTCATTCAGCAGTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCCGAAGCGCAGtttaat                                                                                                                                       <  1:470357/150‑4 (MQ=17)
             aaaGCCTCATTCAGCAGTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGtt                                                                                                                                     <  2:715934/150‑1 (MQ=255)
             aaaGCCTCATTCAGCAGTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGtt                                                                                                                                     >  1:1673828/1‑150 (MQ=255)
              aaGCCTCATTCAGCAGTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGttt                                                                                                                                    <  2:641841/150‑1 (MQ=255)
                gCCTCATTCAGCAGTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTAc                                                                                                                                  >  1:1016264/1‑150 (MQ=255)
                 ccTCATTCAGCAGTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACa                                                                                                                                 <  1:790574/150‑1 (MQ=255)
                 ccTCATTCAGCAGTACGTGAAGCGGACATCTCGGGGTGGGGTATCGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACa                                                                                                                                 <  1:200142/150‑1 (MQ=255)
                      ttCAGCAGTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGAGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATg                                                                                                                            >  2:1412093/1‑150 (MQ=255)
                      ttCAGCAGTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCACGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATg                                                                                                                            <  1:849383/150‑1 (MQ=255)
                        cagcagTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGtt                                                                                                                          >  1:1591185/1‑150 (MQ=255)
                           cagTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGaa                                                                                                                       <  1:398322/150‑1 (MQ=255)
                              tACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAgg                                                                                                                    <  1:685969/150‑1 (MQ=255)
                                    tAGCGGACAGCTCGGGGTGGGGTATAGCCCTGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAAt                                                                                                              <  1:1199015/150‑1 (MQ=255)
                                       cGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCt                                                                                                           >  2:832206/1‑150 (MQ=255)
                                       cGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCt                                                                                                           >  1:1410882/1‑150 (MQ=255)
                                         gACATGTCGGGGTGCGCTATAGCGATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTAGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGCGGTTGAAAGGGAAAATTCTgg                                                                                                         <  2:1285337/150‑1 (MQ=25)
                                         gACATCTCGGGGTGTGGTATAGCCATGGTGCCGGGGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTgg                                                                                                         <  1:1246950/150‑1 (MQ=255)
                                          aCATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTggg                                                                                                        <  2:98496/150‑1 (MQ=255)
                                           cATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTgggg                                                                                                       >  2:1565586/1‑150 (MQ=255)
                                                      ggggtATATCCATGGGTCCGGTGTGGTATCCCCGTTGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATAAACGATTAAACTTCGCTTCGGCGCATTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGcc                                                                                            <  2:718794/150‑1 (MQ=37)
                                                        ggtATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCg                                                                                          >  1:1264826/1‑150 (MQ=255)
                                                         gtATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGc                                                                                         >  2:1321655/1‑150 (MQ=255)
                                                          tataGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGcc                                                                                        <  1:1606731/150‑1 (MQ=255)
                                                            taGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAg                                                                                      >  1:358875/1‑150 (MQ=255)
                                                             aGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAGCCCGCCAGtt                                                                                    >  2:803550/1‑150 (MQ=255)
                                                                   ggTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTAcggc                                                                               <  1:1065272/150‑1 (MQ=255)
                                                                    gTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCATTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTAcggcg                                                                              <  2:162999/150‑1 (MQ=255)
                                                                    gTGCCGGTGTGGAATCCCCGTGGAACGCGGGGGCTGGTGTGAACTGGTCATTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTAcggcg                                                                              <  2:986295/150‑1 (MQ=255)
                                                                      gCCGGTGGGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCggg                                                                            <  1:1577852/150‑1 (MQ=255)
                                                                       ccGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGa                                                                           <  2:1673828/150‑1 (MQ=255)
                                                                        cGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGaa                                                                          <  2:1402670/150‑1 (MQ=255)
                                                                        cGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGaa                                                                          >  1:1002682/1‑150 (MQ=255)
                                                                         ggTGTGGAATCCCCGTGGAACGCGGTGGCTGGGGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGaaa                                                                         <  1:1565586/150‑1 (MQ=255)
                                                                          gtgtGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAAc                                                                        >  1:141437/1‑150 (MQ=255)
                                                                            gtgGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCt                                                                      <  1:1085338/150‑1 (MQ=255)
                                                                            gtgGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCAGCCAGTTACGGCGGGAAACCt                                                                      >  2:840292/1‑150 (MQ=255)
                                                                             tgGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCAGCCAGTTACGGCGGGAAACCTc                                                                     >  1:331275/1‑150 (MQ=255)
                                                                                 aTCCCCGTGGAACGCGGGGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATcc                                                                 <  1:643076/150‑1 (MQ=255)
                                                                                    cccGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTCCGGCGGGAACCCTCATCCTAt                                                              >  2:292312/1‑150 (MQ=255)
                                                                                      cGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATgg                                                            <  1:925866/150‑1 (MQ=255)
                                                                                       gTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCGAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATggg                                                           <  2:993352/150‑1 (MQ=255)
                                                                                       gTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATggg                                                           <  1:1004386/150‑1 (MQ=255)
                                                                                         ggAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTAAAGGATGATGAAAGTGAAAATTCTGGGTCAAAAAAAGCCCGCAAGTTACGGCGGGAAACCTCATCCTATGGgag                                                         >  1:1107689/1‑150 (MQ=39)
                                                                                           aaCGCGGTGGCTGGTGGGAACTGGTCGGTCTGACCATCAACGATTAAAATGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGaa                                                       <  2:1128577/150‑1 (MQ=255)
                                                                                                    gCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGaata                                              <  2:1003888/150‑1 (MQ=255)
                                                                                                     cTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGaataa                                             >  1:622749/1‑150 (MQ=255)
                                                                                                       ggTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATg                                           <  1:886514/150‑1 (MQ=255)
                                                                                                       ggTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATg                                           >  1:225115/1‑150 (MQ=255)
                                                                                                       ggTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATg                                           <  1:847240/150‑1 (MQ=255)
                                                                                                          gtgAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGaaa                                        <  1:91048/150‑1 (MQ=255)
                                                                                                          gtgAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGaaa                                        >  1:452808/1‑150 (MQ=255)
                                                                                                           tgAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAAt                                       <  2:1212582/150‑1 (MQ=255)
                                                                                                            gAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAAtt                                      >  1:392464/1‑150 (MQ=255)
                                                                                                               cTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCg                                   <  1:923583/150‑1 (MQ=255)
                                                                                                               cTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCg                                   >  1:701589/1‑150 (MQ=255)
                                                                                                               cTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAAATCTGGGGCAAAAAAAGCCATCCAGTTAATGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCg                                   >  1:1671742/1‑150 (MQ=255)
                                                                                                                      ttCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTat                            <  2:1264826/150‑1 (MQ=255)
                                                                                                                       tCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATCAAATTGCGGCGATatc                           >  2:1638913/1‑150 (MQ=255)
                                                                                                                        cTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCTGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTatca                          <  2:331275/150‑1 (MQ=255)
                                                                                                                        cTGACCATCAACGATTAAACTGCGCGTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAACAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATATTGAAATTGCGGGGTTatca                          <  2:312562/150‑1 (MQ=255)
                                                                                                                            ccATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCAtctc                      >  2:1656105/1‑150 (MQ=255)
                                                                                                                            ccATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCAtctc                      <  1:1474417/150‑1 (MQ=255)
                                                                                                                             cATCAACGATTAAACTGCGATTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAACCCTCATCCTATGGGAGAACAATGAATAATGATATTGCGGGGTTATCATCTcc                     >  2:141168/1‑150 (MQ=255)
                                                                                                                                   cGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGtat               <  2:1016264/150‑1 (MQ=255)
                                                                                                                                    gATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGtata              <  2:1376995/150‑1 (MQ=255)
                                                                                                                                    gATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGtata              >  2:1260731/1‑150 (MQ=255)
                                                                                                                                     cgaaaaCTGCGCCGAAGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTatat             >  2:1712638/4‑150 (MQ=255)
                                                                                                                                     cgaaaaCTGCGCCGAAGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTAACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTatat             >  1:86906/4‑150 (MQ=255)
                                                                                                                                       tAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATAGCGGGGTTATCATCTCCCAGTATATcc           <  2:1410882/150‑1 (MQ=255)
                                                                                                                                        aaaCTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCa          >  2:26402/1‑150 (MQ=255)
                                                                                                                                               gcTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCACCCAGTATATCCATACTAAc   <  2:849005/150‑1 (MQ=255)
                                                                                                                                                cTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACa  >  1:1604940/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CCAGGCAGCTATTAAAGCCTCATTCAGCAGTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAGCCCGCCAGTTACGGCGGGAAACCTCATCCTATGGGAGAACAATGAATAATGAAATTGCGGGGTTATCATCTCCCAGTATATCCATACTAACA  >  NC_000913/1223784‑1224077

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