Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,742,795 A→T 5.3% E138V (GAA→GTA)  yfiN → putative membrane‑anchored diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,742,7950AT 5.3% 337.1 / 2.3 116E138V (GAA→GTA) yfiNputative membrane‑anchored diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base A (57/50);  new base T (2/5);  total (61/55)
Fisher's exact test for biased strand distribution p-value = 2.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 121 total aligned reads displayed.

ATTTTCAACTGCAGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAG  >  NC_000913/2742646‑2742941
                                                                                                                                                     |                                                                                                                                                   
aTTATCAACAGCAGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGCATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCACAGAACCCATCATTCAGCCGATTACTCACAATGgaga                                                                                                                                                     >  2:1255279/1‑150 (MQ=255)
 ttttCAACTGCAGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGaa                                                                                                                                                    <  2:1250460/150‑1 (MQ=255)
  tttCAACTGCAGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATCGTGaaa                                                                                                                                                   >  2:1591517/1‑150 (MQ=255)
     cAACTGCAGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACCCCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGGTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAAACa                                                                                                                                                >  2:280501/1‑150 (MQ=255)
      aaCTGCAGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCAt                                                                                                                                               >  1:692896/1‑150 (MQ=255)
        cTGCAGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTg                                                                                                                                             <  1:1557616/150‑1 (MQ=255)
          gcagcaGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACATTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCTCCCATCATTCAGCCGATTCGTCACAATGGTTTCACCATTgtg                                                                                                                                           <  2:1781311/150‑3 (MQ=255)
          gCAGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGc                                                                                                                                           <  2:1304401/150‑1 (MQ=255)
            cgAAGTACGTGATAAGCAGCATAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTTTCACCATTgtgac                                                                                                                                         <  2:180848/149‑5 (MQ=255)
            aGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTTTCACCATTgtgac                                                                                                                                         <  1:997733/150‑5 (MQ=255)
            aGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTTTCACCATTgtgac                                                                                                                                         <  1:1575494/150‑5 (MQ=255)
            aGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGa                                                                                                                                         <  1:1221036/150‑1 (MQ=255)
             gAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGaa                                                                                                                                        >  2:736595/1‑150 (MQ=255)
              aaGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAg                                                                                                                                       <  2:1834549/150‑1 (MQ=255)
              aaGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAg                                                                                                                                       <  1:1834712/150‑1 (MQ=255)
               aGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGt                                                                                                                                      >  1:522045/1‑150 (MQ=255)
                gTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTa                                                                                                                                     >  1:795795/1‑150 (MQ=255)
                gTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTa                                                                                                                                     >  1:795789/1‑150 (MQ=255)
                   cGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCACTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGc                                                                                                                                  >  2:1238304/1‑150 (MQ=255)
                     tGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCtt                                                                                                                                >  1:1835388/1‑150 (MQ=255)
                     tGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCACCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCtt                                                                                                                                >  2:1748906/1‑150 (MQ=255)
                          agcagcAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCATAAccc                                                                                                                           >  1:944520/1‑149 (MQ=255)
                          agcagaAAAATATTCTTGATTCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAATAAAATGGATCTTCCATGACCCCATCATTCAGCCGATTCTTCACAATGGTGAAACCATTGGCGAAGTACGCTTAAcct                                                                                                                           <  1:243263/150‑2 (MQ=255)
                           gcagcaAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCCTCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGc                                                                                                                          >  1:1533591/1‑150 (MQ=255)
                           gcagcaAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACAATTGGCGAAGTACGCTTAACCGc                                                                                                                          >  1:429687/1‑150 (MQ=255)
                             agcaAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTc                                                                                                                        <  1:59648/150‑1 (MQ=255)
                              gcaAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTcg                                                                                                                       >  1:187925/1‑150 (MQ=255)
                               caAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCAGCATTCAGCCGATTCGTCACAATGGTGCAACCATTGGCGAAGTACGCTTAACCGCTcgc                                                                                                                      >  1:1237527/1‑150 (MQ=255)
                                aaaaTATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTcgcg                                                                                                                     >  1:258660/1‑150 (MQ=255)
                                  aaTATTCTGGCATCCTGGCATTACAGCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGAc                                                                                                                   <  1:1345699/150‑1 (MQ=255)
                                    tatTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCCTCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAg                                                                                                                 >  1:386246/1‑150 (MQ=255)
                                     atTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCATCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGt                                                                                                                >  2:1825300/1‑150 (MQ=255)
                                      ttCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATACGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGtt                                                                                                               >  2:1144213/1‑150 (MQ=255)
                                           gCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCACTc                                                                                                          >  1:1383440/1‑150 (MQ=255)
                                             aTCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAg                                                                                                        >  2:430425/1‑150 (MQ=255)
                                               ccTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGAAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGcc                                                                                                      >  1:288416/1‑150 (MQ=255)
                                                  ggCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGGACGCTTACCCGCTCGCGACAGTTCAATCAGCCAtt                                                                                                   >  2:1185218/1‑150 (MQ=255)
                                                       tACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATc                                                                                              <  2:381804/150‑1 (MQ=255)
                                                           cccGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGt                                                                                          >  2:676746/1‑150 (MQ=255)
                                                           cccGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGt                                                                                          >  1:813959/1‑150 (MQ=255)
                                                            ccGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTACAATCAGCCATTTTATCTGGtt                                                                                         >  1:938849/1‑150 (MQ=255)
                                                             cGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCAATGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGttt                                                                                        >  1:651501/1‑150 (MQ=255)
                                                                     tCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTtcgctcgc                                                                                <  2:1403623/150‑1 (MQ=255)
                                                                     tCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTtcgctcgc                                                                                <  2:692896/150‑1 (MQ=255)
                                                                     tCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTtcgctcgc                                                                                <  1:240725/150‑1 (MQ=255)
                                                                      ccAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGATGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGcc                                                                               >  2:1564286/1‑150 (MQ=255)
                                                                      ccAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGcc                                                                               <  1:673576/150‑1 (MQ=255)
                                                                       cAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCg                                                                              <  2:1629252/150‑1 (MQ=255)
                                                                            gATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTg                                                                         <  1:181000/150‑1 (MQ=255)
                                                                              tACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATACGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGAc                                                                       >  2:1137675/1‑150 (MQ=255)
                                                                               aCTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGAcc                                                                      <  2:1691710/150‑1 (MQ=255)
                                                                                cTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCACCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACACTTCAATCAGCCATTTTATCTGGTTTTCGATCGCCGTACTGACCg                                                                     >  2:557856/1‑150 (MQ=255)
                                                                                 tttCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCgg                                                                    >  1:1221920/1‑150 (MQ=255)
                                                                                 tttCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCgg                                                                    >  1:382400/1‑150 (MQ=255)
                                                                                    cAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTg                                                                 <  2:258660/150‑1 (MQ=255)
                                                                                    cAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTg                                                                 <  2:187925/150‑1 (MQ=255)
                                                                                       tAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTAt                                                              <  1:1585627/150‑1 (MQ=255)
                                                                                            tCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACAGGTTGTATTctgc                                                         >  2:698309/1‑150 (MQ=255)
                                                                                             cATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTctgct                                                        >  2:1362466/1‑150 (MQ=255)
                                                                                             cATAAGCAACTGGCTCTTCCCTGCCCCCATCATTCAGACGATTCGTCACAATGGTGAAAACATTTGCGAATTACGCTTAACCGCTGGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACAGGTTGTATTctgct                                                        >  2:1599314/1‑150 (MQ=255)
                                                                                              aTAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTctgctg                                                       <  1:493603/150‑1 (MQ=255)
                                                                                               tAAGCCACTGGCTCTTCCCTGCCACCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTgg                                                      >  1:1790250/1‑150 (MQ=255)
                                                                                                  gCCACTGGCTCTTCCCTGCCCCCATCCTTCAGCCGATTCGTCACAATGGTGTAACGATTGGCGCAGTACGCTTAACCGCTCGCGACAGGTCAATAAGCCATTTTATCTGGTTTTCGATCGCCGTACTGACCGGTTGTATTCTGCTGGCAt                                                   >  1:307329/1‑150 (MQ=255)
                                                                                                  gCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCAt                                                   >  1:1370620/1‑150 (MQ=255)
                                                                                                   ccACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATc                                                  <  2:1730214/150‑1 (MQ=255)
                                                                                                     aCTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGCCCGGTTGTATTCTGCTGGCagaag                                                >  2:471197/1‑146 (MQ=255)
                                                                                                     aCTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAg                                                <  1:604469/150‑1 (MQ=255)
                                                                                                     aCTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAg                                                <  2:527006/150‑1 (MQ=255)
                                                                                                      cTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAgg                                               <  2:240271/150‑1 (MQ=255)
                                                                                                        ggCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCa                                             >  2:1851133/1‑150 (MQ=255)
                                                                                                        ggCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGACCTGACCGGGTGTATTCTGCTGGCATCAGGCa                                             >  1:143654/1‑150 (MQ=255)
                                                                                                          ctctTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATc                                           >  2:1032560/1‑150 (MQ=255)
                                                                                                             ttCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCa                                        >  2:580703/1‑150 (MQ=255)
                                                                                                              tCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCaa                                       >  1:1085195/1‑150 (MQ=255)
                                                                                                              tCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCaa                                       >  2:1827525/1‑150 (MQ=255)
                                                                                                               cccTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAAt                                      <  1:761369/150‑1 (MQ=255)
                                                                                                                ccTGCCCCCATCCTTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCCATCAGCCATTTTATCTGGGTTTCGCTCGCCGTACTGACGGGTTGTATTCTGCTGGCATCAGGCATCGCACtt                                     >  2:1455866/1‑150 (MQ=255)
                                                                                                                    cccccATCATTCAGCCGATTCGTCACAATGGTGAACCCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTAccc                                 <  2:903333/150‑1 (MQ=255)
                                                                                                                     ccccGTCATTCAGCCGATTCGTCACAATGGTGAAAGCATTGGCGTAGTACGCTTAACCGCTCGCGACAGTGCATTCGGCCTTTTTATCTGGTTTTCGCTCGCCGTACTGAGCGGTTGGATTCTGCTGGCGTCAGGCATCGCAATTAcccg                                >  2:238037/1‑149 (MQ=255)
                                                                                                                        catcatTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTg                                                                 >  2:26241/1‑114 (MQ=255)
                                                                                                                        catcatTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTg                                                                 <  1:26241/114‑1 (MQ=255)
                                                                                                                             ttCAGCCGATTCGTCACACTGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGcc                        <  2:795795/150‑1 (MQ=255)
                                                                                                                             ttCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGcc                        <  2:386246/150‑1 (MQ=255)
                                                                                                                              tCAGCCGCTTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCa                       <  2:799002/150‑1 (MQ=255)
                                                                                                                              tCAGCCGATTCGTCACAATGGTGCAACCATTGGCGAAGTACGCTTACCCGCTCGCCTCTCTTCAATAAGCAATTTTTTATGGTTTTCGCTCGCCGTACTGAACGGTTGTATGCTGCTGGCATCAGGCAGCGCGATTACCCGCAacctcga                       >  1:297967/1‑146 (MQ=255)
                                                                                                                              tCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCa                       <  1:1569916/150‑1 (MQ=255)
                                                                                                                               cAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATATTATCTGGTTTTCGCTCGCCGTACTGCCCGGTTGTATTCTGCTGGCATCAGGCATCGCATTTTCCCTCACCCGCCAt                      >  1:232585/1‑150 (MQ=255)
                                                                                                                                 gCCGATTCGTCACAATGTTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCAttt                   <  1:1135221/150‑1 (MQ=255)
                                                                                                                                 gCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCAttt                   >  2:1216139/1‑150 (MQ=255)
                                                                                                                                   cGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGACATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCAGTTGc                 >  2:97379/1‑150 (MQ=255)
                                                                                                                                   cGAATCGTCACAATGGGGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGc                 <  2:1142287/150‑1 (MQ=255)
                                                                                                                                     aTTCGTCACACTGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGcac               <  1:980703/150‑1 (MQ=255)
                                                                                                                                     aTTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGGTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGcac               <  1:261090/150‑1 (MQ=255)
                                                                                                                                      ttCGTCACAATGGTGTAACCATTGGCGAAGTACGCTTAACCGCTCTCGACAGTTCAATCATCCATTTTATCTGGTTTTCTCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGcaca              <  1:1591517/150‑1 (MQ=255)
                                                                                                                                          tCACAATGGTGAAACCATTGGCGAAGTACGCTTACCCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATgg          >  2:1327943/1‑150 (MQ=255)
                                                                                                                                          gcacaATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTGGCGAGAGGTCAATAAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATAAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATgg          <  2:1128590/149‑1 (MQ=255)
                                                                                                                                           gccaATGGTTTCACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCAATTGCACAATGGc         >  2:658632/3‑150 (MQ=255)
                                                                                                                                            acaATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGcc        >  2:376663/1‑150 (MQ=255)
                                                                                                                                              aaTGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTGTGCAAATGGcggg      <  2:1328532/150‑4 (MQ=255)
                                                                                                                                                  gTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAg  <  2:478291/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
ATTTTCAACTGCAGAAGTACGTGATAAGCAGCAAAATATTCTGGCATCCTGGCATTACACCCGTAAGGATCCAGGCGATACTTTCAGTAATTTCATAAGCCACTGGCTCTTCCCTGCCCCCATCATTCAGCCGATTCGTCACAATGGTGAAACCATTGGCGAAGTACGCTTAACCGCTCGCGACAGTTCAATCAGCCATTTTATCTGGTTTTCGCTCGCCGTACTGACCGGTTGTATTCTGCTGGCATCAGGCATCGCAATTACCCTCACCCGCCATTTGCACAATGGCCTGGTAG  >  NC_000913/2742646‑2742941

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: