Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,084,178 C→G 5.3% L578L (CTG→CTC) ‡ speA ← arginine decarboxylase, biosynthetic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,084,1780CG 5.3% 259.7 / 4.2 76L578L (CTG→CTC) ‡speAarginine decarboxylase, biosynthetic
Reads supporting (aligned to +/- strand):  ref base C (39/33);  new base G (2/2);  total (41/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.57e-01

GATCGCGGAACTGGGTTAACAGCGTTTTCGGATCGAGCTGTACATATTGCAGCATGTCCGCCACGGTATCGCCTTCGTCAGACAGTTCTACTTCTACGCTACCGTCAGGGAAGACGAACACGTCAACCGCTTCGGTATCACCGAACAGGTTGTGCATGTTGCCGAGGATCTCCTGATATGCGCCGACCATAAAGAAACCGAGCATCGGCGGATTCTCTGGATCGTACTCCGGCATTGGCATTGTCGTGGCAATACCGTCACCATCAATATAGTGGTCGATAGCACCGTCAG  >  NC_000913/3084033‑3084323
                                                                                                                                                 |                                                                                                                                                 
gATCGCGGAACTGGGTTAACAGCGTTTTCGGATCGAGCTGTACATATTGCAGCATGTCCGCCACGGTATCGCCTTCGTCAGACAGTTCTACTTCTACGCTACCGTCAGGGAAGACGAACACGTCAACCGCTTCGGTATCACCGAACAGGt                                                                                                                                               <  1:370551/150‑1 (MQ=255)
 aTCGCGGAACTGGGTTAACAGCGTTTTCGGATCGAGCTGTACATATTGCAGCATGTCCGCCACGGTATCGCCTTCGTCAGACAGTTCTACTTCTACGCTACCGTCAGGGAAGACGAACACGTCAACCTGTTCGGTGATACCGAAGCGGtt                                                                                                                                              <  1:1049490/150‑1 (MQ=255)
 aTCGCGGAACTGGGTTAACAGCGTTTTCGGATCGAGCTGTACATATTGCAGCATGTCCGCCACGGTATCGCCTTCGTCAGACAGTTCTACTTCTACGCTACCGTCAGGGAAGACGAACACGTCAACCGCTTCGGTATCACCGAACAGGtt                                                                                                                                              >  2:1047435/1‑150 (MQ=255)
   cgcgGAACTGGGTTAACAGCGTTTTCGGATCGAGCTGTACATATTGCAGCATGTCCGCCACGGTATCGCCTTCGTCAGACAGTTCTACTTCTACGCTACCGTCAGGGAAGACGAACACGTCAACCGCTTCGGTATCACCGAACAGGTtgt                                                                                                                                            <  2:556065/150‑1 (MQ=255)
    gcgGAACTGGGTTAACAGCGTTTTCGGATCGAGCTGTACATATTGCAGCATGTCCGCCACGGTATCGCCTTCGTCAGACAGTTCTACTTCTACGCTACCGTCAGGGAAGACGAACACGTCAACCGCTTCGGTATCACCGAACAGGTtgtg                                                                                                                                           >  1:1392585/1‑150 (MQ=255)
      ggAACTGGGTTAACAGCGTTTTCGGATCGAGCTGTACATATTGCAGCATGTCCGCCACGGTATCGCCTTCGTCAGACAGTTCTACTTCTACGCTACCGTCAGGGAAGACGAACACGTCAACCGCTTCGGTATCACCGAACAGGTTGTGCa                                                                                                                                         <  2:1030671/150‑1 (MQ=255)
          cTGGGTTAACAGCGTTTTCGGATCGAGCTGTACATATTGCAGCATGTCCGCCACGGTATCGCCTTCGTCAGACAGTTCTACTTCTACGCTACCTTCAGGGAAGACGAACACGTCAACCGCTTCGGTATCACCGAACAGGTTGTGCATGtt                                                                                                                                     >  2:767696/1‑150 (MQ=255)
           tGGGTTAACAGCGTTTTCGGATCGAGCTGTACATATTGCAGCATGTCCGCCACGGTATCGCCTTCGTCAGACAGTTCTACTTCTACGCTACCGTCAGGGAAGACGAACACGTCAACCGCTTCGGTATCACCGAACAGGTTGTGCATGTTg                                                                                                                                    <  2:742203/150‑1 (MQ=255)
            gggTTAACAGCGTTTTCGGATCGAGCTGTACATATTGCAGCATGTCCGCCACGGTATCGCCTTCGTCAGACAGTTCTACTTCTACGCTACCGTCAGGGAAGACGAACACGTCAACCGCTTCGGTATCACCGAACAGGTTGTGCATGTTGc                                                                                                                                   >  2:1367988/1‑150 (MQ=255)
              gTTAACAGCGTTTTCGGATCGAGCTGTACATATTGCAGCATGTCCGCCACGGTATCGCCTTCGTCAGACAGTTCTACTTCTACGCTACCGTCATGGAAGACGAACACGTCAACCGCTTCGGTATCACCGAACAGCTTGAGCATGTTGCCg                                                                                                                                 >  2:555561/1‑150 (MQ=255)
                tAACAGCGTTTTCGGATCGAGCTGTACATATTGCAGCATGTCCGCCACGGTATCGCCTTCGTCAGACAGTTCTACTTCTACGCTACCGTCAGGGAAGACGAACACGTCAACCGCTTCGGTATCACCGAACAGGTTGTGCATGTTGCCGAg                                                                                                                               >  2:338399/1‑150 (MQ=255)
                tAACAGCGTTTTCGGATCGAGCTGTACATATTGCAGCATGTCCGCCACGGTATCGCCTTCGTCAGACAGTTCTACTTCTACGCTACCGTCAGGGAAGACGAACACGTCAACCGCTTCGGTATCACCGAACAGGTTGTGCATGTTGCCGAg                                                                                                                               <  1:675817/150‑1 (MQ=255)
                tAACAGCGTTTTCGGATCGAGCTGTACATATTGCAGCATGTCCGCCACGGTATCGCCTTCGTCAGACAGTTCTACTTCTACGCTACCGTCAGGGAAGACGAACACGTCAACCGCTTCGGTATCACCGAACAGGTTGTGCATGTTGCCGAg                                                                                                                               <  1:1603753/150‑1 (MQ=255)
                  aCAGCGTTTTCGGATCGAGCTGTACATATTGCAGCATGTCCGCCACGGTATCGCCTTCGTCAGACAGTTCTACTTCTACGCTACCGTCAGGGAAGACGAACACGTCAACCGCTTCGGTATCACCGAACAGGTTGTGCATGTTGCCGAGGa                                                                                                                             >  1:1178826/1‑150 (MQ=255)
                   cAGCGTTTTCGGATCGAGCTGTACATATTGCAGCATGTCCGCCACGGTATCGCCTTCGTCAGACAGTTCTACTTCTACGCTACCGTCAGGGAAGACGAACACGTCAACCGCTTCGGTATCACCGAACAGGTTGTGCATGTTGCCGAGGAt                                                                                                                            >  2:433530/1‑150 (MQ=255)
                         tttCGGATCGAGCTGTACATATTGCAGCATGTCCGCCACGGTATCGCCTTCGTCAGACAGTTCTACTTCTACGCTACCGGCAGGGAAGACGAACACGTCAACCGCTTCGGTATCACCGAACAGGTTGTGCATGTTGCCGAGGATCTCCTg                                                                                                                      <  2:762267/150‑1 (MQ=255)
                           tCGGATCGAGCTGTACATATTGCAGCATGTCCGCCACGGTATCGCCTTCGTCAGACAGTTCTACTTCTACGCTACCGTCAGGGAAGACGAACACGTCAACCGCTTCGGTATCACCGAACAGGTTGTGCATGTTGCCGAGGATCTCCTGat                                                                                                                    >  2:1293845/1‑150 (MQ=255)
                            cGGATCGAGCTGTACATATTGCAGCATGTCCGCCACGGTATCGCCTTCGTCAGACAGTTCTACTTCTACGCTACCGTCAGGGAAGACGAACACGTCAACCGCTTCGGTATCACCGAACAGGTTGTGCTTGTTGCCGAGGATCTCCTGata                                                                                                                   >  2:1224564/1‑150 (MQ=255)
                             ggATCGAGCTGTACATATTGCAGCATGTCCGCCACGGTATCGCCTTCGTCAGACAGTTCTACTTCTACGCTACCGTCAGGGAAGACGAACACGTCAACCGCTTCGGTATCACCGAACAGGTTGTGCATGTTGCCGAGGATCTCCTGatat                                                                                                                  >  2:71101/1‑150 (MQ=255)
                             ggATCGAGCTGTACATATTGCAGCATGTCCGCCACGGTATCGCCTTCGTCAGACAGTTCTACTTCTACGCTACCGTCAGGGAAGACGAACACGTCAACCGCTTCGGTATCACCGAACAGGTTGTGCATGTTGCCGAGGATCTCCTGatat                                                                                                                  >  2:405777/1‑150 (MQ=255)
                                tCGAGCTGTACATATTGCAGCATGTCCGCCACGGTATCGCCTTCGTCAGACAGTTCTACTTCTACGCTACCGTCAGGGAAGACGAACACGTCAACCGCTTCGGTATCACCGAACAGGTTGTGCATGTTGCCGAGGATCTCCTGATATgcg                                                                                                               >  1:126689/1‑150 (MQ=255)
                                 cGAGCTGTACATATTGCAGCATGTCCGCCACGGTATCGCCTTCGTCAGACAGTTCTACTTCTACGCTACCGTCAGGGAAGACGAACACGTCAACCGCTTCGGTATCACCGAACAGGTTGTGCATGTTGCCGAGGATCTCCTGATATgcgc                                                                                                              <  1:1598798/150‑1 (MQ=255)
                                  gAGCTGTACATATTGCAGCATGTCCGCCACGGTATCGCCTTCGTCAGACAGTTCTACTTCTACGCTACCGTCAGGGAAGACGAACACGTCAACCGCTTCGGTATCACCGCACAGGTTGTGCATGTTGCCGAGGATCTCCTGATATGCGcc                                                                                                             >  1:1457220/1‑150 (MQ=255)
                                  gAGCTGTACATATTGCAGCATGTCCGCCACGGTATCGCCTTCGTCAGACAGTTCTACTTCTACGCTACCGTCAGGGAAGACGAACACGTCAACCGCTTCGGTATCACCGAACAGGTTGTGCATGTTGCCGAGGATCTCCTGATATGCGcc                                                                                                             >  2:921335/1‑150 (MQ=255)
                                    gCTGTACATATTGCAGCATGTCCGCCACGGTATCGCCTTCGTCAGACAGTTCTACTTCTACGCTACCGTCAGGGAAGACGAACACGTCAACCGCTTCGGTATCACCGAACAGGTTGTGCATGTTGCCGAGGATCTCCTGATATGCGCCGa                                                                                                           <  2:1359797/150‑1 (MQ=255)
                                      tGTACATATTGCAGCATGTCCGCCACGGTATCGCCTTCGTCAGACAGTTCTACTTCTACGCTACCGTCAGGGAAGACGAACACGTCAACCGCTTCGGTATCACCGAACAGGTTGTGCATGTTGCCGAGGATCTCCTGATATGCGCCGAcc                                                                                                         <  2:1392585/150‑1 (MQ=255)
                                       gTACATATTGCAGCATGTCCGCCACGGTATCGCCTTCGTCAGACAGTTCTACTTCTACGCTACCGTCAGGGAAGACGAACACGTCAACCGCTTCGGTATCACCGAACAGGTTGTGCATGTTGCCGAGGATCTCCTGATATGCGCCGACCa                                                                                                        <  1:363106/150‑1 (MQ=255)
                                         aCATATTGCAGCATGTCCGCCACGGTATCGCCTTCGTCAGACAGTTCTACTTCTACGCTACCGTCAGGGAAGACGAACACGTCAACCGCTTCGGTATCACCGAACAGGTTGTGCATGTTGCCGAGGATCTCCTGATATGCGCCGACCATa                                                                                                      <  1:1621693/150‑1 (MQ=255)
                                            tatTGCAGCATGTCCGCCACGGTATCGCCTTCGTCAGACAGTTCTACTTCTACGCTACCGTCAGGGAAGACGAACACGTCAACCGCTTCGGTATCACCGAAc                                                                                                                                                   >  2:879673/1‑102 (MQ=255)
                                            tatTGCAGCATGTCCGCCACGGTATCGCCTTCGTCAGACAGTTCTACTTCTACGCTACCGTCAGGGAAGACGAACACGTCAACCGCTTCGGTATCACCGAAc                                                                                                                                                   <  1:879673/102‑1 (MQ=255)
                                             atTGCAGCATGTCCGCCACGGTATCGCCTTCGTCAGACAGTTCTACTTCTACGCTACCGTCAGGGAAGACGAACACGTCAACCTGTTCGGTGATACCGAAGCGGTtg                                                                                                                                             >  2:511103/1‑107 (MQ=38)
                                             atTGCAGCATGTCCGCCACGGTATCGCCTTCGTCAGACAGTTCTACTTCTACGCTACCGTCAGGGAAGACGAACACGTCAACCTGTTCGGTGATACCGAAGCGGTtg                                                                                                                                             <  1:511103/107‑1 (MQ=38)
                                             atTGCAGCATGTCCGCCACGGTATCGCCTTCGTCAGACAGTTCTACTTCTACGCTACCGTCAGGGAAGACGAACACGTCAACCGCTTCGGTATCACCGAACAGGTTGTGCATGTTGCCGAGGATCTCCTGATATGCGCCGACCATAAAGa                                                                                                  <  2:912684/150‑1 (MQ=255)
                                                      tGTCCGCCACGGTATCGCCTTCGTCAGACAGTTCTACTTCTACGCTACCGTCAGGGAAGACGAACACGTCAACCGCTTCGGTATCACCGAACAGGTTGTGCATGTTGCCGAGGATCTCCTGATATGCGCCGACCATAAAGAAACCGAGCa                                                                                         <  1:622719/150‑1 (MQ=255)
                                                          cGCCACGGTATCGCCTTCGTCAGACAGTTCTACTTCTACGCTACCGTCAGGGAAGACGAACACGTCAACCGCTTCGGTATCACCGAACAGGTTGTGCATGTTGCCGAGGATCTCCTGATATGCGCCGACCATAAAGAAACCGAGCATcgg                                                                                     <  1:730270/150‑1 (MQ=255)
                                                            ccACGGTATCGCCTTCGTCAGACAGTTCTACTTCTACGCTACCGTCAGGGAAGACGAACACGTCAACCGCTTCGGTATCACCGAACAGGTTGTGCATGTTGCCGAGGATCTCCTGATATGCGCCGACCATAAAGAAACCGAGCATcggcg                                                                                   <  1:1293845/150‑1 (MQ=255)
                                                                        cTTCGTCAGACAGTTCTACTTCTACGCTACCGTCAGGGAAGACGAACACGTCAACCGCTTCGGTATCACCGAACAGGTTGTGCATGTTGCCGAGGATCTCCTGATATGCGCCGACCATAAAGAAACCGAGCATCGGCGGATTCTCTGGAt                                                                       <  1:921335/150‑1 (MQ=255)
                                                                             tCAGACAGTTCTACTTCTACGCTACCGTCAGGGAAGACGAACACGTCAACCGCTTCGGTATCACCGAACAGGTTGTGCATGTTGCCGAGGATCTCCTGATATGCGCCGACCATAAAGAAACCGAGCATCGGCGGATTCTCTGGATCGTAc                                                                  >  2:1552235/1‑150 (MQ=255)
                                                                               aGACAGTTCTACTTCTACGCTACCGTCAGGGAAGACGAACACGTCAACCGCTTCGGTATCACCGAACAGGTTGTGCATGTTGCCGAGGATCTCCTGATATGCGCCGACCATAAAGAAACCGAGCATCGGCGGATTCTCTGGATCGTACTc                                                                >  2:104123/1‑150 (MQ=255)
                                                                                   aGTTCTACTTCTACGCTACCGTCAGGGAAGACGAACACGTCAACCGCTTCGGTATCACCGAACAGGTTGTGCATGTTGCCGAGGATCTCCTGATATGCGCCGACCATAAAGAAACCGAGCATCGGCGGATTCTCTGGATCGTACTCCggc                                                            <  2:1178826/150‑1 (MQ=255)
                                                                                    gTTCTACTTCTACGCTACCGTCAGGGAAGACGAACACGTCAACCGCTTCGGTATCACCGAACAGGTTGTGCATGTTGCCGAGGATCTCCTGATATGCGCCGACCATAAAGAAACCGAGCATCGGCGGATTCTCTGGATCGTACTCCggca                                                           >  1:580613/1‑150 (MQ=255)
                                                                                        tacttctacGCTACCGTCAGGGAAGACGAACACGTCAACCGCTTCGGTATCACCGAACAGGTTGTGCATGTTGCCGAGGATCTCCTGATATGCGCCGACCATAAAGAAACCGAGCATCGGCGGATTCTCTGGATCGTACTCCggcattgg                                                       <  2:126689/150‑1 (MQ=255)
                                                                                        tacttctacGCTACCGTCAGGGAAGACGAACACGTCAACCGCTTCGGTATCACCGAACAGGTTGTGCATGTTGCCGAGGATCTCCTGATATGCGCCGACCATAAAGAAACCGAGCATCGGCGGATTCTCTGGATCGTACTCCggcattgg                                                       <  1:599095/150‑1 (MQ=255)
                                                                                         ccttctacGCTCCCGTCAGGGAAGACGAACACGTCAACCGCTTCGGTATCACCGAACAGGTTGTGCATGTTGCCGAGGATCTCCTGATATGCGCCGACCATAAAGAAACCGAGCATCGGCGGATTCTCTGGATCGTACTCCggcattggc                                                      <  1:1188753/149‑1 (MQ=255)
                                                                                          cttctacGCTACCGTCAGGGAAGACGAACACGTCAACCGCTTCGGTATCACCGAACAGGTTGTGCATGTTGCCGAGGATCTCCTGATATGCGCCGACCATAAAGAAACCGAGCATCGGCGGATTCTCTGGATCGTACTCCggcattggca                                                     >  2:1189090/1‑150 (MQ=255)
                                                                                             ctacGCTACCGTCAGGGAAGACGAACACGTCAACCGCTTCGGTATCACCGAACAGGTTGTGCATGTTGCCGAGGATCTCCTGATATGCGCCGACCATAAAGAAACCGAGCATCGGCGGATTCTCTGGATCGTACTCCGgcattggcattg                                                  >  1:978519/1‑150 (MQ=255)
                                                                                             ctacGCTACCGTCAGGGAAGACGAACACGTCAACCGCTTCGGTATCACCGAACAGGTTGTGCATGTTGCCGAGGATCTCCTGATATGCGCCGACCATAAAGAAACCGAGCATCGGCGGATTCTCTGGATCGTACTCCGgcattggcattg                                                  <  2:371207/150‑1 (MQ=255)
                                                                                                cGCTACCGTCAGGGAAGACGAACACGTCAACCGCTTCGGTATCACCGAACAGGTTGTGCATGTTGCCGAGGATCTCCTGATATGCGCCGACCATAAAGAAACCGAGCATCGGCGGATTCTCTGGATCGTACTCCGGCATTGGCATTGTCg                                               <  1:71101/150‑1 (MQ=255)
                                                                                                  cTACCGTCAGGGAAGACGAACACGTCAACCGCTTCGGTATCACCGAACAGGTTGTGCATGTTGCCGAGGATCTCCTGATATGCGCCGACCATAAAGAAACCGAGCATCGGCGGATTCTCTGGATCGTACTCCGGCATTGGCATTGTCGTg                                             >  1:579906/1‑150 (MQ=255)
                                                                                                    aCCGTCACGGAATACGAACACGTCAACCGCTTCGGTATCACCGAACAGGTTGTGCATGTTGCCGAGGATCTCCTGATATGCGCCGACCATAAAGAAACCGAGCATCGGCGGATTCTCTGGATCGTACTCCGGCATTGGCATTGTCGTGGc                                           >  2:56004/1‑150 (MQ=255)
                                                                                                     ccGTCAGGGAAGACGAACACGTCAACCGCTTCGGTATCACCGAACAGGTTGTGCTTGTTGCCGAGGATCTCCTGATATGCGCCGACCATAAAGAAACCGAGCATCGGCGGATTCTCTGGATCGTACTCCGGCATTGGCATTGTCGTGGCa                                          <  1:1224564/150‑1 (MQ=255)
                                                                                                     ccGTCAGGGAAGACGAACACGTCAACCGCTTCGGTATCACCGAACAGGTTGTGCATGTTGCCGAGGATCTCCTGATATGCGCCGACCATAAAGAAACCGAGCATCGGCGGATTCTCTGGATCGTACTCCGGCATTGGCATTGTCGTGGCa                                          >  2:189952/1‑150 (MQ=255)
                                                                                                         cAGGGAAGACGAACACGTCAACCGCTTCGGTATCACCGAACAGGTTGTGCATGTTGCCGAGGATCTCCTGATATGCGCCGACCATAAAGAAACCGAGCATCGGCGGATTCTCTGGATCGTACTCCGGCATTGGCATTGTCGTGGCAATAc                                      >  1:88005/1‑150 (MQ=255)
                                                                                                           gggAAGACGAACACGTCAACCGCTTCGGTATCACCGAACAGGTTGTGCATGTTGCCGAGGATCTCCTGATATGCGCCGACCATAAAGAAACCGAGCATCGGCGGATTCTCTGGATCGTACTCCGGCATTGGCATTGTCGTgg                                            <  2:1548160/142‑1 (MQ=255)
                                                                                                           gggAAGACGAACACGTCAACCGCTTCGGTATCACCGAACAGGTTGTGCATGTTGCCGAGGATCTCCTGATATGCGCCGACCATAAAGAAACCGAGCATCGGCGGATTCTCTGGATCGTACTCCGGCATTGGCATTGTCGTgg                                            >  1:1548160/1‑142 (MQ=255)
                                                                                                           gggAAGACGAACACGTCAACCGCTTCGGTATCACCGAACAGGTTGTGCATGTTGCCGAGGATCTCCTGATATGCGCCGACCATAAAGAAACCGAGCATCGGCGGATTCTCTGGATCGTACTCCGGCATTGGCATTGTCGTGGCAATACCg                                    >  2:1401796/1‑150 (MQ=255)
                                                                                                            ggAAGACGAACACGTCAACCGCTTCGGTATCACCGAACAGGTTGTGCATGTTGCCGAGGATCTCCTGATATGCGCCGACCATAAAGAAACCGAGCATCGGCGGATTCTCTGGATCGTACTCCGGCATTGGCATTGTCGTGGCAATACCGt                                   <  2:579906/150‑1 (MQ=255)
                                                                                                             gAAGACGAACACGTCAACCGCTTCGGTATCACCGAACAGGTTGTGCATGTTGCCGAGGATCTCCTGATATGCGCCGACCATAAAGAAACCGAGCATCGGCGGATTCTCTGGATCGTACTCCGGCATTGGCATTGTCGTGGCAATACCGTc                                  <  1:1047435/150‑1 (MQ=255)
                                                                                                                  cGAACACGTCAACCGCTTCGGTATCACCGAACAGGTTGTGCATGTTGCCGAGGATCTCCTGATATGCGCCGACCATAAAGAAACCGAGCATCGGCGGa                                                                                 >  2:461806/1‑98 (MQ=255)
                                                                                                                  cGAACACGTCAACCGCTTCGGTATCACCGAACAGGTTGTGCATGTTGCCGAGGATCTCCTGATATGCGCCGACCATAAAGAAACCGAGCATCGGCGGa                                                                                 <  1:461806/98‑1 (MQ=255)
                                                                                                                    aaCACGTCAACCGCTTCGGTATCACCGAACAGGGTGTGCATGTTGCCGAGGATCTCCTGATCTGCGCCGACCCTAAAGAAAGGGAGCAGCGGCGGATTCACTGGATCGTACTCCGGCATTGGCATTGTCGTGGCAATACCGTCACCATCa                           <  1:821086/150‑1 (MQ=255)
                                                                                                                     catgcacAACCTGTTCGGTGATACCGAAGCGGTTGTGCATGTTGCCGAGGATCTCCTGATATGCGCCGACCATAAAGAAACCGAGCATCGGCGGATTCTCTGGATCGTACTCCGGCATTGGCATTGTCGTGGCAATACCGTCACCATCaa                          >  2:42652/7‑150 (MQ=255)
                                                                                                                        cGTCAACCGCTTCGGTATCACCGAACAGGTTGTGCATGTTGCCGAGGATCTCCTGATATGCGCCGACCATAAAGAAACCGAGCATCGGCGGATTGTCTGGATCGTACTCCGGCATTGGCATTGTCGTGGCAATACCGTCACCATCAatat                       >  1:392256/1‑150 (MQ=255)
                                                                                                                        cGTCAACCGCTTCGGTATCACCGAACAGGTTGTGCATGTTGCCGAGGATCTCCTGATATGCGCCGACCATAAAGAAACCGAGCATCGGCGGATTCTCTGGATCGTACTCCGGCATTGGCATTGTCGTGGCAATACCGTCACCATCAatat                       >  1:1243696/1‑150 (MQ=255)
                                                                                                                          tCAACCGCTTCGGTATCACCGAACAGGTTGTGCATGTTGCCGAGGATCTCCGGATATGCGCCGACCATAATGAAACCGAGCATCGGCGGATTCTCTGGATCGTACTCCGGCATTGGCATTGTCGTGGCAATACCGTCACCATCAATATAg                     >  1:874792/1‑150 (MQ=255)
                                                                                                                           cAACCGCTTCGGTATCACCGAACAGGTTGTGCATGTTGCCGAGGATCTCCTGATATGCGCCGACCATAAAGAAACCGAGCATCGGCGGATTCTCTGGATCGTACTCCGGCATTGGCATTGTCGTGGCAATACCGTCACCATCAATATAGt                    >  2:236324/1‑150 (MQ=255)
                                                                                                                           cAACCGCTTCGGTATCACCGAACAGGTTGTGCATGTTGCCGAGGATCTCCTGATATGCGCCGACCATAAAGAAACCGAGCATCGGCGGATTCTCTGGATCGTACTCCGGCATTGGCATTGTCGTGGCAATACCGTCACCATCAATATAGt                    <  1:581792/150‑1 (MQ=255)
                                                                                                                             aCCGCTTCGGTATCACCGAACAGGTTGTGCATGTTGCCGAGGATCTCCTGATATGCGCCGACCATAAAGAAACCGAGCATCGGCGGATTCTCTGGATCGTACTCCGGCATTGGCATTGTCGTGGCAATACCGTCACCATCAATATAGTgg                  >  1:634299/1‑150 (MQ=255)
                                                                                                                               cGCTTCGGTATCACCGAACAGGTTGTGCATGTTGCCGAGGATCTCCTGATATGCGCCGACCATAAAGAAACCGAGCATCGGCGGATTCTCTGGATCGTACTCCGGCATTGGCATTGTCGTGGCAATACCGTCACCATCAATATAGTGGTc                <  2:1190995/150‑1 (MQ=255)
                                                                                                                                 cTTCGGTATCACCGAACAGGTTGTGCATGTTGCCGAGGATCTCCTGATATGCGCCGACCATAAAGAAACCGAGCATCGGCGGATTCTCTGGATCGTACTCCGGCATTGGCATTGTCGTGGCAATACCGTCACCATCAATATAGTGGTCGa              >  1:355835/1‑150 (MQ=255)
                                                                                                                                   tCGGTATCACCGAACAGGTTGTGCATGTTGCCGAGGATCTCCTGATATGCGCCGACCATAAAGAAACCGAGCATCGGCGGATTCTCTGGATCGTACTCCGGCATTGGTATTGTCGTGGCAATACCGTCACCATCAATATAGTGGTCGATa            >  2:1207805/1‑150 (MQ=255)
                                                                                                                                     ggTATCACCGAACAGGTTGTGCATGTTGCCGAGGATCTCCTGATATGCGCCGACCATAAAGAAACCGAGCATCGGCGGATTCTCTGGATCGTACTCCGGCATTGGCATTGTCGTGGCAATACCGTCACCATCAATATAGTGGTCGATAGc          <  2:580613/150‑1 (MQ=255)
                                                                                                                                      gTATCACCGAACAGGTTGGGCATGTTGCCGAGGATCTCCTGATATGCGCCGACCATAAAGAAACCGAGCATCGGCGGATTCTCTGGATCGTACTCCGGCATTGGCATTGTCGTGGCAATACCGTCACCATCAATATAGTGGTCGATAGcc         >  1:677837/1‑149 (MQ=255)
                                                                                                                                         tCACCGAACAGGTTGTGCATGTTGCCGAGCATCTCCTGATATGCGCCGACCATAAAGAAACCGAGCATCGGCGGATTCTCTGGATCGTACTCCGGCATTGGCATTGGCGTGGCAATACCGTCACCATCAATATAGTGGTCGATAGCACCg      <  2:464006/150‑1 (MQ=255)
                                                                                                                                          cACCGAACAGGTTGTGCATGTTGCCGAGGATCTCCTGATATGCGCCGACCATAAAGAAACCGAGCATCGGCGGATTCTCTGGATCGTACTCCGGCATTGGCATTGTCGTGGCAATACCGTCACCATCAATATAGTGGTCGATAGCACCGt     >  2:401984/1‑150 (MQ=255)
                                                                                                                                            ccGAACAGGTTGTGCATGTTGCCGAGGATCTCCTGATATGCGCCGACCATAAAGAAACCGAGCATCGGCGGATTCTCTGGATCGTACTCCGGCATTGGCATTGTCGTGGCAATACCGTCACCATCAATATAGTGGTCGATAGCACCGTCa   >  2:153119/1‑150 (MQ=255)
                                                                                                                                             cGAACAGGTTGTGCATGTTGCCGAGGATCTCCTGATATGCGCCGACCATAAAGAAACCGAGCATCGGCGGATTCTCTGGATCGTACTCCGGCATTGGCATTGTCGTGGCAATACCGTCACCATCAATATAGTGGTCGATAGCACCGTCag  >  1:881020/1‑150 (MQ=255)
                                                                                                                                             cGAACAGGTTGTGCATGTTGCCGAGGATCTCCTGATATGCGCCGACCATAAAGAAACCGAGCATCGGCGGATTCTCTGGATCGTACTCCGGCATTGGCATTGTCGTGGCAATACCGTCACCATCAATATAGTGGTCGATAGCACCGTCag  >  2:1573238/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
GATCGCGGAACTGGGTTAACAGCGTTTTCGGATCGAGCTGTACATATTGCAGCATGTCCGCCACGGTATCGCCTTCGTCAGACAGTTCTACTTCTACGCTACCGTCAGGGAAGACGAACACGTCAACCGCTTCGGTATCACCGAACAGGTTGTGCATGTTGCCGAGGATCTCCTGATATGCGCCGACCATAAAGAAACCGAGCATCGGCGGATTCTCTGGATCGTACTCCGGCATTGGCATTGTCGTGGCAATACCGTCACCATCAATATAGTGGTCGATAGCACCGTCAG  >  NC_000913/3084033‑3084323

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