Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,618,945 A→T 5.6% F406I (TTC→ATC)  uraA ← uracil permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,618,9450AT 5.6% 618.8 / 19.5 179F406I (TTC→ATC) uraAuracil permease
Reads supporting (aligned to +/- strand):  ref base A (77/92);  new base T (2/8);  total (79/100)
Fisher's exact test for biased strand distribution p-value = 1.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.96e-01

Only 100 of 184 total aligned reads displayed.

TCACAAAATCGCGAGAAGTTTACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAGGTTCTGTGCTTTGTTGTAATCCACTTTC  >  NC_000913/2618796‑2619093
                                                                                                                                                     |                                                                                                                                                     
tCACAAAATCGCGAGAAGTTTACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGaa                                                                                                                                                       <  1:4118823/150‑1 (MQ=255)
  acaAAATCGCGAGAAGTTTACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGGGCAGCACACTGATCAGTTTGAAgt                                                                                                                                                     >  1:3854009/1‑149 (MQ=255)
  acaAAATCGCGAGAAGTTTACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGa                                                                                                                                                     >  2:2912400/1‑150 (MQ=255)
  acaAAATCGCGAGAAGTTTACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGa                                                                                                                                                     >  1:2037588/1‑150 (MQ=255)
  acaAAATCGCGAGAAGTTTACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGa                                                                                                                                                     >  1:2037617/1‑150 (MQ=255)
   caAAATCGCGAGAAGTTTACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGAt                                                                                                                                                    <  2:3740463/150‑1 (MQ=255)
    aaaaTCGCGAGAAGTTTACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATc                                                                                                                                                   <  1:924401/150‑1 (MQ=255)
     aaaTCGCGAGAAGTTTACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCa                                                                                                                                                  <  2:112569/150‑1 (MQ=255)
     aaaTCGCGAGAAGTTTACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCa                                                                                                                                                  <  2:1762176/150‑1 (MQ=255)
        tCGCGAGAAGTTTACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAAc                                                                                                                                               <  1:1347661/150‑1 (MQ=255)
            gagaAGTTTACCTCAGTATCCCATTATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCTCGTCTTCTGCGTCCATCTACACTTCTTCAGGGCGCAGCACACTGATCAGTTTGATCTTCAAACtgat                                                                                                                                           <  1:489410/150‑4 (MQ=255)
            gagaAGTTTACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATCTCCGCGTCTTCTGCGTCCAGCACCTCTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTaa                                                                                                                                           <  2:4129849/150‑1 (MQ=255)
            gagaAGTTTACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGATCTTCAAACtgat                                                                                                                                           <  1:83921/150‑4 (MQ=255)
            gagaAGTTTACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTaa                                                                                                                                           <  1:697897/150‑1 (MQ=255)
            gagaAGTTTACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTaa                                                                                                                                           <  1:2781960/150‑1 (MQ=255)
            gagaAGTTTACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTAGATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGATCTTCAAACtgat                                                                                                                                           <  1:1578509/150‑4 (MQ=255)
              gaAGTTTACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGc                                                                                                                                         <  1:3510538/150‑1 (MQ=255)
                aGTTTACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCg                                                                                                                                       >  1:2138720/1‑150 (MQ=255)
                aGTTTACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCg                                                                                                                                       >  1:3792658/1‑150 (MQ=255)
                  tttACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGAt                                                                                                                                     >  1:1006623/1‑150 (MQ=255)
                   ttACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATa                                                                                                                                    <  1:3496340/150‑1 (MQ=255)
                   ttACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATa                                                                                                                                    <  2:2592377/150‑1 (MQ=255)
                       cacaGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACcga                                                                                                                                <  2:2646748/150‑1 (MQ=255)
                         caGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACcgacg                                                                                                                              >  2:1222243/1‑150 (MQ=255)
                         caGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACcgacg                                                                                                                              >  2:1222189/1‑150 (MQ=255)
                          aGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACcgacga                                                                                                                             <  2:1483962/150‑1 (MQ=255)
                          aGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACcgacga                                                                                                                             <  2:1015500/150‑1 (MQ=255)
                          aGAAGCCCATGATGCTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTTTTTGACTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTACTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACcgacga                                                                                                                             <  2:1714351/150‑1 (MQ=255)
                                 cATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGc                                                                                                                      >  1:2197550/1‑150 (MQ=255)
                                    gatgTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTTTCTGTTATATAAGCGTCTTATGCGTACAGCACCACATCTTAAGGGCGCAGCACCAAGATAATTTTGAAGATAAAAATTAAGCCGATACCGAAAATGTTCGCCac                                                                                                                   >  1:354020/1‑149 (MQ=255)
                                      tgTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAgcg                                                                                                                 <  2:1314312/150‑1 (MQ=255)
                                      tgTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAgcg                                                                                                                 <  2:1852628/150‑1 (MQ=255)
                                       gTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAgcgc                                                                                                                >  1:2153312/1‑150 (MQ=255)
                                         ttGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCa                                                                                                              <  2:3321115/150‑1 (MQ=255)
                                         ttGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCa                                                                                                              <  1:3691841/150‑1 (MQ=255)
                                           gAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATa                                                                                                            <  2:3595799/150‑1 (MQ=255)
                                             aCCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATAcc                                                                                                          >  1:3018067/1‑150 (MQ=255)
                                                gggCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑ttt                                                                                                      <  1:2451514/150‑1 (MQ=255)
                                                 ggCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑tttt                                                                                                     <  2:703318/150‑1 (MQ=255)
                                                 ggCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑tttt                                                                                                     <  1:1025133/150‑1 (MQ=255)
                                                  gcagcaACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGAGGGTCGCCAGCGCCATACC‑TTTTa                                                                                                    <  2:216005/150‑1 (MQ=255)
                                                  gcagcaACACTGCCCGGTCGGGACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTa                                                                                                    >  1:3584189/1‑150 (MQ=255)
                                                  gcagcaACACTGCCAGGTCGGTACCATATTTGTGTGTTATATCCGCGTATTCAGCGTCCAGCACCACTTCTTCAGGGCGCCGCACACTGATCAGTTTGAAGATAAAACTTAAGCCGATACCCACGATGGTAGCCAGCGCCAAACC‑Ttata                                                                                                    >  1:3584225/1‑150 (MQ=255)
                                                                      gtACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGttt                                                                                >  1:3975829/1‑150 (MQ=255)
                                                                            tattTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTc                                                                          >  1:1912946/1‑150 (MQ=255)
                                                                            tattTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTc                                                                          <  2:3707372/150‑1 (MQ=255)
                                                                             attTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTcg                                                                         <  2:180612/150‑1 (MQ=255)
                                                                              cttGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTcgc                                                                        <  1:95381/149‑1 (MQ=255)
                                                                                tGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCgcgc                                                                      <  1:158594/150‑1 (MQ=255)
                                                                                tGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTGCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCgcgc                                                                      <  2:2853813/150‑1 (MQ=255)
                                                                                  tCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCa                                                                    <  1:3344394/150‑1 (MQ=255)
                                                                                   cTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCCGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCAc                                                                   <  1:1133470/150‑1 (MQ=255)
                                                                                    tGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACt                                                                  >  2:2505394/1‑150 (MQ=255)
                                                                                      ttATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGa                                                                <  2:3406219/150‑1 (MQ=255)
                                                                                      ttATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGa                                                                <  2:141700/150‑1 (MQ=255)
                                                                                      ttATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCAGATGTTAACCTACGCGCCACTGa                                                                >  2:1463700/1‑150 (MQ=255)
                                                                                       tataTCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGAc                                                               <  2:2731129/150‑1 (MQ=255)
                                                                                                 tcttctGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATc                                                     <  2:3850278/150‑1 (MQ=255)
                                                                                                   ttctGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCaa                                                   >  1:3900579/1‑150 (MQ=255)
                                                                                                     ctGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCaaaa                                                 >  1:2741053/1‑150 (MQ=255)
                                                                                                     ctGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCaaaa                                                 >  2:1940320/1‑150 (MQ=255)
                                                                                                      tGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAAt                                                <  2:3062822/150‑1 (MQ=255)
                                                                                                      tGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAAt                                                >  1:3768965/1‑150 (MQ=255)
                                                                                                        cGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCa                                              <  2:964276/150‑1 (MQ=255)
                                                                                                        cGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCa                                              <  2:1021049/150‑1 (MQ=255)
                                                                                                             aGCACCACTTCTTCAGGGCGCAGCACTCTGATCTGTTTGTAGATCAAAGTTATGCCGTTAGCGACGATGGTCGCCAGCGCCATAGC‑TTTTATCTCCGCCGCGCCGTTGTTTTCCTTCGCGCCTCTGACGGCGATGTTCTTATTCTCGGaa                                         <  2:86254/150‑1 (MQ=255)
                                                                                                             aGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCTCCAGCGCCATACCATTTAAACACCCCCGCCCCAATTTTAACGTACGCCCCACTGACCCCAAGAAGCAAAATCACGGa                                          >  1:2484300/1‑150 (MQ=255)
                                                                                                             aGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGaa                                         >  1:2688114/1‑150 (MQ=255)
                                                                                                             aGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGaa                                         >  1:2574449/1‑150 (MQ=255)
                                                                                                                accacTTCTTCCGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTc                                      <  1:2443773/150‑1 (MQ=255)
                                                                                                                 ccacTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCa                                     <  2:2257140/150‑1 (MQ=255)
                                                                                                                 ccacTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAAACACGGAAGTCa                                     >  2:3101924/1‑150 (MQ=255)
                                                                                                                  cacTTCTTCAGGGCGCAACACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAg                                    >  2:391781/1‑150 (MQ=255)
                                                                                                                   acTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAgg                                   >  1:1207140/1‑150 (MQ=255)
                                                                                                                    cttcttCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGa                                  >  1:445222/1‑150 (MQ=255)
                                                                                                                     ttcttcAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTGTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGag                                 >  2:2396711/1‑149 (MQ=255)
                                                                                                                        ttcAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCAGGGAAGTCAGGATCAg                              >  2:730719/1‑150 (MQ=255)
                                                                                                                        ttcAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAg                              >  1:2512902/1‑150 (MQ=255)
                                                                                                                        ttcAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAg                              >  1:3528342/1‑150 (MQ=255)
                                                                                                                        ttcAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAg                              >  2:40732/1‑150 (MQ=255)
                                                                                                                        ttcAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAg                              >  1:3188374/1‑150 (MQ=255)
                                                                                                                        ttcAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAg                              >  2:106336/1‑150 (MQ=255)
                                                                                                                         tcAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAgg                             >  2:3857133/1‑150 (MQ=255)
                                                                                                                         tcAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCAGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAgg                             >  2:208552/1‑150 (MQ=255)
                                                                                                                               cgcAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAGGTTCtgt                       <  1:2369493/150‑1 (MQ=255)
                                                                                                                               cgcAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAGGTTCtgt                       <  1:1694359/150‑1 (MQ=255)
                                                                                                                               cgcAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAGGTTCtgt                       <  1:1694634/150‑1 (MQ=255)
                                                                                                                                gcagcaCACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGAGGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAGGTTCtgtg                      <  2:3284587/150‑1 (MQ=255)
                                                                                                                                 cagcaCACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAGGTTCTGTGc                     <  1:2780412/150‑1 (MQ=255)
                                                                                                                                 cagcaCACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAGGTTCTGTGc                     <  2:2700614/150‑1 (MQ=255)
                                                                                                                                   gcaCACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAGGTTCTGTGCtt                   >  1:1251180/1‑150 (MQ=255)
                                                                                                                                    cacaCTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAGGTTCTGTGCttt                  <  1:1501792/150‑1 (MQ=255)
                                                                                                                                       acTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAGGTTCTGTGCTttgtt               >  2:1297960/1‑150 (MQ=255)
                                                                                                                                       acTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAGGTTCTGTGCTttgtt               >  2:217563/1‑150 (MQ=255)
                                                                                                                                          cttaaGTTTGATCTTCAAACTTAAGCCGATACCGACGATGGTCGCAAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAGGTTCTGTGCTTTGTTGTa            >  2:4043879/5‑150 (MQ=255)
                                                                                                                                            tCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAGGTTCTGTGCTTTGTTGTAAt          <  2:3404868/150‑1 (MQ=255)
                                                                                                                                                 ttGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAGGTTCTGTGCTTTGTTGTAATCCACt     >  1:2966780/1‑150 (MQ=255)
                                                                                                                                                  tgtaGATCAAAGTTAAGCCGATACCGACGTTGGTCGCCTGCGCCATACG‑TTTTAACTGGGCCGCGCCGATGTTTACCTTCGCGCGACTGACGCCGATGATCAAATTCTGGGTTGTCAGGATCTGGTTCTGTGCTTTGTTGTATTCCAGtt    >  2:3097507/4‑150 (MQ=255)
                                                                                                                                                  tGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAGGTTCTGTGCTTTGTTGTAATCCACtt    >  2:1665114/1‑150 (MQ=255)
                                                                                                                                                    aaGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAGGTTCTGTGCTTTGTTGTAATCCACTttc  <  2:3975829/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCACAAAATCGCGAGAAGTTTACCACAGAATCCCATGATGTTTGAACCGGGCAGCAACACTGCCCGGTCGGTACATTATTTGTCTGTTATATCCGCGTCTTCTGCGTCCAGCACCACTTCTTCAGGGCGCAGCACACTGATCAGTTTGAAGATCAAACTTAAGCCGATACCGACGATGGTCGCCAGCGCCATACC‑TTTTAACTCCGCCGCGCCGATGTTTACCTTCGCGCCACTGACGCCGATGATCAAAATCACGGAAGTCAGGATCAGGTTCTGTGCTTTGTTGTAATCCACTTTC  >  NC_000913/2618796‑2619093

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑