Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,136,632 C→A 6.7% intergenic (‑261/+31) pitB ← / ← gss phosphate transporter/glutathionylspermidine amidase and glutathionylspermidine synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,136,6320CA 6.7% 212.8 / 6.4 83intergenic (‑261/+31)pitB/gssphosphate transporter/glutathionylspermidine amidase and glutathionylspermidine synthetase
Reads supporting (aligned to +/- strand):  ref base C (38/38);  new base A (2/4);  total (41/42)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.24e-01

GTATCAGACAATAAGCATTCATTCAAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCACCGTAGTTGCCGCCAACGGTGAAGGTACATACCTG  >  NC_000913/3136484‑3136779
                                                                                                                                                    |                                                                                                                                                   
gtATCAGACAATAAGCATTCATTCAAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTATTCa                                                                                                                                                    >  2:1152181/1‑150 (MQ=255)
  aTCAGACAATAAGCATTCATTCAAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATa                                                                                                                                                  >  1:680722/1‑150 (MQ=255)
         aaTAAGCATTCATTCAAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTa                                                                                                                                           <  1:1298224/150‑1 (MQ=255)
          aTAAGCATTCATTCAAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTAc                                                                                                                                          >  2:269332/1‑150 (MQ=255)
           tAAGCATTCATTCAAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATGAAGAATAATCAGTACg                                                                                                                                         <  1:433814/150‑1 (MQ=255)
           tAAGCATTCATTCAAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGCAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACg                                                                                                                                         >  1:400044/1‑150 (MQ=255)
           tAAGCATTCATTCAAAACCTTTAGAAATAATCAAAACACAGGTATTACACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACg                                                                                                                                         <  1:36922/150‑1 (MQ=255)
           aaaGCATTCATTCAAAACCTTTAGAAATAATCAAAACACAGGTAATCCACTTGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACg                                                                                                                                         <  2:473417/149‑1 (MQ=255)
            aaGCATTCATTCAAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACgg                                                                                                                                        >  2:766469/1‑150 (MQ=255)
              gCATTCATTCAAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACgggg                                                                                                                                      <  2:248083/150‑1 (MQ=255)
              gCATTCATTCAAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATGAAGAATAATCAGTACgggg                                                                                                                                      <  1:484826/150‑1 (MQ=255)
              gCATTCATTCAAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATGAAGAATAATCAGTACgggg                                                                                                                                      <  2:1196177/150‑1 (MQ=255)
                attcattcAAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATGAAGAATAATCAGTACggggat                                                                                                                                    <  1:167424/150‑3 (MQ=39)
                actcattcaAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGtt                                                                                                                                    >  1:1237082/3‑150 (MQ=255)
                   cattcaAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCt                                                                                                                                 >  2:1110649/1‑150 (MQ=255)
                             cTTTAGAAATAATCAAAACACATGTATTCCAATTGTTTATATATTAACAAATGAAATGAATAACTATTGGAATATATAAATATTGAATCTTTAGATTAATCCCCGTACGGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTa                                                                                                                       <  1:191352/150‑1 (MQ=255)
                             cTTTAGAAATAATCAAAACACAGGTATTCCACTGGGTAATATAATCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTa                                                                                                                       <  2:882912/150‑1 (MQ=255)
                             atttAGAAATAATCAAAACACAGGTATGCCACTGGTTAATATAATCACATATGAAATGAATAAATATTGGAATCTATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTa                                                                                                                       <  2:802213/149‑1 (MQ=255)
                               ttAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTAtt                                                                                                                     <  1:1058631/150‑1 (MQ=255)
                                tAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTa                                                                                                                    <  2:1101393/150‑1 (MQ=255)
                                     ataatCAAAACACAGGTATTCCACTGGTTAATATAATCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACtttt                                                                                                               <  2:604420/150‑1 (MQ=255)
                                       aatCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTc                                                                                                             <  2:620498/150‑1 (MQ=255)
                                         tCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGGAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTcac                                                                                                           >  1:1258306/1‑150 (MQ=255)
                                                 caGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACTGGGTTGCTACAACATGTATTACTTTTTCACAACAATTa                                                                                                   >  2:91066/1‑150 (MQ=255)
                                                   ggTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACTGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAAc                                                                                                 >  1:157481/1‑150 (MQ=255)
                                                   ggTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTACATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAAc                                                                                                 >  2:1300725/1‑150 (MQ=255)
                                                      aTTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGt                                                                                              <  1:896826/150‑1 (MQ=255)
                                                        tCCACTGGTTAATATATTCACATATGAAATGAATATATATTGGAATATATAAATATTGAATATTTTGATTAATCACCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGAATTACTTTTTCACCACAATTAACGGtta                                                                                            >  2:641860/1‑149 (MQ=255)
                                                        tCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGTATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTc                                                                                            <  1:195296/150‑1 (MQ=255)
                                                        tCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTc                                                                                            <  1:612471/150‑1 (MQ=255)
                                                        tCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTc                                                                                            >  2:309843/1‑150 (MQ=255)
                                                         ccACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCa                                                                                           >  2:1192902/1‑150 (MQ=255)
                                                          cACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCaa                                                                                          <  1:1139132/150‑1 (MQ=255)
                                                                ttAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTTAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCTATATCAc                                                                                    <  2:1099046/150‑1 (MQ=255)
                                                                ttAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCAc                                                                                    >  1:509396/1‑150 (MQ=255)
                                                                 tAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCAct                                                                                   >  2:913971/1‑150 (MQ=255)
                                                                  aaTATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCActc                                                                                  <  1:7611/150‑1 (MQ=255)
                                                                         tCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCttttt                                                                           >  2:213288/1‑150 (MQ=255)
                                                                         tCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCtttat                                                                           >  2:1273567/1‑148 (MQ=255)
                                                                          cacaTATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCtttttt                                                                          >  2:518871/1‑150 (MQ=255)
                                                                            caTATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCaa                                                                                          <  2:782959/132‑1 (MQ=255)
                                                                            caTATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCaa                                                                                          >  1:782959/1‑132 (MQ=255)
                                                                            caTATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTtga                                                                        >  2:230458/1‑150 (MQ=255)
                                                                             atatGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTtgat                                                                       <  1:913971/150‑1 (MQ=255)
                                                                              tatGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCACTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCCCTCTCTTTTTtgatg                                                                      >  1:933550/1‑150 (MQ=255)
                                                                              tatGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTAATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTtgatt                                                                      >  1:831196/1‑149 (MQ=255)
                                                                                tgaaatgaaTAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGAc                                                                    <  2:360106/150‑1 (MQ=255)
                                                                                     tgaaTAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATTACCAGTg                                                               >  2:922122/1‑150 (MQ=255)
                                                                                            ttatTGTAATATATAAATATTGAATATTTTGATTAATCCCCGTACTTATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATtcatc                                                        <  1:831654/149‑1 (MQ=255)
                                                                                                ttgAATATATAAATATTGAATATTTTGATTAATGCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCTCCACAATTAACTGTTCAATATCACTCTCTTTTTTGATGACCATTGATTCATCACcg                                                    <  2:296684/148‑1 (MQ=255)
                                                                                                tGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACcg                                                    <  2:485548/150‑1 (MQ=255)
                                                                                                     tatataAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTGAGAATCAGTACGGGGTTGCTACAACATGAATTAGTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTGGATGACCAGTGAGTCATGACCTCGCaa                                               >  2:187779/1‑150 (MQ=255)
                                                                                                     tatataAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACTACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCaa                                               <  2:163058/150‑1 (MQ=255)
                                                                                                       tataAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAac                                             <  1:927545/150‑1 (MQ=255)
                                                                                                       tataAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAac                                             <  1:676906/150‑1 (MQ=255)
                                                                                                       tataAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCAATCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAac                                             >  1:772152/1‑150 (MQ=255)
                                                                                                           aaTATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATAACAGCGCAAACACGt                                         >  1:373406/1‑150 (MQ=255)
                                                                                                             tatTGAATATTTTTATTAATCCCCCTACTGATTATTCTTCATAATCAGTACGGGGTAGCTACAACATGTATTACTTTTTCACCACACTTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTcc                                       <  1:214286/150‑1 (MQ=255)
                                                                                                               ttGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTcccc                                     >  2:735951/1‑149 (MQ=255)
                                                                                                                tGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACgg                                                                                               >  1:941554/1‑91 (MQ=255)
                                                                                                                tGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACgg                                                                                               <  2:941554/91‑1 (MQ=255)
                                                                                                                tGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCAc                                    >  1:847114/1‑150 (MQ=255)
                                                                                                                     atTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCACCGTAg                               <  2:157481/150‑1 (MQ=255)
                                                                                                                        ttGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCACCGTAGTTg                            >  2:464742/1‑150 (MQ=255)
                                                                                                                         tGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGTGGTTTCTACAAAATGTATTACTTTTTCACCACTATTAACGGTGCAATATCACTCTCTTTTTTGATGACCAGCGATGCATCACCGCGCAAACACGTCCCACCGTAGTTgc                           <  1:857493/150‑1 (MQ=255)
                                                                                                                         tGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCACCGTAGTTgc                           >  1:695593/1‑150 (MQ=255)
                                                                                                                         tGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCACCGTAGTTgc                           >  1:1237267/1‑150 (MQ=255)
                                                                                                                         tGATTAATCCCCGTACTGATTATTCTTCATAATAAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATTACCATTGATTCATCACCGCGCAAACACGTCCCACCGTAGTTgc                           <  1:641293/150‑1 (MQ=255)
                                                                                                                            ttAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCACCGTAGTTgccgc                        <  1:404429/150‑1 (MQ=255)
                                                                                                                            agcaaccccGTACTGATTATGAAGAATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATAACCGCGCAAACACGTCCCACCGTAGTTgccgc                        >  2:277378/6‑150 (MQ=38)
                                                                                                                                    cGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCTCCGTAGTTGCCGCCAACGGTg                >  1:1694/1‑150 (MQ=255)
                                                                                                                                    cGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCACCGTAGTTGCCGCCAACGGTg                <  1:766469/150‑1 (MQ=255)
                                                                                                                                    cGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGGTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCACCGTAGTTGCCGCCAACGGTg                <  2:1056862/150‑1 (MQ=255)
                                                                                                                                       aCTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCACCGTAGTTGCCGCCAACGGTGAAg             >  2:380128/1‑150 (MQ=255)
                                                                                                                                        cTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCACCGTAGTTGCCGCCAACGGTGAAgg            >  2:437371/1‑150 (MQ=255)
                                                                                                                                        atGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCACCGTAGTTGCCGCCAACGGTGAAgg            <  2:247471/149‑1 (MQ=255)
                                                                                                                                         tGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCACCGTAGTTGCCGCCAACGGTGAAGGt           <  1:882762/150‑1 (MQ=255)
                                                                                                                                          gATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAACCACGTCCCACCGTAGTTGCCGCCAACGGTGAAGGTa          >  1:366513/1‑150 (MQ=255)
                                                                                                                                          gATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCACCGTAGTTGCCGCCAACGGTGAAGGTa          <  2:499679/150‑1 (MQ=255)
                                                                                                                                          gATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCAATCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCACCGTAGTTGCCGCCAACGGTGAAGGTa          >  1:1059978/1‑150 (MQ=255)
                                                                                                                                           attattCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCACCGTAGTTGCCGCCAACGGTGAAGGTAc         >  2:1061991/1‑150 (MQ=255)
                                                                                                                                           attattCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCACCGTAGTTGCCGCCAACGGTGAAGGTAc         <  1:574445/150‑1 (MQ=255)
                                                                                                                                             tattCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAAATAACGGTTCAATATCACTGTCTTTTTTGATGACAAGTGATTCATCACCGCGCTACCGCGTCCCACCGTAGTTGCCGCCTACGGTGCAGGTACAt       >  2:560855/1‑150 (MQ=255)
                                                                                                                                                 cttcATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTAGATGACCAGTGATTCATCACAGCGCAAACACGTCCCACCGTAGTTGCCGCCAAAGGTGAAGGTACATACCt   >  1:1297326/1‑150 (MQ=255)
                                                                                                                                                  ttcATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCACCGTAGTTGCCGCCAACGGTGAAGGTACATACCTg  <  2:795165/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GTATCAGACAATAAGCATTCATTCAAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCACCGTAGTTGCCGCCAACGGTGAAGGTACATACCTG  >  NC_000913/3136484‑3136779

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: