Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,303,984 A→C 6.6% A42A (GCT→GCG mdtN ← membrane fusion protein of efflux pump

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,303,9840AC 6.6% 260.8 / 2.5 99A42A (GCT→GCGmdtNmembrane fusion protein of efflux pump
Reads supporting (aligned to +/- strand):  ref base A (45/45);  new base C (6/2);  total (51/48)
Fisher's exact test for biased strand distribution p-value = 2.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTACTGCGAGGAGCTTTTTTCGGCGTACTTTCCATGAGCGCGTCTCAAATAGATTA  >  NC_000913/4303837‑4304130
                                                                                                                                                   |                                                                                                                                                  
aGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGc                                                                                                                                                  >  2:388034/1‑150 (MQ=255)
 gCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCg                                                                                                                                                 <  2:636398/150‑1 (MQ=255)
  cTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGt                                                                                                                                                >  2:1444729/1‑150 (MQ=255)
       gCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCTTACGCGTCagca                                                                                                                                           <  2:633290/150‑4 (MQ=255)
        ccAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAg                                                                                                                                          >  2:783422/1‑150 (MQ=255)
          aGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACCATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGGCATTAgat                                                                                                                                        >  2:634749/1‑148 (MQ=255)
          aGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACACTGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGtt                                                                                                                                        >  2:1025798/1‑150 (MQ=255)
          aGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGtt                                                                                                                                        >  1:755527/1‑150 (MQ=255)
          aGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGtt                                                                                                                                        >  2:942200/1‑150 (MQ=255)
          aGATTGGCTTCGTACGGGCGCGGGGCGATGCGGAACAGCAAAGCGCCCTGTTTGACTGCCTGGGTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGtt                                                                                                                                        <  1:1057827/150‑1 (MQ=255)
          aGATTGGCTGCGTACGGGCGCGGGTCGATGCGGAAGAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCTGCCGCTGACTTCCGGTACCACATCAATGGTATATGCTGACGCGTAAGCGTCATTAGtt                                                                                                                                        <  1:739411/150‑1 (MQ=255)
               ggCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATgg                                                                                                                                   >  2:765858/1‑150 (MQ=255)
                    cGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACCTCAATGGTATCTGCTGACGCGTAAGCTTCATTAGTTGATGGCGCAc                                                                                                                              >  2:306151/1‑150 (MQ=255)
                        cGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTc                                                                                                                          >  2:1565521/1‑150 (MQ=255)
                              cgGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGACTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTTAGCGTCATTAGTTGATGGCGCACTGTCTAcgcg                                                                                                                    >  2:870715/1‑150 (MQ=255)
                                ggTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGcc                                                                                                                  <  1:1177021/150‑1 (MQ=255)
                                ggTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGGCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGcc                                                                                                                  <  2:1527662/150‑1 (MQ=255)
                                 ggcGATGCGGAACAGCAAAACGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCAATAGTTGATGGCGCACTGTCTACGCGCCa                                                                                                                 <  2:837651/148‑1 (MQ=255)
                                 gTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCa                                                                                                                 <  1:1007035/150‑1 (MQ=255)
                                 gTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCa                                                                                                                 <  1:742771/150‑1 (MQ=255)
                                  tCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAg                                                                                                                <  2:1716000/150‑1 (MQ=255)
                                  gcGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGGGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGGATCTGCTGACGCGTACGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAg                                                                                                                <  1:1249314/149‑1 (MQ=255)
                                    gATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTTATGGCGCACTGTCTACGCGCCAGAt                                                                                                              >  1:197054/1‑150 (MQ=255)
                                    gATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGAt                                                                                                              >  1:618979/1‑150 (MQ=255)
                                     aTGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATa                                                                                                             >  2:872965/1‑150 (MQ=255)
                                        cGGAACAGCAACTCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCAGGTACCACATAAATGGTATCTGCTGACGCGTAAGCGTCATATGTTGATGGAGCACTTTCTACGCGACAGAtcaag                                                                                                          >  1:115186/1‑146 (MQ=255)
                                        cGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACg                                                                                                          <  2:1584566/150‑1 (MQ=255)
                                        cGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACg                                                                                                          <  2:1518881/150‑1 (MQ=255)
                                         ggAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGa                                                                                                         <  2:690164/150‑1 (MQ=255)
                                         ggAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGa                                                                                                         <  2:1686129/150‑1 (MQ=255)
                                          gAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGGGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGaa                                                                                                        <  2:1431048/150‑1 (MQ=255)
                                             cAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAAc                                                                                                     >  2:1599774/1‑150 (MQ=255)
                                                  aaTCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAggg                                                                                                <  2:550957/150‑1 (MQ=255)
                                                  aaTCGCCCTGTTTGACTGCCTGATTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCTTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAggg                                                                                                >  2:1503117/1‑150 (MQ=255)
                                                   aTCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACAGTCTACGCGCCAGATAACGAAAACAAGGGc                                                                                               >  1:788507/1‑150 (MQ=255)
                                                      gCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACATTGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAAc                                                                                            <  1:838439/150‑1 (MQ=255)
                                                      gCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAAc                                                                                            <  2:1451190/150‑1 (MQ=255)
                                                       cccTGTTTGAATGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTACGGTACAAAATCAATGGAATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAAAGAAAACAAGGGCTAcc                                                                                           >  1:1246935/1‑150 (MQ=255)
                                                                 cTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACg                                                                                 >  1:1118450/1‑150 (MQ=255)
                                                                   gCCTGGTTGTCGGTGACCGCCTGTTCTACAATGCGGCCGCTGACTTCCGTTACCACTTCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCa                                                                               <  1:1342235/150‑1 (MQ=255)
                                                                     cTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAAc                                                                             <  1:1157733/150‑1 (MQ=255)
                                                                        gTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCAc                                                                          <  1:941359/150‑1 (MQ=255)
                                                                          tGTCGGTGACCGCCAGATCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGAATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTa                                                                        <  1:893945/150‑1 (MQ=255)
                                                                              ggTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCAGTAGTTGATGGCGCACTGTCTACGCGCCAGATATCGAAAACAAGGGCAACCAGCGCCATCGCACCCACTAACag                                                                    >  2:209448/1‑150 (MQ=255)
                                                                                 gACCGCCAGTTCTACAAGGCGGCCGCTGCCTTCCGGTACCTCATCAATGGTATCTGCTGACGCGTAACCGTCATTAGTTGATGGCGCCCTGTCTACGCGCCAGAGAACGAAAACAAGGGCACCCAGCGCCAACGCAACCACTAACAGAGc                                                                 >  2:883842/1‑150 (MQ=255)
                                                                                   ccGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATTCCGAAAACAAGGGCAACCAGCGCCACCGCAACCTCGAACAGAGCAg                                                               >  2:150698/1‑150 (MQ=255)
                                                                                     gCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGCTAACGAATACAAGGGCAACCCGCGCCACCGCACCCACTAACGGAGCAggg                                                             >  1:282383/1‑150 (MQ=255)
                                                                                     gCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAggg                                                             <  1:181563/150‑1 (MQ=255)
                                                                                       cAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTCCGCGCCAGATACCGAAACCCAGGGCCACCCGCGCCACCGCAACCACGAACAGAGCAGGGaa                                                           >  2:1563925/1‑150 (MQ=255)
                                                                                       cAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTAATTAGTTGATGGCGCACTGTCTACGAGCCAGATAACGAAAACTAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGaa                                                           <  1:1467635/150‑1 (MQ=255)
                                                                                         gTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAAtt                                                         <  2:755527/150‑1 (MQ=255)
                                                                                          ttCTACATTGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGGAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAAttt                                                        <  1:870715/150‑1 (MQ=255)
                                                                                            cTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTAc                                                      >  2:540144/1‑150 (MQ=255)
                                                                                            cTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTAc                                                      >  2:1383030/1‑150 (MQ=255)
                                                                                            cTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTAc                                                      >  1:484369/1‑150 (MQ=255)
                                                                                            cTACAATGCGGCCACTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGTATTTAc                                                      >  2:764861/1‑150 (MQ=255)
                                                                                              aCAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTACTg                                                    >  1:355662/1‑150 (MQ=255)
                                                                                                 aTGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTACTGCgt                                                 >  1:1008676/1‑149 (MQ=255)
                                                                                                   gCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTTATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTACTGCgagg                                               <  2:197054/150‑1 (MQ=255)
                                                                                                   gCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTACTGCgagg                                               <  1:1388224/150‑1 (MQ=255)
                                                                                                    cGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTACTGCgagga                                              >  2:161140/1‑150 (MQ=255)
                                                                                                       ccGCTGACTTCCGGTACCACATCAATGGTATCTGCTTACGCGTAAGCGTCATTAGTTGATGGCGCACTTTCAACGCGCCAGATAACGAAAACAAGGGCAACCATAGCCAACGCAACCACTAACAGAGCAGGGAATTTACTGCGAGGAGCt                                           >  2:623681/1‑150 (MQ=255)
                                                                                                          cTGCCTTCCGGTACCACATCAATGGTATCTGCTGACGCGGAAGCGTCATTAGTTGATGGCGCACGGGCTACGCGCCAGATAACGGACTCAAGGGCGCCCAGCGCCACCGGCACGACTATGAGAGCAGGGAATTTACGGCGTGGAGCtggt                                        >  2:1064014/1‑147 (MQ=255)
                                                                                                            gACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTACTGCGAGGAGCtttttt                                      >  2:1362058/1‑150 (MQ=255)
                                                                                                             aCTTCCGGTACCACATCAATGGTATCTGCTTACGCGTAAGAGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTACTGCGAGGAGCTTTTTTc                                     >  1:1020433/1‑150 (MQ=255)
                                                                                                              cTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATaa                                                                                                            <  2:1394546/78‑1 (MQ=255)
                                                                                                              cTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATaa                                                                                                            >  1:1387784/1‑78 (MQ=255)
                                                                                                              cTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATaa                                                                                                            >  1:1394546/1‑78 (MQ=255)
                                                                                                              cTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATaa                                                                                                            <  2:1387784/78‑1 (MQ=255)
                                                                                                               ttCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTACTGCGAGGAGCTTTTTTCgg                                   <  2:363459/150‑1 (MQ=255)
                                                                                                                tCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCCACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTACTGCGAGGAGCTTTTTTCGGc                                  >  1:136704/1‑150 (MQ=255)
                                                                                                                 ccGGTACCACATCAATGGTATCTTCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTACTGCGAGGAGCTTTTTTCGGCg                                 <  1:1382695/150‑1 (MQ=255)
                                                                                                                 ccGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTACTGCGAGGAGCTTTTTTCGGCg                                 <  1:1127833/150‑1 (MQ=255)
                                                                                                                     tACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTACTGCGAGGAGCTTTTTTCGGCGTACt                             <  1:942200/150‑1 (MQ=255)
                                                                                                                       ccACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTACTGCGAGGAGCTTTTTTCGGCGTACttt                           >  1:1057313/1‑150 (MQ=255)
                                                                                                                        cacaTCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTACTGCGAGGAGCTTTTTTc                                     <  1:579678/139‑1 (MQ=255)
                                                                                                                        cacaTCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTACTGCGAGGAGCTTTTTTc                                     >  2:579678/1‑139 (MQ=255)
                                                                                                                        cacaTCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTACTGCGAGGAGCTTTTTTCGGCGTACTTTc                          >  2:142710/1‑150 (MQ=255)
                                                                                                                          caTCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTACTGCGAGGAGCTTTTTTCGGCGTACTTTCCa                        >  2:336753/1‑150 (MQ=255)
                                                                                                                           aTCAATGGTATCTGCTGACGCGTATGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAACCAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTACTGCGAGGAGCTTTTTTCGGCGTACTTTCCAt                       <  2:1317849/150‑1 (MQ=255)
                                                                                                                              aaTGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTACTGCGAGGAGCTTTTTTCGGCGTACTTTCCATGAg                    >  2:801991/1‑150 (MQ=255)
                                                                                                                              aaTGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTACTGCGAGGAGCTTTTTTCGGCGTACTTTCCATGAg                    <  1:1444729/150‑1 (MQ=255)
                                                                                                                                tGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTACTGCGAGGAGCTTTTTTCGGCGTACTTTCCATGAgcg                  <  1:945804/150‑1 (MQ=255)
                                                                                                                                    aactaatgACGCTTACGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTACTGCGAGGAGCTTTTTTCGGCGTAGTTTCCATGAGCGAGtc              >  1:1530436/7‑150 (MQ=255)
                                                                                                                                    aTCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTACTGCGAGGAGCTTTTTTCGGTGTACTTTCCATGAGCGCGtc              <  1:1364713/150‑1 (MQ=255)
                                                                                                                                     actaatgACGCTTACGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTACTGCGAGGAGCTTTTTTCGGCGTACTTTCCATGAGCGCGtct             >  1:229152/6‑150 (MQ=255)
                                                                                                                                      ctaatgACGCTTACGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAATCAAGGGCAACCAGCGCCAGCGCCCCCACTAACAGAGCAGGGAATTTACTGCGAGTAGCTGTTTTCGGCGTACTTTCAATGAGCGCGtctc            >  2:978170/5‑150 (MQ=255)
                                                                                                                                       tgctgACGCGTAAGCGTCATTAGTTGTTGGCGCACTGGCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTACTGCGAGGAGCTTTTTTCGGCGTACTTTCCATGAGCGCGTCTCa           <  2:520141/150‑1 (MQ=255)
                                                                                                                                       tgctgACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTACTGCGAGGAGCTTTTTTCGGCGTACTTTCCATGAGCGCGTCTCa           <  1:388034/150‑1 (MQ=255)
                                                                                                                                       tgctgACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTACTGCGAGGAGCTTTTTTCGGCGTACTTTCCATGAGCGCGTCTCa           <  2:484369/150‑1 (MQ=255)
                                                                                                                                       tgctgACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTACTGCGAGGAGCTTTTTTCGGCGTACTTTCCATGAGCGCGTCTCa           <  1:1421913/150‑1 (MQ=255)
                                                                                                                                       tgctgACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTACTGCGAGGAGCTTTTTTCGGCGTACTTTCCATGAGCGCGTCTCa           <  2:525957/150‑1 (MQ=255)
                                                                                                                                        aatgACGCTTACGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTACTGCGAGGAGCTTTTTTCGGCGTACTTTCCATGAGCGCGTCTCaa          >  1:1155458/3‑150 (MQ=255)
                                                                                                                                        aatgACGCTTACGCGTCATTAGTTGATGGCGCACTGTCTAAGCGCCAGATAAAAAAAACAAGTTCAACCAGCGCCAACGCAAAAACTAAAATATCAGGGAATTTACTGCGAGGAGCTTTTTTCGGCGTACTTTCCATGAGCGCGTCTCaa          >  2:218924/3‑150 (MQ=255)
                                                                                                                                        aatgACGCTTACGCGTCATTAGTTGATGGCGCACGGTCTACGCGCCAGATAACGAAAACAAGTGCAACCAGCGCCAACGCAACCACTAACAGAGCAGCGAATTTACTGCGAGGAGCTTTTATCGGCGTACTTTCCATGAGCGCGTCTCaa          >  1:286147/3‑150 (MQ=255)
                                                                                                                                         ctgACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACCCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTACTGCGAGGAGCTTTTTTCGGCGTACTTTCCATGAGCGCGTCTCaaa         >  2:1533747/1‑150 (MQ=255)
                                                                                                                                          tgACGCGTAAGCGTGATTATTTGATGGCGCAGTGGCTACGCGGCTGATAACGAAAACAAGGTCAACCAGCTCCATCGCAACCACTAACAGTGCTGGGAATTTTCTGCGAGTAGCTTTTTTCGGCGTACTTTCCATGAGCGCGTCTCAAAt        <  1:1011996/150‑1 (MQ=255)
                                                                                                                                              gcgTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGTGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTTCTGCGATGAGCTTTTTTCGGCGTACTTTCCATGAGCGCGTCTCAAATAGAt    <  1:960520/150‑1 (MQ=255)
                                                                                                                                                gTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTACTGCGAGGAGCTTTTTTCGGCGTACTTTCCATGAGCGCGTCTCAAATAGATTa  <  1:1424944/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATGGCGCACTGTCTACGCGCCAGATAACGAAAACAAGGGCAACCAGCGCCAACGCAACCACTAACAGAGCAGGGAATTTACTGCGAGGAGCTTTTTTCGGCGTACTTTCCATGAGCGCGTCTCAAATAGATTA  >  NC_000913/4303837‑4304130

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