Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,355,187 A→G 5.4% V95A (GTG→GCG) ‡ sapB ← putrescine ABC exporter membrane subunit SapB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,355,1870AG 5.4% 337.5 / 5.3 109V95A (GTG→GCG) ‡sapBputrescine ABC exporter membrane subunit SapB
Reads supporting (aligned to +/- strand):  ref base A (45/58);  new base G (4/2);  total (49/60)
Fisher's exact test for biased strand distribution p-value = 4.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 110 total aligned reads displayed.

GAAACGGGCAACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAATGCTGGAAACACCAAAATCCCAGTGGATCAGGCCGTTAAACCAGAAC  >  NC_000913/1355040‑1355336
                                                                                                                                                   |                                                                                                                                                     
gAAACGGGCAACCAACCTAGCCTGAGTGAACAAAACAGCGTCAACAGCAGCGCCAGCCAGAACACCGCGATTGAAAAGCCAAGTAACGCGATGGCGTTGATCAAGTAATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACg                                                                                                                                                     >  2:1162068/1‑150 (MQ=255)
  aaCGGGCAACCAACCTAGGGTCTTTTAACAAAACAGCGTCAACAGAAGCCCCAGCCAGAACACCGGCATTCAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTTTCGTGCCAGTTTTGGCGCGTAATCCGGCAATCATCCCCGTgggg                                                                                                                                                  <  2:158101/150‑1 (MQ=255)
  aaCGGGCAACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACggg                                                                                                                                                   <  2:1386723/150‑1 (MQ=255)
  aaCGGGCAACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACggg                                                                                                                                                   >  1:1172074/1‑150 (MQ=255)
   aCGGGCAACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACgggg                                                                                                                                                  <  2:1176351/150‑1 (MQ=255)
   aCGGGCAACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACgggg                                                                                                                                                  <  1:404962/150‑1 (MQ=255)
   aCGGGCAACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACgggg                                                                                                                                                  <  2:564054/150‑1 (MQ=255)
     gggCAACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGAt                                                                                                                                                >  2:719137/1‑150 (MQ=255)
     gggCAACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGAt                                                                                                                                                <  1:169137/150‑1 (MQ=255)
        caaccaacCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGcc                                                                                                                                             <  2:1284828/150‑1 (MQ=255)
        caaccaacCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGcc                                                                                                                                             <  2:1343319/150‑1 (MQ=255)
        caaccaacCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGcc                                                                                                                                             <  1:785396/150‑1 (MQ=255)
            caaccTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCcaaca                                                                                                                                         >  2:1248328/1‑150 (MQ=255)
             aaccTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCcaacaa                                                                                                                                        <  1:870509/150‑1 (MQ=255)
               ccTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATc                                                                                                                                      <  2:502895/150‑1 (MQ=255)
                cTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCa                                                                                                                                     >  1:1609917/1‑150 (MQ=255)
                 tAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCACTCATCCCCACGGGGATGCCAACAATCAg                                                                                                                                    >  1:828905/1‑150 (MQ=255)
                   gCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGTGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAgcg                                                                                                                                  <  1:1383607/150‑1 (MQ=255)
                     gtgagtgaACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCa                                                                                                                                >  1:504129/1‑150 (MQ=255)
                       gagtgaACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCaaa                                                                                                                              >  2:319848/1‑150 (MQ=255)
                            aaCAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGa                                                                                                                         >  1:833446/1‑150 (MQ=255)
                              cAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAg                                                                                                                       >  2:1057328/1‑150 (MQ=255)
                                aaaCAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGAGTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGc                                                                                                                     <  2:914597/150‑1 (MQ=255)
                                 aaCAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGGCACTTATGGCGAGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCg                                                                                                                    <  2:11333/150‑1 (MQ=255)
                                 aaCAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCg                                                                                                                    <  2:833446/150‑1 (MQ=255)
                                  aCAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGa                                                                                                                   <  2:609525/150‑1 (MQ=255)
                                        tCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGc                                                                                                             <  1:1569169/150‑1 (MQ=255)
                                        tCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGc                                                                                                             <  1:435109/150‑1 (MQ=255)
                                           aCAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGcaca                                                                                                          >  2:740890/1‑150 (MQ=255)
                                           aCAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGcaca                                                                                                          <  2:35788/150‑1 (MQ=255)
                                            cAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAg                                                                                                         >  2:538943/1‑150 (MQ=255)
                                            cAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAg                                                                                                         >  2:1180267/1‑150 (MQ=255)
                                               aaGAGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTc                                                                                                      <  1:858296/150‑1 (MQ=255)
                                                aGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTcc                                                                                                     >  1:145294/1‑150 (MQ=255)
                                                   gccagccaGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAACACCGAAGGCGAGGATGCACAGTTCCAtc                                                                                                  >  2:1489123/1‑150 (MQ=255)
                                                   gccagccaGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCAtc                                                                                                  >  1:1114179/1‑150 (MQ=255)
                                                    ccagccagAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCAtcg                                                                                                 <  2:269831/150‑1 (MQ=255)
                                                       gccagAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCAtcgtcg                                                                                              >  2:1054891/1‑150 (MQ=255)
                                                          agAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCg                                                                                           >  1:1536987/1‑150 (MQ=255)
                                                           gAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCgg                                                                                          <  1:144757/150‑1 (MQ=255)
                                                            tacacCGGGATTGAACAGCCCAGTAAGCCGATGGCGTTGATCAACTTATCCTGCCACTTATGGCGCGTATTCCCGGCATTCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCggg                                                                                         <  2:685916/149‑1 (MQ=255)
                                                             acacCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGa                                                                                        <  1:749992/150‑1 (MQ=255)
                                                               acCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAg                                                                                      <  2:1361230/150‑1 (MQ=255)
                                                                ccGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGa                                                                                     >  1:1300042/1‑150 (MQ=255)
                                                                ccGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGa                                                                                     >  1:1155046/1‑150 (MQ=255)
                                                                   ggATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCcaac                                                                                                                                          >  1:1438853/1‑94 (MQ=255)
                                                                        gAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCtttt                                                                             >  2:1309777/1‑150 (MQ=255)
                                                                        gAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCtttt                                                                             >  2:759919/1‑150 (MQ=255)
                                                                         aaaaGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTa                                                                            <  1:1057328/150‑1 (MQ=255)
                                                                          ataGCCCAGTAAGGCGATGGCGTTGATCAAGTTATGCTGCCACTTATGGCGCGTAATCCCGGAAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCGATCCTCGCCGGGAAGACCTCTTTTaa                                                                           <  1:1316052/148‑1 (MQ=255)
                                                                           aaGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAAc                                                                          <  1:740890/150‑1 (MQ=255)
                                                                           aaGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAAc                                                                          <  1:336331/150‑1 (MQ=255)
                                                                            aGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACt                                                                         <  2:35692/150‑1 (MQ=255)
                                                                            aGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACt                                                                         >  1:1495181/1‑150 (MQ=255)
                                                                              cccAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGt                                                                       >  1:814293/1‑150 (MQ=255)
                                                                                cAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTc                                                                     <  1:1272541/150‑1 (MQ=255)
                                                                                  gTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCgg                                                                   >  2:918148/1‑150 (MQ=255)
                                                                                  gTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCgg                                                                   <  1:269103/150‑1 (MQ=255)
                                                                                   tAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGc                                                                  <  1:918148/150‑1 (MQ=255)
                                                                                    aaGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCa                                                                 <  1:1115852/150‑1 (MQ=255)
                                                                                         gATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCgg                                                            >  1:1156050/1‑150 (MQ=255)
                                                                                                 tGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGt                                                    >  1:1186644/1‑150 (MQ=255)
                                                                                                  gATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGtt                                                   <  2:279857/150‑1 (MQ=255)
                                                                                                  gATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGtt                                                   >  1:581950/1‑150 (MQ=255)
                                                                                                   aacAAGTTATCCTGCCACTTATTGCGCGTAATACCAGCAATCATACCAACCGGTATCCAAACAAACAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTa                                                  >  2:160710/3‑150 (MQ=255)
                                                                                                    tCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTaa                                                 >  2:653936/1‑150 (MQ=255)
                                                                                                     cAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAAt                                                >  2:751633/1‑150 (MQ=255)
                                                                                                        gTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAATGCt                                             <  2:1305084/150‑1 (MQ=255)
                                                                                                         ttATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAATGCTg                                            <  1:354769/150‑1 (MQ=255)
                                                                                                           aTCCTGCCACTTATGCCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAATGCTGGa                                          <  1:299699/150‑1 (MQ=255)
                                                                                                             ccTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAATGCTGGaaa                                        >  2:391784/1‑150 (MQ=255)
                                                                                                              cTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAATGCTGGAAac                                       <  2:1414628/150‑1 (MQ=255)
                                                                                                                        tGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAATGCTGGAAACACCAAAATcc                             <  1:759919/150‑1 (MQ=255)
                                                                                                                        tGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAATGCTGGAAACACCAAAATcc                             <  1:1495226/150‑1 (MQ=255)
                                                                                                                         ggCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAATGCTgg                                           >  1:172728/1‑135 (MQ=255)
                                                                                                                         ggCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAATGCTgg                                           <  2:172728/135‑1 (MQ=255)
                                                                                                                            gcgTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAATGCTGGAAACACCAAAATCCCAGt                         >  1:305954/1‑150 (MQ=255)
                                                                                                                                 aTCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAATAGCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAATGCTGGAAACACCAAAATCCCAGTGGATc                    <  2:504129/150‑1 (MQ=255)
                                                                                                                                   cccGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAATGCTGGAAACACCAAAATCCCAGTGGATCAg                  >  1:632156/1‑150 (MQ=255)
                                                                                                                                   cccGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAATGCTGGAAACACCAAAATCCCAGTGGATCAg                  >  2:739644/1‑150 (MQ=255)
                                                                                                                                     ttgttggcatcCCCGTGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAATGCTGGAAACACCAAAATCCCAGTGGATCAGGc                >  1:869825/8‑150 (MQ=255)
                                                                                                                                      tgttggcatcCCCGTGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTGTTAACTGTTCGGCAATCGGTTGACCGTTAATGCTGGAAACACCAAAATCCCAGTGGATCAGGcc               >  1:1151375/7‑150 (MQ=255)
                                                                                                                                      ggCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAATGCTGGAAACACCAAAATCCCAGTGGATCAGGcc               <  2:1114179/150‑1 (MQ=255)
                                                                                                                                        ttggcatcCCCGTGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAATGCTGGAAACACCAAAATCCCAGTGGATc                    >  1:388960/5‑143 (MQ=255)
                                                                                                                                        ttggcatcCCCGTGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAATGCTGGAAACACCAAAATCCCAGTGGATc                    <  2:388960/139‑1 (MQ=255)
                                                                                                                                        cAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTACCTGTTCGGCAATCGGTTGACCGTTAATGCTGGAAACACCAAAATCCCAGTGGATCAGGCCGt             >  2:1027110/1‑150 (MQ=255)
                                                                                                                                        cAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATAGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAATGCTGGAAACACCAAAATCCCAGTGGATCAGGCCGt             >  2:1281174/1‑150 (MQ=255)
                                                                                                                                         agccatcCCCGTGGCGATGACAACAATAACCGCAAAACCGAAGCCCAGGATGCACAGTTCCATCCTCCCCGGGAAGACCTCTTTTAACTGTTAGACAAAAGGTTCACCCTTAATGCACGAAAAACCAAAATACCAGACCATCACGCCGtt            >  2:951704/4‑150 (MQ=255)
                                                                                                                                          atcatcCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAATGCTGGAAACACCAAAATCCCAGTGGATCAGGCCGTTa           <  1:45753/150‑1 (MQ=255)
                                                                                                                                          atcatcCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGCGAAGACCTCTTTTAACTGGTCGGCAATCGGTTGACCGTTAATGCTGGAAACCCCAAAATCCCAGTGGATCAGGCCGTTa           >  1:1128824/1‑150 (MQ=255)
                                                                                                                                           tcatcCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTAGGCAATCGGTTGACCGTTAATGCTGGAAACACCAAAATCCCAGTGGATCAGGCCGTTaa          >  2:833782/1‑150 (MQ=255)
                                                                                                                                              tcCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAATGCTGGAAACACCAAAATCCCAGTGGATCAGGCCGTTAAAcc       <  2:1156050/150‑1 (MQ=255)
                                                                                                                                              tcCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAATGCTGGAAACACCAAAATCCCAGTGGATCAGGCCGTTAAAcc       <  1:719137/150‑1 (MQ=255)
                                                                                                                                              tcCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAATGCTGGAAACACCAAAATCCCAGTGGATCAGGCCGTTAAAcc       <  1:1054891/150‑1 (MQ=255)
                                                                                                                                                cccACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAATGCTGGAAACACCAAAATCCCAGTGGATCAGGCCGTTAAACCAg     <  1:1162068/150‑1 (MQ=255)
                                                                                                                                                cccACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAATGCTGGAAACACCAAAATCCCAGTGGATCAGGCCGTTAAACCAg     <  2:1300042/150‑1 (MQ=255)
                                                                                                                                                cccACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAATGCTGGAAACACCAAAATCCCAGTGGATCAGGCCGTTAAACCAg     >  2:910366/1‑150 (MQ=255)
                                                                                                                                                 ccACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAATGCTGGAAACACCAAAATCCCAGTGGATCAGGCCGTTAAACCAGa    <  1:484348/150‑1 (MQ=255)
                                                                                                                                                  cACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAATGCTGGAAACACCAAAATCCCAGTGGATCAGGCCGTTAAACCAGaa   >  2:342472/1‑150 (MQ=255)
                                                                                                                                                   aCGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCAGAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAATGCTGGAAACACGAAAATCGGAGTGGATCAGGCCGTTAAAGCTGAac  <  1:52528/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GAAACGGGCAACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAATGCTGGAAACACCAAAATCCCAGTGGATCAGGCCGTTAAACCAGAAC  >  NC_000913/1355040‑1355336

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