Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,289,573 C→T 5.3% A110V (GCG→GTG)  rssA → putative patatin‑like family phospholipase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,289,5730CT 5.3% 393.7 / 7.5 112A110V (GCG→GTG) rssAputative patatin‑like family phospholipase
Reads supporting (aligned to +/- strand):  ref base C (56/50);  new base T (3/3);  total (59/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 113 total aligned reads displayed.

TTGGGTGACCTCTTTCAGTTATTGGGATGTTTTACGCCTGATGGATCTCTCCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCCGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGTATTCCAGGACTCATGGCACCTGTTGCACATAACGGCTACTGGCTGGTTGATGGAGCAGTC  >  NC_000913/1289424‑1289718
                                                                                                                                                     |                                                                                                                                                 
ttGGGTGACCTCTTTCAGTTATTGGGATGTTTTACGCCTGATGGATCTCTCCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCCGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGc                                                                                                                                                   >  1:1630324/1‑150 (MQ=255)
      gACCTCTTTCAGTTATTGGGATGTTTTACGCCTGATGGATCTCTCCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAAAACCGCGAAATAATGCCGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGt                                                                                                                                             >  1:1716952/1‑150 (MQ=255)
           ctTTCAGTTATTGGGATGTTTTACGCCTGATGGATCTCTCCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCCGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGcca                                                                                                                                        <  2:1003402/150‑1 (MQ=255)
                aGTTATTGGGATGTTTTACGCCTGATGGATCTCTCCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCCGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAAt                                                                                                                                   >  2:1101683/1‑150 (MQ=255)
                aGTTATTGGGATGTTTTACGCCTGATGGATCTCTCCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCCGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAAt                                                                                                                                   >  2:2184395/1‑150 (MQ=255)
                        ggATGTTTTACGCCTGATGGATCTCTCCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCCGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTAc                                                                                                                           >  2:2216689/1‑150 (MQ=255)
                        ggATGTTTTACGCCTGATGGATCTCTCCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCCGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTAc                                                                                                                           >  2:380774/1‑150 (MQ=255)
                        ggATGTTTTACGCCTGATGGATCTCTCCTGGCAGCGCGGAGGGTTACTGCGCGGCGAGCGAGTCTTCAATCAATACAGCGAAATAATGCCGGAAACAGTGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTAAAGTAc                                                                                                                           >  1:587049/1‑150 (MQ=255)
                         gATGTTTTACGCCTGATGGATCTCTCCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCCGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACg                                                                                                                          >  1:1947107/1‑150 (MQ=255)
                         gATGTTTTACGCCTGATGGATCTCTCCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCCGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACg                                                                                                                          <  1:2153272/150‑1 (MQ=255)
                         gATGTTTTACGCCTGATGGATCTCTCCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCCGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACg                                                                                                                          <  1:2152804/150‑1 (MQ=255)
                             ttttACGCCTGATGGATCTCTCCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCCGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGAc                                                                                                                      >  2:1726923/1‑150 (MQ=255)
                             ttttACGCCTGATGGATCTCTCCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCCGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGAc                                                                                                                      >  2:1164107/1‑150 (MQ=255)
                                   gCCTGATGGATCTCTCCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCCGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAAt                                                                                                                <  2:1101026/150‑1 (MQ=255)
                                     cTGATGGATCTCTCCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCCGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTa                                                                                                              >  1:2110083/1‑150 (MQ=255)
                                      tGATGGATCTCTCCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCCGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTAt                                                                                                             >  1:1559615/1‑150 (MQ=255)
                                      tGATGGATCTCTCCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCCGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTAt                                                                                                             <  1:1209799/150‑1 (MQ=255)
                                      tGATGGATCTCTCCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCCGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTAt                                                                                                             <  1:336893/150‑1 (MQ=255)
                                             tctctcCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCCGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTAc                                                                                                      <  2:448129/150‑1 (MQ=255)
                                             tctctcCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCCGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTAc                                                                                                      <  2:616963/150‑1 (MQ=255)
                                               tctcCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCCGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTg                                                                                                    >  1:1807754/1‑150 (MQ=255)
                                               tctcCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCCGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTg                                                                                                    >  1:950968/1‑150 (MQ=255)
                                                 tcCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCCGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGaa                                                                                                  >  2:2200782/1‑150 (MQ=255)
                                                 tcCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCCGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGaa                                                                                                  <  2:146264/150‑1 (MQ=255)
                                                 tcCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCCGGAAACAGAGATCGAAAATTGTTCCCGTCGCTATGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGaa                                                                                                  >  2:2200805/1‑150 (MQ=255)
                                                 tcCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCCGGAAACAGAGATCGAAAATTGTTCCCGACGCTTTGCGGCTGTTGCCACCACTTTAAGTACGGGACGTGAAGTATGGTTTACTGaa                                                                                                  >  2:2200783/1‑150 (MQ=255)
                                                   cTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCCGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAgg                                                                                                >  1:1268991/1‑150 (MQ=255)
                                                   cTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCCGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAgg                                                                                                >  1:1268930/1‑150 (MQ=255)
                                                      gCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCCGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGa                                                                                             >  1:675073/1‑150 (MQ=255)
                                                      gCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCCGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGa                                                                                             >  1:675093/1‑150 (MQ=255)
                                                      gCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCCGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGa                                                                                             >  1:2072227/1‑150 (MQ=255)
                                                      gCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCCGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGa                                                                                             <  1:1468456/150‑1 (MQ=255)
                                                        aGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCCGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGAtc                                                                                           <  1:2184395/150‑1 (MQ=255)
                                                        aGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCCGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGAtc                                                                                           <  2:1574369/150‑1 (MQ=255)
                                                        aGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCCGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGAtc                                                                                           <  2:1329063/150‑1 (MQ=255)
                                                                   ttACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCCGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCt                                                                                >  1:444192/1‑150 (MQ=255)
                                                                       tGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCCGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTc                                                                            <  1:954662/150‑1 (MQ=255)
                                                                         cgcgGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCCGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTcgc                                                                          >  1:1788597/1‑150 (MQ=255)
                                                                              cgAGCGTGTCTTCAATCAATACCGCGAAATAATGCCGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTCGCGcatc                                                                     <  2:1499777/150‑1 (MQ=255)
                                                                               gAGCGTGTCTTCAATCAATACCGCGAAATAATGCCGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTCGCGcatca                                                                    >  2:1026485/1‑150 (MQ=255)
                                                                                  cGTGTCTTCAATCAATACCGCGAAATAATGCCGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGc                                                                 <  2:1006716/150‑1 (MQ=255)
                                                                                    tgtCTTCAATCAATACCGCGAAATAATGCCGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAg                                                               <  2:1947107/150‑1 (MQ=255)
                                                                                        ttCAATCAATACCGCGAAATAATGCCGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGTAtt                                                           >  2:464542/1‑150 (MQ=255)
                                                                                         tcaatcaaTACCGCGAAATAATGCCGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGTATTc                                                          <  2:643105/150‑1 (MQ=255)
                                                                                         tcaatcaaTACCGCGAAATAATGCCGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGTATTc                                                          <  2:686661/150‑1 (MQ=255)
                                                                                          caatcaatACCGCGAAATAATGCCGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGTATTcc                                                         >  2:1138829/1‑150 (MQ=255)
                                                                                          caatcaatACCGCGAAATAATGCCGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGTATTcc                                                         >  2:1910213/1‑150 (MQ=255)
                                                                                          caatcaatACCGCGAAATAATGCCGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGTATTcc                                                         >  2:1973554/1‑150 (MQ=255)
                                                                                           aatcaatACCGCGAAATAATGCCGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGTATTCCa                                                        <  2:1954030/150‑1 (MQ=255)
                                                                                           aatcaatACCGCGAAATAATGCCGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGTATTCCa                                                        <  2:969533/150‑1 (MQ=255)
                                                                                           aatcaatACCGCGAAATAATGCCGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGTATTCCa                                                        <  1:1839712/150‑1 (MQ=255)
                                                                                           aatcaatACCGCGAAATAATGCCGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGTATTCCa                                                        <  1:379276/150‑1 (MQ=255)
                                                                                           aatcaatACCGCGAAATAATGCCGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGTATTCCa                                                        <  1:1839461/150‑1 (MQ=255)
                                                                                            atcaatACCGCGAAATAATGCCGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGTATTCCAg                                                       >  2:2021791/1‑150 (MQ=255)
                                                                                              caatACCGCGAAATAATGCCGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGTATTCCAGGa                                                     <  2:1616944/150‑1 (MQ=255)
                                                                                              caatACCGCGAAATAATGCCGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGTATTCCAGGa                                                     <  2:1616889/150‑1 (MQ=255)
                                                                                                  aCCGCGAAATAATGCCGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGTATTCCAGGACTCa                                                 >  2:53061/1‑150 (MQ=255)
                                                                                                  aCCGCGAAATAATGCCGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGTATTCCAGGACTCa                                                 >  2:52946/1‑150 (MQ=255)
                                                                                                  aCCGCGAAATAATGCCGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGTATTCCAGGACTCa                                                 >  2:52931/1‑150 (MQ=255)
                                                                                                    cgcgAAATAATGCCGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATTGTTTACTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGTATTCCAGGACTCATg                                               <  1:227993/150‑1 (MQ=255)
                                                                                                    cgcgAAATAATGCCGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGTATTCCAGGACTCATg                                               <  1:146514/150‑1 (MQ=255)
                                                                                                      cgAAATAATGCCGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGTATTCCAGGACTCATGGc                                             <  2:2231891/150‑1 (MQ=255)
                                                                                                      cgAAATAATGCCGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGTATTCCAGGACTCATGGc                                             <  2:1477557/150‑1 (MQ=255)
                                                                                                        aaaTAATGCCGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGTATTCCAGGACTCATGGCAc                                           <  1:270065/150‑1 (MQ=255)
                                                                                                        aaaTAATGCCGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGTATTCCAGGACTCATGGCAc                                           <  1:270029/150‑1 (MQ=255)
                                                                                                        aaaTAATGCCGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGTATTCCAGGACTCATGGCAc                                           <  2:1179786/150‑1 (MQ=255)
                                                                                                        aaaTAATGCCGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGTATTCCAGGACTCATGGCAc                                           <  2:497692/150‑1 (MQ=255)
                                                                                                        aaaTAATGCCGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGTATTCCAGGACTCATGGCAc                                           <  1:1001873/150‑1 (MQ=255)
                                                                                                         aataatGCCGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTAGTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGTATTCCAGGACTCTTGGCAcc                                          >  2:803880/1‑150 (MQ=255)
                                                                                                         aataatGCCGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGTATTCCAGGACTCATGGCAcc                                          >  2:1975186/1‑150 (MQ=255)
                                                                                                              tGCCGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGTATTCCAGGACTCATGGCACCTGTTg                                     >  2:1130566/1‑150 (MQ=255)
                                                                                                                ccGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGTATTCCAGGACTCATGGCACCTGTTGca                                   <  1:2039245/150‑1 (MQ=255)
                                                                                                                 cGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGTATTCCAGGACTCATGGCACCTGTTGcac                                  >  1:20296/1‑150 (MQ=255)
                                                                                                                     aaCAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGTATTCCAGGACTCATGGCACCTGTTGCACATaa                              <  1:884599/150‑1 (MQ=255)
                                                                                                                     aaCAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGTATTCCAGGACTCATGGCACCTGTTGCACATaa                              <  1:1693166/150‑1 (MQ=255)
                                                                                                                       cAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGTATTCCAGGACTCATGGCACCTGTTGCACATAACg                            >  1:69685/1‑150 (MQ=255)
                                                                                                                       cAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGTATTCCAGGACTCATGGCACCTGTTGCACATAACg                            >  1:1590434/1‑150 (MQ=255)
                                                                                                                        agagATCGAAAATTGTTCCCGTCGCTTTGTGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGTATTCCAGGACTCATGGCACCTGTTGCACATAACgg                           >  2:663485/1‑150 (MQ=255)
                                                                                                                        agagATCGAAAATTGTTCCCGTCGCTTTGTGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGTATTCCAGGACTCATGGCACCTGTTGCACATAACgg                           >  2:29103/1‑150 (MQ=255)
                                                                                                                        agagATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGTATTCCAGGACTCATGGCACCTGTTGCACATAACgg                           >  1:873040/1‑150 (MQ=255)
                                                                                                                        agagATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGTATTCCAGGACTCATGGCACCTGTTGCACATAACgg                           <  1:1726923/150‑1 (MQ=255)
                                                                                                                           gaTCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGTATTCCAGGACTCATGGCACCTGTTGCACATAACGGCTa                        >  2:1866542/1‑150 (MQ=255)
                                                                                                                           gaTCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGTATTCCAGGACTCATGGCACCTGTTGCACATAACGGCTa                        <  1:2256772/150‑1 (MQ=255)
                                                                                                                           gaTCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGTATTCCAGGACTCATGGCACCTGTTGCACATAACGGCTa                        >  2:670572/1‑150 (MQ=255)
                                                                                                                             tCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGTATTCCAGGACTCATGGCACCTGTTGCACATAACGGCTAct                      >  2:1034571/1‑150 (MQ=255)
                                                                                                                               gAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATTGTTTACTGAAGGCGATCTCCATCTTGCTATTAGCGCATCATGCAGTATTACAGGACTCAAGGCACCAGTTGCACATAACGGCTActgg                    >  1:1024835/1‑150 (MQ=255)
                                                                                                                               gAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGTATTCCAGGACTCATGGCACCTGTTGCACATAACGGCTActgg                    >  2:1166381/1‑150 (MQ=255)
                                                                                                                               gAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGTATTCCAGGACTCATGGCACCTGTTGCACATAACGGCTActgg                    >  2:2219645/1‑150 (MQ=255)
                                                                                                                               gAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGTATTCCAGGACTCATGGCACCTGTTGCACATAACGGCTActgg                    >  2:1805003/1‑150 (MQ=255)
                                                                                                                               gAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTAGCGCATAATGCAGTATTCCAGGACTCATGGCACCAGTTTCACATAACGGCTActgg                    >  1:1138194/1‑150 (MQ=255)
                                                                                                                               gAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAACTACGGGACGTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGTATTCCAGGACTCAAGGCACCTGTTGCACATAACGGCTActgg                    >  1:1290768/1‑150 (MQ=255)
                                                                                                                                 aaaTTGGTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGTATTCCAGGACTCATGGCACCTGTTGCACATAACGGCTActggct                  <  1:2041090/150‑1 (MQ=255)
                                                                                                                                   aTTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGTATTCCAGGACTCATGGCACCTGTTGCACATAACGGCTActggctgg                <  2:1077125/150‑1 (MQ=255)
                                                                                                                                     tggtCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGTATTCCAGGACTCATGGCACCTGTTGCACATAACGGCTACTGGCTGGtt              <  2:1716952/147‑1 (MQ=255)
                                                                                                                                         cccGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGTATTCCAGGACTCATGGCACCTGTTGCACATAACGGCTACTGGCTGGTTGATg          <  1:2303841/150‑1 (MQ=255)
                                                                                                                                           cGTCGCTTTGTGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGTATTCCAGGACTCATGGCACCTGTTGCACATAACGGCTACTGGCTGGTTGATGGa        <  1:29133/150‑1 (MQ=255)
                                                                                                                                           cGTCGCTTTGTGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGTATTCCAGGACTCATGGCACCTGTTGCACATAACGGCTACTGGCTGGTTGATGGa        <  1:663485/150‑1 (MQ=255)
                                                                                                                                           cGTCGCTTTGTGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGTATTCCAGGACTCATGGCACCTGTTGCACATAACGGCTACTGGCTGGTTGATGGa        <  1:29103/150‑1 (MQ=255)
                                                                                                                                           cGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAATGCGATCTCCATCTTTCTATTCGCGCATCATGCAGTATTCCAGGACTCATGGCACCTGTTGCACATAACGGCTACTGGCTGGTTGATGGa        >  1:1122848/1‑150 (MQ=255)
                                                                                                                                                 tttGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGTATTCCAGGACTCATGGCACCTGTTGCACATAACGGCTACTGGCTGGTTGATGGAGCAGTc  >  2:796575/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
TTGGGTGACCTCTTTCAGTTATTGGGATGTTTTACGCCTGATGGATCTCTCCTGGCAGCGCGGTGGGTTACTGCGCGGCGAGCGTGTCTTCAATCAATACCGCGAAATAATGCCGGAAACAGAGATCGAAAATTGTTCCCGTCGCTTTGCGGCTGTTGCCACCAATTTAAGTACGGGACGTGAATTATGGTTTACTGAAGGCGATCTCCATCTTGCTATTCGCGCATCATGCAGTATTCCAGGACTCATGGCACCTGTTGCACATAACGGCTACTGGCTGGTTGATGGAGCAGTC  >  NC_000913/1289424‑1289718

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