Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,258,252 T→A 6.2% intergenic (‑242/+33) cyuP ← / ← tdcG putative D/L‑serine transporter/L‑serine deaminase III

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,258,2520TA 6.2% 311.3 / 7.8 97intergenic (‑242/+33)cyuP/tdcGputative D/L‑serine transporter/L‑serine deaminase III
Reads supporting (aligned to +/- strand):  ref base T (49/42);  new base A (4/2);  total (53/44)
Fisher's exact test for biased strand distribution p-value = 6.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.36e-01

GAATAAAGAAATGATACGCGCGAGGTCAGAGAAAATCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATACATCGTCTCGATCACTTTATCGAGTGAAACACG  >  NC_000913/3258103‑3258398
                                                                                                                                                     |                                                                                                                                                  
gAATAAAGAAATGATACGCGCGAGGTCAGAGAAAATCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAACGGTGCACAtt                                                                                                                                                    >  1:918076/1‑150 (MQ=255)
  aTAAAGAAATGATACGCGCGAGGTCAGAGAAAATCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAACGGTGCACATTtg                                                                                                                                                  >  2:1767143/1‑150 (MQ=255)
  aTAAAGAAATGATACGCGCGAGGTCAGAGAAAATCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTtg                                                                                                                                                  <  2:552858/150‑1 (MQ=255)
  aTAAAGAAATGATACGCGCGAGGTCAGAGAAAATCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTtg                                                                                                                                                  <  2:1420525/150‑1 (MQ=255)
     aaGAAATGATACGCGCGAGGTCAGAGAAAATCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGc                                                                                                                                               <  2:289903/150‑1 (MQ=255)
     aaGAAATGATACGAGCGAGGTCAGAGAAAATCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGAAGGAGAATATTCATTCTATATTTGCTGAAGATAACAAAAAAAAGGTGCACATTTGTGc                                                                                                                                               >  1:506194/1‑150 (MQ=255)
      aGAAATGATACGCGCGAGGTCAGAGAAAATCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCa                                                                                                                                              >  2:1756326/1‑150 (MQ=255)
      aGAAATGATACGCGCGAGGTCAGAGAAAATCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACAAATGTGCa                                                                                                                                              <  1:1962593/150‑1 (MQ=255)
       taaaTGATACGCGCGAGGTCAGAGAAAATCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCAc                                                                                                                                             <  2:1839699/149‑1 (MQ=255)
       gAAATGATACGCGCGAGGTCAGAGAAAATCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCAc                                                                                                                                             >  2:1514301/1‑150 (MQ=255)
          aTGATACGCGCGAGGTCAGAGAAAATCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACAAATGTGCAcctt                                                                                                                                          <  2:1125535/150‑3 (MQ=255)
            gATACGCGCGAGGTCAGAGAAAATCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAg                                                                                                                                        <  1:994731/150‑1 (MQ=255)
            gATACGCGCGAGGTCAGAGAAAATCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAg                                                                                                                                        >  2:174639/1‑150 (MQ=255)
             aTACGCGCGATGTCAGAGAAAATCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAATAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAgg                                                                                                                                       <  2:852629/150‑1 (MQ=255)
                cgcgcgAGGTCAGAGAAAATCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATg                                                                                                                                    <  2:1063292/150‑1 (MQ=255)
                cgcgcgAGGTCAGAGAAAATCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATg                                                                                                                                    <  1:1281650/150‑1 (MQ=255)
                   gcgAGGTCAGAGAAAATCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGaaa                                                                                                                                 <  1:60908/150‑1 (MQ=255)
                         tCAGAGAAAATCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGAc                                                                                                                           <  1:72778/150‑1 (MQ=255)
                          cAGAGAAAATCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACa                                                                                                                          <  1:1489581/150‑1 (MQ=255)
                          cAGAGAAAATCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACa                                                                                                                          <  2:160774/150‑1 (MQ=255)
                            gagaAAATCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGc                                                                                                                        <  2:1026028/150‑1 (MQ=255)
                                aaaTCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATg                                                                                                                    >  2:1260565/1‑150 (MQ=255)
                                aaaTCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATg                                                                                                                    <  2:1631819/150‑1 (MQ=255)
                                  aTCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTc                                                                                                                  >  2:71917/1‑150 (MQ=255)
                                      tCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGcc                                                                                                              >  1:1642431/1‑150 (MQ=255)
                                      tCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGcc                                                                                                              >  2:69142/1‑150 (MQ=255)
                                         caATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAATCTGACAGCAATGTCAGCCGCa                                                                                                           >  1:1062373/1‑150 (MQ=255)
                                          aaTCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCGCATTTGTGCACCCAAGGATGAACGCTGACAGCAATGTCAGCCGCAg                                                                                                          >  2:322913/1‑150 (MQ=255)
                                           aTCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAgc                                                                                                         >  2:672712/1‑149 (MQ=255)
                                           aTCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGa                                                                                                         >  1:507789/1‑150 (MQ=255)
                                            tCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCATGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGAc                                                                                                        <  2:376283/150‑1 (MQ=255)
                                            tCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGAc                                                                                                        >  1:390795/1‑150 (MQ=255)
                                              aactcTTTGCAATAACGGACCTTGCAAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCa                                                                                                      <  1:1012100/146‑1 (MQ=255)
                                              aaGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCa                                                                                                      >  2:661050/1‑150 (MQ=255)
                                                  cTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACttt                                                                                                  <  2:1594648/150‑1 (MQ=255)
                                                     tGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAAt                                                                                               >  2:1443882/1‑150 (MQ=255)
                                                     tGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACAACTTTAAt                                                                                               >  2:129168/1‑150 (MQ=255)
                                                      gAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATg                                                                                              >  2:1704503/1‑150 (MQ=255)
                                                          tAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCa                                                                                          <  1:1756326/150‑1 (MQ=255)
                                                           aaCGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAg                                                                                         >  2:584454/1‑150 (MQ=255)
                                                               gagaTGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGtcct                                                                                     >  1:601904/1‑150 (MQ=255)
                                                               gagaTGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGtcct                                                                                     >  1:1574501/1‑150 (MQ=255)
                                                                   tGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGAcc                                                                                                       <  1:1923158/128‑1 (MQ=255)
                                                                   tGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGAcc                                                                                                       >  2:1923158/1‑128 (MQ=255)
                                                                   tGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCcgc                                                                                 >  1:226155/1‑150 (MQ=255)
                                                                    ggATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGGTCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCcgcg                                                                                >  2:1682313/1‑150 (MQ=255)
                                                                    ggATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGGTCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCcgcg                                                                                >  2:1689262/1‑150 (MQ=255)
                                                                      aTAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTg                                                                              >  2:64762/1‑150 (MQ=255)
                                                                       tAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTTACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGa                                                                             <  2:765103/150‑1 (MQ=255)
                                                                         aGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATg                                                                           >  2:1670760/1‑150 (MQ=255)
                                                                           aaTTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGtt                                                                         <  1:174639/150‑1 (MQ=255)
                                                                           aaTTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGtt                                                                         >  1:898429/1‑150 (MQ=255)
                                                                           aaTTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGtt                                                                         <  2:507789/150‑1 (MQ=255)
                                                                              tttCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTcg                                                                      <  1:1597181/150‑1 (MQ=255)
                                                                                tCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTcgcg                                                                    <  2:1642431/150‑1 (MQ=255)
                                                                                  tactaAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGt                                                                  >  2:728363/1‑150 (MQ=255)
                                                                                   actaAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTa                                                                 >  1:1254975/1‑150 (MQ=255)
                                                                                    ctaAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTAt                                                                >  1:108843/1‑150 (MQ=255)
                                                                                    ctaAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTAt                                                                >  1:612440/1‑150 (MQ=255)
                                                                                     taAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTAtt                                                               >  1:700590/1‑150 (MQ=255)
                                                                                       aaTTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTa                                                             <  1:1704503/150‑1 (MQ=255)
                                                                                        aTTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAATAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACGACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTAt                                                            <  1:129168/150‑1 (MQ=255)
                                                                                         ttAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATc                                                           <  1:1181561/150‑1 (MQ=255)
                                                                                            aTCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGtt                                                        <  1:584454/150‑1 (MQ=255)
                                                                                              cGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCa                                                      >  2:1281405/1‑150 (MQ=255)
                                                                                               gCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCat                                                     >  2:18443/1‑150 (MQ=255)
                                                                                                 aGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCatat                                                   <  1:1232912/150‑1 (MQ=255)
                                                                                                    agaATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTTTTGTTTCGCGGTATTTATCGTTCATATCtt                                                >  2:891523/1‑150 (MQ=255)
                                                                                                    agaATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCtt                                                >  2:1452500/1‑150 (MQ=255)
                                                                                                       atatTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGc                                             >  2:219554/1‑150 (MQ=255)
                                                                                                       atatTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCATTGc                                             <  1:586237/150‑1 (MQ=255)
                                                                                                        tatTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCGGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGcc                                            >  2:1177109/1‑150 (MQ=255)
                                                                                                        tatTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGcc                                            >  2:28164/1‑150 (MQ=255)
                                                                                                         attcattcTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCg                                           <  1:147584/150‑1 (MQ=255)
                                                                                                          ttcattcTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTGGCGGTATTTATCGTTCATATCTTTGCCgg                                          >  1:59649/1‑150 (MQ=255)
                                                                                                            cattcTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAATCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGtt                                        <  2:1062373/150‑1 (MQ=255)
                                                                                                                ctctATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCAt                                    <  1:1260565/150‑1 (MQ=255)
                                                                                                                    aTTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATACAt                                <  1:196004/150‑1 (MQ=255)
                                                                                                                      ttGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATACATCg                              <  1:18443/150‑1 (MQ=255)
                                                                                                                         cGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATACATCGtct                           >  2:1617364/1‑150 (MQ=255)
                                                                                                                          ggCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATACATCGtctc                          >  2:1046659/1‑150 (MQ=255)
                                                                                                                             aGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATACATCGTCTCGAt                       <  2:1574501/150‑1 (MQ=255)
                                                                                                                              gATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATACATCGTCTCGATc                      >  2:529220/1‑150 (MQ=255)
                                                                                                                               aTCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATACATCGTCTCGATCa                     <  2:601904/150‑1 (MQ=255)
                                                                                                                                 cacaAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATACATCGTCTCGATCACt                   >  1:255419/1‑150 (MQ=255)
                                                                                                                                      tccttgggTGCACAAATGTGCACCTAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATACATCGTCTCGATCACTTTATc              >  1:671370/7‑150 (MQ=255)
                                                                                                                                      tccttgggTGCACAAATGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATACATCGTCTCGATCACTTTATc              >  2:1456304/7‑150 (MQ=255)
                                                                                                                                       aaaaaGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATACATCGTCTCGATCACTTTATCg             <  1:1177109/150‑1 (MQ=255)
                                                                                                                                        cttgggTGCACAAATGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTCTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATACATCGTCTCGATCACTTTATCGa            >  1:1343615/5‑150 (MQ=255)
                                                                                                                                          aaGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGTCGGTTTCATACATCGTCTCGATCACTTTATCGAGt          <  2:898429/150‑1 (MQ=255)
                                                                                                                                            ggTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATACATCGTCTCGATCACTTTATCGAGTGa        <  1:1452500/150‑1 (MQ=255)
                                                                                                                                            ggTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCAGCAGACCACTTTAATGGCCAGTCCTCCGCTTGTTGTTTCGCTGTATTTATCGTTCATATCTTTGCCGGTTTCATACATCGTCTCGATCACTTTATCGAGTGa        <  1:891523/150‑1 (MQ=255)
                                                                                                                                            ggTGCACAAATGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATACATCGTCTCGATCACTTTATCGAGTGa        >  2:1336582/1‑150 (MQ=255)
                                                                                                                                             gTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATACATCGTCTCGATCACTTTATCGAGTGaa       <  1:728363/150‑1 (MQ=255)
                                                                                                                                               gCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATACATCGTCTCGATCACTTTATCGAGTGAAac     >  2:368027/1‑150 (MQ=255)
                                                                                                                                                  caTTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATACATCGTCTCGATCACTTTATCGAGTGAAacacc  >  2:522668/1‑149 (MQ=255)
                                                                                                                                                  caTTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATACATCGTCTCGATCACTTTATCGAGTGAAACAcg  >  2:253280/1‑150 (MQ=255)
                                                                                                                                                  caTTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATACATCGTCTCGATCACTTTATCGAGTGAAACAcg  >  2:171147/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GAATAAAGAAATGATACGCGCGAGGTCAGAGAAAATCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATACATCGTCTCGATCACTTTATCGAGTGAAACACG  >  NC_000913/3258103‑3258398

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