Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,136,632 C→A 6.3% 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.3% 442.0 / 13.6 159intergenic (‑261/+31)pitB/gssphosphate transporter/glutathionylspermidine amidase and glutathionylspermidine synthetase
Reads supporting (aligned to +/- strand):  ref base C (81/66);  new base A (7/4);  total (88/71)
Fisher's exact test for biased strand distribution p-value = 7.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.20e-01

Only 100 of 173 total aligned reads displayed.

TGTATCAGACAATAAGCATTCATTCAAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTAC‑TTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCACCGTAGTTGCCGCCAACGGTGAAGGTACATACCTGA  >  NC_000913/3136483‑3136780
                                                                                                                                                     |                                                                                                                                                     
tgtATCAGACAATAAGCATTCATTCAAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTAttcttc                                                                                                                                                       <  2:427516/150‑1 (MQ=255)
tgtATCAGACAATAAGCATTCATTCAAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTAttcttc                                                                                                                                                       <  1:284259/150‑1 (MQ=255)
tgtATCAGACAATAAGCATTCATTCAAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTAttcttc                                                                                                                                                       <  2:245824/150‑1 (MQ=255)
 gtATCAGACAATAATCATTCATTCAAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCa                                                                                                                                                      <  1:609405/150‑1 (MQ=255)
 gtATCAGACAATAAGCATTCATTCAAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCa                                                                                                                                                      <  1:1192548/150‑1 (MQ=255)
 gtATCAGACAATAAGCATTCATTCAAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCa                                                                                                                                                      <  2:527823/150‑1 (MQ=255)
  tATCAGACAATAAGCATTCATTCAAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCAt                                                                                                                                                     <  2:1290344/150‑1 (MQ=255)
   aTCAGACAATAAGCATTCATTCAAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATa                                                                                                                                                    >  1:1520897/1‑150 (MQ=255)
      aGACAATAAGCATTCATTCAAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATc                                                                                                                                                 >  2:156681/1‑150 (MQ=255)
       gACAATAAGCATTCATTCAAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCa                                                                                                                                                >  2:1717742/1‑150 (MQ=255)
        aCAATAAGCATTCATTCAAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAg                                                                                                                                               <  1:557637/150‑1 (MQ=255)
          aaTAAGCATTCATTCAAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTa                                                                                                                                             <  2:1152362/150‑1 (MQ=255)
            tAAGCATTCATTCAAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACg                                                                                                                                           >  1:1388544/1‑150 (MQ=255)
              aGCATTCATTCAAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACggg                                                                                                                                         <  2:27824/150‑1 (MQ=255)
              aGCATTCATTCAAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACggg                                                                                                                                         >  1:645658/1‑150 (MQ=255)
              aGCATTCATTCAAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACggg                                                                                                                                         >  1:1013440/1‑150 (MQ=255)
              aGCATTCATTCAAAACCTTTAGAAAAAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATGAAGAATAATCAGTACggg                                                                                                                                         <  2:347441/150‑1 (MQ=255)
               gCATTCATTCAAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATGAAGAATAATCAGTACgggg                                                                                                                                        <  2:1541383/150‑1 (MQ=255)
                 attcattcAAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGCAATATATAAATATTGAATATTTTGATTAATCCCCTTACTGATTATTCTTCATAATCAGTACGTGGtt                                                                                                                                      >  2:669406/1‑150 (MQ=255)
                 attcattaAAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATGAAGAATAATCAGTACggggat                                                                                                                                      <  2:125480/150‑3 (MQ=39)
                  ttcattcaAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATGAAGAATAATCAGTACggggatt                                                                                                                                     <  2:764561/150‑4 (MQ=39)
                   tcattcaAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGc                                                                                                                                    <  2:1520897/150‑1 (MQ=255)
                    cattcaAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGTGTTGCt                                                                                                                                   >  1:1541996/1‑150 (MQ=255)
                    aattcaAAACCTTTAGAAATAATCAAAACACATTTATTCCACTGTTTAATATATTCACATATGAAATGCATAAATATTGGAATATATAATTATGGTCTATTTTGATTAATCCCCGTACTGATTAGTCTTCATAATCAGTACGGGGTTTCt                                                                                                                                   <  1:280474/149‑1 (MQ=255)
                    aattcaAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCt                                                                                                                                   <  2:793198/149‑1 (MQ=255)
                     attcaAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCACCGTACTGATTATTCTTAATAATCAGTACGGGGTTGCTa                                                                                                                                  >  1:575779/1‑150 (MQ=255)
                       tcaAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGTGTTGCTaca                                                                                                                                >  1:261126/1‑150 (MQ=255)
                         ataaCCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTacaac                                                                                                                              <  2:232980/148‑1 (MQ=255)
                         aaaaCCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTacaac                                                                                                                              <  1:1526569/150‑1 (MQ=255)
                         aaaaCCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACAAATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTacaac                                                                                                                              <  2:312630/150‑1 (MQ=255)
                          aaaCCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTacaaca                                                                                                                             <  1:1472419/150‑1 (MQ=255)
                           aaCCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACAt                                                                                                                            <  2:1549817/150‑1 (MQ=255)
                               tttAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCTGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTAt                                                                                                                        >  2:1271898/1‑150 (MQ=255)
                                   gAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTAC‑t                                                                                                                   >  2:1177679/1‑150 (MQ=255)
                                      ataatCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTAC‑tttt                                                                                                                <  2:1688885/150‑1 (MQ=255)
                                         atCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTAC‑TTTTTca                                                                                                             >  2:209052/1‑150 (MQ=255)
                                         atCAAAACACACTTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATAATGGAATATATAAATATTGAATATTTTGATTAAACCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTAC‑TTTTTca                                                                                                             <  1:413178/150‑1 (MQ=255)
                                           cAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTAC‑TTTTTcacc                                                                                                           <  2:538672/150‑1 (MQ=255)
                                            aaaaCACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTAC‑TTTTTcacca                                                                                                          >  2:1130399/1‑150 (MQ=255)
                                            aaaaCACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTAC‑TTTTTcacca                                                                                                          >  1:303913/1‑150 (MQ=255)
                                             aaaCACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTAC‑TTTTTcaccac                                                                                                         >  2:1582308/1‑150 (MQ=255)
                                                    ggTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTAC‑TTTTTCACCACAATTAAc                                                                                                  <  1:993086/150‑1 (MQ=255)
                                                     gTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTAC‑TTTTTCACCACAATTAACg                                                                                                 <  1:465742/150‑1 (MQ=255)
                                                          ccACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTAC‑TTTTTCACCACAATTAACGGTTCa                                                                                            >  2:1622474/1‑150 (MQ=255)
                                                          ccACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTAC‑TTTTTCACCACAATTAACGGTTCa                                                                                            >  2:808559/1‑150 (MQ=255)
                                                          ccACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTAC‑TTTTTCACCACAATTAACGGTTCa                                                                                            >  2:481232/1‑150 (MQ=255)
                                                            aCTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTAC‑TTTTTCACCACAATTAACGGTTCAat                                                                                          >  2:1424051/1‑150 (MQ=255)
                                                            aCTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTAC‑TTTTTCACCACAATTAACGGTTCAat                                                                                          <  1:859535/150‑1 (MQ=255)
                                                             cTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTAC‑TTTTTCACCACAATTAACGGTTCAata                                                                                         <  1:481909/150‑1 (MQ=255)
                                                              tGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTAAGGTGTTGCTACAACATTTATTAC‑TTTTTCACAAAAATTAACGGTTCAatat                                                                                        >  1:317614/1‑150 (MQ=255)
                                                                gTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTAC‑TTTTTCACCACAATTAACGGTTCAATATCa                                                                                      <  1:597455/150‑1 (MQ=255)
                                                                 ttAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTAC‑TTTTTCACCACAATTAACGGTTCAATATCAc                                                                                     >  1:611652/1‑150 (MQ=255)
                                                                   aaTATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTAATAATAAGTACTGGGTAGCTACAACATGTATTAC‑TTTTTCACCACAATTAACGTTTCAATATCaata                                                                                   >  1:102793/1‑147 (MQ=255)
                                                                     tatatTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTAC‑TTTTTCACCACAATTAACGGTTCAATATCActctc                                                                                 <  1:38041/150‑1 (MQ=255)
                                                                      atatTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACTGGTTTGCTAAAACATGTATTAC‑TTTTTCACCACAATTAACTGTTCAATATAActctct                                                                                >  1:101729/1‑150 (MQ=255)
                                                                      atatTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCAGTACTGATTATTCTTCATAATCAGTACTGGGTTGCTACAAAATGTATTAA‑TTTTTCAACACAAATAACGGTTCAATATCActctct                                                                                >  1:106266/1‑150 (MQ=255)
                                                                       tatTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTAC‑TTTTTCACCACAATTAACGGTTCAATATCACTCTCtt                                                                               >  2:800321/1‑150 (MQ=255)
                                                                            acaTATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTAC‑TTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTtg                                                                          <  2:1496436/150‑1 (MQ=255)
                                                                              atatGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTAC‑TTTTTCACCACAATTAACGGTTCAATATCACTCTCttttt                                                                            >  1:719228/1‑146 (MQ=255)
                                                                              atatGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTAC‑TTTTTCACCACAATTAACGGTTCAATATCACTCTCttttt                                                                            <  2:719228/146‑1 (MQ=255)
                                                                              atatGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTAC‑TTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTtgat                                                                        <  1:1417986/150‑1 (MQ=255)
                                                                               tatGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCTGTACGGGGTTGCTACAACATGTATTAC‑TTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTtgatg                                                                       >  1:1141160/1‑150 (MQ=255)
                                                                                 tgaaatgaaTAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTAC‑TTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGAc                                                                     <  2:1401590/150‑1 (MQ=255)
                                                                                  gaaatgaaTAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTAC‑TTTTTCACCACAATTAACGGTTCAATATCACTCTCTATTTTGATGAca                                                                    >  1:705069/1‑149 (MQ=255)
                                                                                      tgaaTAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTAC‑TTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTg                                                                >  2:583955/1‑150 (MQ=255)
                                                                                        aataaataTTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTAC‑TTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGAt                                                              >  2:1380474/1‑150 (MQ=255)
                                                                                        aataaataTTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTAC‑TTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGAt                                                              <  2:303913/150‑1 (MQ=255)
                                                                                         ataaataTTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTAC‑TTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGAtt                                                             >  1:1563989/1‑150 (MQ=255)
                                                                                         ataaataTTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTAC‑TTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGAtt                                                             <  1:209052/150‑1 (MQ=255)
                                                                                                ttGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTAC‑TTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCAcc                                                      >  2:1345913/1‑150 (MQ=255)
                                                                                                 tGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTAC‑TTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACcg                                                     <  2:1121034/150‑1 (MQ=255)
                                                                                                     atatatAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTAC‑TTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCa                                                 <  1:177627/150‑1 (MQ=255)
                                                                                                                  gAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTAC‑TTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCAcc                                    >  2:1714542/1‑150 (MQ=255)
                                                                                                                     tatTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTAC‑TTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCACCGTa                                 <  1:537934/150‑1 (MQ=255)
                                                                                                                       ttttGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTAC‑TTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATAACCGCGCAAACACGTCCCACCGTAGt                               >  2:1342133/1‑150 (MQ=255)
                                                                                                                         ttGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTAC‑TTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCACCGTAGTTg                             >  2:400642/1‑150 (MQ=255)
                                                                                                                         ttGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTAC‑TTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACAAGTGATTCATCACCGCGCAAACACGTCCCACCGTAGTTg                             >  1:1529021/1‑150 (MQ=255)
                                                                                                                          tGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTAC‑TTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCACCGTAGTAgc                            >  2:427639/1‑150 (MQ=255)
                                                                                                                            aTTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTAC‑TTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCACCGTAGTTgccg                          >  1:661156/1‑150 (MQ=255)
                                                                                                                            aTTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTAC‑TTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCACCGTAGTTgccg                          >  2:1393178/1‑150 (MQ=255)
                                                                                                                               aaTCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTAC‑TTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCACCGTAGTTGCCGCCa                       >  2:1078440/1‑150 (MQ=255)
                                                                                                                                aaccccGTACTGATTATGAAGAATAATCAGTACGGGGTTGCTACAACATGTATTAC‑TTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCACCGTAGTTGCCGCCaa                      >  1:1062647/3‑150 (MQ=39)
                                                                                                                                 accccGTACAGATTATTCTTCATAATCAGTAAGGGGTTGCTACAACATGTATTAC‑TTAATCACCACAATTACCGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCACCGTAGTTGCCGCCAAc                     <  1:1681960/149‑1 (MQ=255)
                                                                                                                                   cccGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTAC‑TTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCACCGTAGTTGCCGCCAACgg                   >  1:940353/1‑150 (MQ=255)
                                                                                                                                   cccGTACTGATTATTCTTCATAATCAGTACGGGGTTGCAACAACATGTATTAC‑TTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCACCGTAGTTGCCGCCAACgg                   >  2:656708/1‑150 (MQ=255)
                                                                                                                                   cccGTACTGATTATGAAGAATAATCAGTACGGGGTTGCTACAACATGTATTAC‑TTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCACCGTAGTTGCCGCCAACgg                   >  1:1122476/1‑150 (MQ=255)
                                                                                                                                     cGTACTGATTATGAAGAATAATCAGTACGGGGTTGCTACAACATGTATTAC‑TTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCACCGTAGTTGCCGCCAACGGTg                 >  1:1680813/1‑150 (MQ=255)
                                                                                                                                      gTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTAC‑TTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCACCGTAGTTGCCGCCAACGGTGa                >  2:1010170/1‑150 (MQ=255)
                                                                                                                                        aCTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTAC‑TTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCACCGTAGTTGCCGCCAACGGTGAAg              <  1:1412569/150‑1 (MQ=255)
                                                                                                                                        aCTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTAATAC‑TTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCACCGTAGTTGCCGCCAACGGTGAAg              <  1:500160/150‑1 (MQ=255)
                                                                                                                                         cTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTAC‑TTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCACCGTAGTTGCCGCCAACGGTGAAgg             <  1:274179/150‑1 (MQ=255)
                                                                                                                                           gATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTAC‑TTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCACCGTAGTTGCCGCCAACGGTGAAGGTa           >  1:428502/1‑150 (MQ=255)
                                                                                                                                            attattCTTCATAATCAGTACGGGGTTGCTACAACATGTATTACTTTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCACCGTAGTTGCCGCCAACGGTGAAGGTa           <  1:229693/150‑1 (MQ=255)
                                                                                                                                               attCTTCATAATCAGTACGGGGTTGCTACAACATGTATTAC‑TTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCACCGTAGTTGCCGCCAACGGTGAAGGTACATa       >  2:1529030/1‑150 (MQ=255)
                                                                                                                                               attCTTCATAATCAGTACGGGGTTGCTACAACATGTATTAC‑TTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCACCGTAGTTGCCGCCAACGGTGAAGGTACATa       >  1:512540/1‑150 (MQ=255)
                                                                                                                                                 tcttcATAATCAGTACTGGGTTGCTACAACATGTATTAC‑TTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCACCGTAGTTGCCGCCAACGGTGAAGGTACATAcc     <  2:267127/150‑1 (MQ=255)
                                                                                                                                                 tcttcATAATCAGTACGGGGTTGCTACAACATGTATTAC‑TTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCACCGTAGTTGCCGCCAACGGTGAAGGTACATAcc     <  2:645658/150‑1 (MQ=255)
                                                                                                                                                  cttcATAATCAGTACGGGGTTGCTACAACATGTATTAC‑TTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCACCGTAGTTGCCGCCAACGGTGAAGGTACATACCt    >  2:970223/1‑150 (MQ=255)
                                                                                                                                                  cttcATAATCAGAACGGGGTTGCTACAACATGTATTAC‑TTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCACCGTAGTTGCCGCCAACGGTGAAGGTACATACCt    <  1:117255/150‑1 (MQ=255)
                                                                                                                                                    tcATAATCAGTACGGGGTTGCTACAACATGTATTAC‑TTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCACCGTAGTTGCCGCCAACGGTGAAGGTACATACCTGa  >  1:869629/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGTATCAGACAATAAGCATTCATTCAAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATAAATATTGGAATATATAAATATTGAATATTTTGATTAATCCCCGTACTGATTATTCTTCATAATCAGTACGGGGTTGCTACAACATGTATTAC‑TTTTTCACCACAATTAACGGTTCAATATCACTCTCTTTTTTGATGACCAGTGATTCATCACCGCGCAAACACGTCCCACCGTAGTTGCCGCCAACGGTGAAGGTACATACCTGA  >  NC_000913/3136483‑3136780

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: ‑