Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,136,632 C→A 6.1% 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.1% 526.3 / 17.1 172intergenic (‑261/+31)pitB/gssphosphate transporter/glutathionylspermidine amidase and glutathionylspermidine synthetase
Reads supporting (aligned to +/- strand):  ref base C (87/74);  new base A (3/8);  total (90/82)
Fisher's exact test for biased strand distribution p-value = 1.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.87e-01

Only 100 of 180 total aligned reads displayed.

TGTATCAGACAATAAGCATTCATTCAAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCACCGTAGTTGCCGCCAACGGTGAAGGTACATACCTG  >  NC_000913/3136483‑3136779
                                                                                                                                                     |                                                                                                                                                   
tgtATCAGACAATAAGGATTCATTCAAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTAttcttc                                                                                                                                                     <  1:822776/150‑1 (MQ=255)
tgtATCAGACAATAAGCATTCATTCAAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTAttcttc                                                                                                                                                     >  2:672288/1‑150 (MQ=255)
   aTCAGACAATAAGCATTCATTCAAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTATTAATa                                                                                                                                                  >  1:1018414/1‑150 (MQ=255)
   aTCAGACAATAAGCATTCATTCAAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCAAATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTATTCATa                                                                                                                                                  >  2:1126857/1‑150 (MQ=255)
      aGACAATAAGCATTCATTCAAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATc                                                                                                                                               >  1:751564/1‑150 (MQ=255)
      aGACAATAAGCATTCATTCAAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATc                                                                                                                                               <  1:1297832/150‑1 (MQ=255)
          aaTAAGCATTCATTCAAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTa                                                                                                                                           <  1:1347175/150‑1 (MQ=255)
          aaTAAGCATTAATTCAAAACCTTTATAAATAATCAAAACACATGTATTACACTTGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCTTACTGATTATGAAGAATAATCAGTa                                                                                                                                           <  1:426211/150‑1 (MQ=255)
           aTAAGCATTCATTCAAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATGAAGAATAATCAGTAc                                                                                                                                          <  2:1508653/150‑1 (MQ=255)
             aaGCATTCATTCAAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACgg                                                                                                                                        <  1:487673/150‑1 (MQ=255)
              aGCATTCATTCAAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACggg                                                                                                                                       >  2:583576/1‑150 (MQ=255)
              aGCATTCATTCAAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATGAAGAATAATCAGTACggg                                                                                                                                       <  2:873979/150‑1 (MQ=255)
               tcattcattCAAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATGAAGAATAATCAGTACgggg                                                                                                                                      <  1:309235/149‑1 (MQ=255)
               gCATTCATTCAAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATGAAGAATAATCAGTACgggg                                                                                                                                      <  2:1713118/150‑1 (MQ=255)
                    cattcaAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCt                                                                                                                                 >  1:869566/1‑150 (MQ=255)
                      ttcaAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTac                                                                                                                               >  2:41991/1‑150 (MQ=255)
                      ttcaAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTac                                                                                                                               >  2:532905/1‑150 (MQ=255)
                      ttcaAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATGAAGAATAATCAGTACggggattaatc                                                                                                                               <  1:601733/150‑8 (MQ=38)
                        caAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTacaa                                                                                                                             <  1:423991/150‑1 (MQ=255)
                         aaaaCCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAAt                                                                                                                                                >  1:329064/1‑130 (MQ=255)
                          aaaCCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGTGGTTGCTacaaca                                                                                                                           >  1:669773/1‑150 (MQ=255)
                          aaaCCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTacaaca                                                                                                                           >  1:1004720/1‑150 (MQ=255)
                            aCCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATg                                                                                                                         <  2:1386762/150‑1 (MQ=255)
                            aCCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATg                                                                                                                         <  2:1201110/150‑1 (MQ=255)
                                ttAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTAtt                                                                                                                     <  1:856833/150‑1 (MQ=255)
                                  aGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAATACTGTATTTATACAACATGTATTAc                                                                                                                   >  1:100927/1‑150 (MQ=255)
                                       taatCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACttttt                                                                                                              <  1:271267/150‑1 (MQ=255)
                                       taatCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACttttt                                                                                                              <  1:1487636/150‑1 (MQ=255)
                                       taatCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACttttt                                                                                                              <  2:857791/150‑1 (MQ=255)
                                         atCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTca                                                                                                            >  2:738060/1‑150 (MQ=255)
                                         atCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTca                                                                                                            <  2:806073/150‑1 (MQ=255)
                                         atCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTca                                                                                                            >  1:975040/1‑150 (MQ=255)
                                         atCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTca                                                                                                            <  1:387543/150‑1 (MQ=255)
                                            aaaaCACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTTATTAATCCCCTTACTGATTATTATTCATAATCAGTACGGGGTTGCTACAACATTTATTAATTTTTcacaa                                                                                                         >  1:1234041/1‑148 (MQ=255)
                                             aaaCACAGGTATTCCACTGGTTAATATATTCACATATGTAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTcaccac                                                                                                        >  1:616772/1‑150 (MQ=255)
                                             aaaCACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTcaccac                                                                                                        <  1:1276421/150‑1 (MQ=255)
                                              aaCACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACcaca                                                                                                       <  1:720380/150‑1 (MQ=255)
                                                cacaGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAAt                                                                                                     >  1:536426/1‑150 (MQ=255)
                                                      tATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACgg                                                                                               <  2:697405/150‑1 (MQ=255)
                                                          ccACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCa                                                                                           >  2:202519/1‑150 (MQ=255)
                                                          ccACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCa                                                                                           >  1:1676714/1‑150 (MQ=255)
                                                            aCTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAat                                                                                         <  2:751564/150‑1 (MQ=255)
                                                             cTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAata                                                                                        <  1:1344134/150‑1 (MQ=255)
                                                              tGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAatat                                                                                       <  2:1004720/150‑1 (MQ=255)
                                                               ggTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATc                                                                                      >  2:1467805/1‑150 (MQ=255)
                                                               cgTTAATATATTCACATATGAAAAGAATAAATATTGGAATATATAAATATTCAATATTTTGATTAAACCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAAAATGAATTACTTTTTCACCACAATTAACTGTTCAaaata                                                                                      >  2:569633/2‑149 (MQ=255)
                                                                gTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCa                                                                                     <  1:1552957/150‑1 (MQ=255)
                                                                 ttAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATTTCAc                                                                                    >  1:644619/1‑150 (MQ=255)
                                                                 ttAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCAc                                                                                    >  2:894307/1‑150 (MQ=255)
                                                                       tatTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCtt                                                                              >  2:1338212/1‑150 (MQ=255)
                                                                           cacaTATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCtttttt                                                                          >  2:1578608/1‑150 (MQ=255)
                                                                           aacaTATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCtttttt                                                                          <  1:1178671/149‑1 (MQ=255)
                                                                             caTATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTtga                                                                        >  2:203778/1‑150 (MQ=255)
                                                                              atatGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTtgat                                                                       <  1:814461/150‑1 (MQ=255)
                                                                               tatGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTtgatg                                                                      >  1:62078/1‑150 (MQ=255)
                                                                               tatGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTtgatg                                                                      <  1:313718/150‑1 (MQ=255)
                                                                                     atgaaTAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGt                                                                <  1:738060/150‑1 (MQ=255)
                                                                                       gaaTAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGa                                                              >  2:810778/1‑150 (MQ=255)
                                                                                        aataaataTTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGAt                                                             >  1:1539020/1‑150 (MQ=255)
                                                                                        aataaataTTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGAt                                                             <  1:237724/150‑1 (MQ=255)
                                                                                        aataaataTTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGAt                                                             <  1:412060/150‑1 (MQ=255)
                                                                                           aaataTTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATtca                                                          <  1:1589916/150‑1 (MQ=255)
                                                                                             atatTGGAATATATAAATATTGAATATTTTGATTAATCCCAGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATtcatt                                                        <  1:179702/150‑2 (MQ=255)
                                                                                              tatTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACAAGTGATtcatca                                                       >  1:1001619/1‑150 (MQ=255)
                                                                                                ttGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCAcc                                                     <  1:1400245/150‑1 (MQ=255)
                                                                                                ttGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCAcc                                                     >  1:979366/1‑150 (MQ=255)
                                                                                                 tGGAATATATAAATATTGAATATTTTGTTTAATACCAGTACTTATTATTCTTCATAATCAGTACGGGGTTGCTAAAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATTACCAGTGATTCATCACcg                                                    <  2:591108/150‑1 (MQ=255)
                                                                                                 tGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACcg                                                    >  2:1205978/1‑150 (MQ=255)
                                                                                                   gAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACcgc                                                   >  2:1116876/1‑149 (MQ=255)
                                                                                                   gAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACcgc                                                   >  2:254144/1‑149 (MQ=255)
                                                                                                   gAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACcgc                                                   <  1:1116876/149‑1 (MQ=255)
                                                                                                   gAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAAATAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACcgcg                                                  >  2:1219703/1‑150 (MQ=255)
                                                                                                    aaTATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGTGGTTGCTACAAAATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCgcgc                                                 >  2:325589/1‑150 (MQ=255)
                                                                                                    aaTATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCgtgc                                                 >  2:1421109/1‑150 (MQ=255)
                                                                                                       atataAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCaaa                                              <  1:1615651/150‑1 (MQ=255)
                                                                                                        tataAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAac                                             >  2:779916/1‑150 (MQ=255)
                                                                                                         ataAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAaca                                            <  2:1676714/150‑1 (MQ=255)
                                                                                                         ataAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAaca                                            <  2:1234041/150‑1 (MQ=255)
                                                                                                         ataAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAaca                                            <  2:1226486/150‑1 (MQ=255)
                                                                                                         ataAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAaca                                            >  2:846608/1‑150 (MQ=255)
                                                                                                              tatTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTcc                                       <  1:1685383/150‑1 (MQ=255)
                                                                                                               atTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTccc                                      >  2:185998/1‑150 (MQ=255)
                                                                                                               atTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTccc                                      >  1:634814/1‑150 (MQ=255)
                                                                                                                  gAATATTTTGATTAATCCCCGTACTGATTATTCTTCATTATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCTCCGCGCAAACACGTCCCAcc                                   >  2:1566412/1‑150 (MQ=255)
                                                                                                                   aaTATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCACCg                                  <  2:965967/150‑1 (MQ=255)
                                                                                                                          tGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCACCGTAGTTgc                           >  2:1690013/1‑150 (MQ=255)
                                                                                                                            aTTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCCCCACCATTAACGGTTCATTATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCACCGTAGTTgccg                         >  1:1496069/1‑150 (MQ=255)
                                                                                                                              tAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCACCGTAGTTgccgcc                       <  2:62078/150‑1 (MQ=255)
                                                                                                                               aaTCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCACCGTAGTTGCCGCCa                      >  2:391909/1‑150 (MQ=255)
                                                                                                                               aaTCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCACCGTAGTTGCCGCCa                      >  1:999450/1‑150 (MQ=255)
                                                                                                                                  ccccGTACTGATTATGAAGAATAATCAGTACGGGGTTGCTTCAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCAGTCTCTTTGTTGATGCCCAGTGAGTCATCACCGCGCAAACCCGGCCCGCCGTAGTTGCCGGAAACg                   >  2:346451/1‑150 (MQ=255)
                                                                                                                                   cccGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGTCCAGTGATTCATCACCGCGCAAACACGTCCCACCGTAGTTGCCGCCAACgg                  <  1:846608/150‑1 (MQ=255)
                                                                                                                                        aCTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCACCGTAGTTGCCGCCAACGGTGAAg             >  1:906694/1‑150 (MQ=255)
                                                                                                                                         cTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCACCGTAGTTGCCGCCAACGGTGAAgg            <  1:1578608/150‑1 (MQ=255)
                                                                                                                                          tGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCACCGTAGTTGCCGCCAACGGTGAAGGt           >  1:1096970/1‑150 (MQ=255)
                                                                                                                                           gATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCACCGTAGTTGCCGCCAACGGTGAAGGTa          >  2:354781/1‑150 (MQ=255)
                                                                                                                                            attattCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCACCGTAGTTGCCGCCAACGGTGAAGGTAc         <  2:1316378/150‑1 (MQ=255)
                                                                                                                                               attCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCACCGTAGTTGCCGCCAACGGTGAAGGTACATa      <  2:1023346/150‑1 (MQ=255)
                                                                                                                                                 tcttcATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCACCGTAGTTGCCGCCAACGGTGAAGGTACATAcc    <  2:1537853/150‑1 (MQ=255)
                                                                                                                                                 tcttcATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCACCGTAGTTGCCGCCAACGGTGAAGGTACATAcc    >  1:964229/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TGTATCAGACAATAAGCATTCATTCAAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCACCGTAGTTGCCGCCAACGGTGAAGGTACATACCTG  >  NC_000913/3136483‑3136779

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