Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 979,099 A→T 5.5% E140V (GAA→GTA)  U712_04785 → Stage V sporulation protein R

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP006881979,0990AT 5.5% 275.7 / 4.4 95E140V (GAA→GTA) U712_04785Stage V sporulation protein R
Reads supporting (aligned to +/- strand):  ref base A (40/47);  new base T (2/4);  total (43/52)
Fisher's exact test for biased strand distribution p-value = 6.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.62e-01

ACTGATCCAAAATAAATTGATTGTCGCACACGTTCTTGCTCATTGCGATTTCTTCAAAAACAACTGCCGTTTCCAAAACACGAAGCGCGATATGGTGGAAAGCATGGCGGCAACAGCTGAACGGATAAAACACTATGAACAGGTTCACGGGATAAAAGAAGTCGAATCTTTTTTAGATGCAATCCTTTCGATTCAGGAGCATATCGATCCGTCGTTGGTCCGACCAAAGCTCCAATGGAGTGTGGATGATGAGGAAGAGGAGATTGAAGAGGC  >  CP006881/978962‑979234
                                                                                                                                         |                                                                                                                                       
aCTGATCCAAAATAAATTGATTGTCGCACACGTTCTTGCTCATTGCGATTTCTTCAAAAACAACTGCCGTTTCCAAACCACGAAGCGCGATATGGTGGAAAGCATGGCGGCAACAGCTGAACGGATAAAACACTATGaa                                                                                                                                        <  8:205780/139‑1 (MQ=255)
aCTGATCCAAAATAAATTGATTGTCGCACACGTTCTTGCTCATTGCGATTTCTTCAAAAACAACTGCCGTTTCCAAAACACGAAGCGCGATATGGTGGAAAGCATGGCGGCAACAGCTGAACGGATAAAACACTATGaa                                                                                                                                        <  7:244796/139‑1 (MQ=255)
    aTCCAAAATAAATTGATTGTCGCACACGTTCTTGCTCATTGCGATTTCTTCTAAAACAACTGCCGTTTCCAAAACACGAAGCGCGATATGGTGGAAAGCATGGCGGCACCAGCTGAACGGATAAAACACTATGAACAgg                                                                                                                                    >  5:432777/1‑139 (MQ=255)
    aTCCAAAATAAATTGATTGTCGCACACGTTCTTGCTCATTGCGATTTCTTCAAAAACAACTGCCGTTTCCAAAACACGAAGCGCGATATGTTGGAAAGCATGGCGGCAACAGCTGAACGGATAAAACACTATGAACAgg                                                                                                                                    >  6:277251/1‑139 (MQ=255)
    aTCCAAAATAAATTGATTGTCGCACACGTTCTTGCTCATTGCGATTTCTTCAAAAACAACTGCCGTTTCCAAAACACGAAGCGCGATATGGTGGAAAGCATGGCGGCAACAGCTGAACGGATAAAACACTATGAACAgg                                                                                                                                    >  4:260619/1‑139 (MQ=255)
        aaaaTAAATTGATTGTCGCACACGTTCTTGCTCATTGCGATTTCTTCAAAAACAACTGCCGTTTCCAAAACACGAAGCGCGATATGGTGGAAAGCATGGCGGCAACAGCTGAACGGATAAAACACTATGAACAGGTTCa                                                                                                                                <  3:368082/139‑1 (MQ=255)
        aaaaTAAATTGATTGTCGCACACGTTCTTGCTCATTGCGATTTCTTCAAAAACAACTGCCGTTTCCAAAACACGAAGCGCGATATGGTGGAAAGCATGGCGGCAACAGCTGAACGGATAAAACACTATGAACAGGTTCa                                                                                                                                <  3:220049/139‑1 (MQ=255)
         aaataaatTGATTGTCGCACACGTTCTTGCTCATTGCGATTTCTTCAAAAACAACTGCCGTTTCCAAAACACGAAGCGCGATATGGTGGAAAGCATGGCGGCAACAGCTGAACGGATAAAACACTATGAACAGGTTCAc                                                                                                                               <  1:10021/139‑1 (MQ=255)
              aatTGATTGTCGCACACGTTCTTGCTCATTGCGATTTCTTCAAAAACACCTGCCGTTTCCAAAACACGAAGCGCGATATGGTGGAACGCATGGCGGCAACAGCTGCACGGATAAAACACTATGAACAGGTTCACGGGa                                                                                                                           >  7:43789/1‑138 (MQ=255)
              aatTGATTGTCGCACACGTTCTTGCTCATTGCGATTTCTTCAAAAACAACTGCCGTTTCCAAAACACGAAGCGCGATATGGTGGAAAGCATGGCGGCAACAGCTGAACGGATACAACACTATGAACAGGTTCACGGGa                                                                                                                           >  4:251424/1‑138 (MQ=255)
               attgattgTCGCACACGTTCTTGCTCATTGCGATTTCTTCAAAAACAACTGCCGTTTCCAAAACACGAAGCGCGATATGGTGGAAAGCATGGCGGCAACAGCTGAACGGATAAAACACTATGAACAGGTTCACGGGAt                                                                                                                          >  4:374136/1‑138 (MQ=255)
                ttgattgTCGCACACGTTCTTGCTCATTTCGATTTCTTCAAAAACAACTGCCGTTTCCAAAACACGAAGCGCGATATGGTGGAAAGCATGGCGGCAACAGCTGAACGGATAAAACACTATGAACAGGTTCACGGGATaa                                                                                                                        <  3:393017/139‑1 (MQ=255)
                   attgTCGCACACGTTCTTGCTCATTGCGATTTCTTCAAAAACAACTGCCGTTTCCAAAACACGAAGCGCGATATGGTGGAAAGCATGGCGGCAACAGCTGAACGGATAAAACACTATGAACAGGTTCACGGGATaa                                                                                                                        >  5:137550/1‑136 (MQ=255)
                   attgTCGCACACGTTCTTGCTCATTGCGATTTCTTCAAAAACAACTGCCGTTTCCAAAACACGAAGCGCGATATGGTGGAAAGCATGGCGGCAACAGCTGAACGGATAAAACACTATGAACAGGTTCACGGGATaa                                                                                                                        >  4:393017/1‑136 (MQ=255)
                      gTCGCACACGTTCTTGCTCATTGCGATTTCTTCAAAAACAACTGCCGTTTCCAAAACACGAAGCGCGATATGGTGGAAAGCATGGCGGCAACAGCTGAACGGATAAAACACTATGAACAGGTTCACGGGATAAaagaag                                                                                                                  >  2:196806/1‑139 (MQ=255)
                                     gCTCATTGCGATTTCTTCAAAAACAACTGCCGTTTCCAAAACACGAAGCGCGATATGGTGGAAAGCATGGCGGCAACAGCTGAACGGATAAAACACTATGAACAGGTTCACGGGATAAAAGAAGTCGAATCtttttt                                                                                                     <  3:57239/137‑1 (MQ=255)
                                         aTTTCGATTTCTTCAAAAACAACTGCCGTTTCCAAAACACGAAGCGCGATATGGTGGAAAGCATGGCGGCAACAGCTGAACGGATAAAACACTATGaa                                                                                                                                        <  6:206958/98‑1 (MQ=255)
                                         aTTGCGATTTCTTCAAAAACAACTGCCGTTTCCAAAACACGAAGCGCGATATGGTGGAAAGCATGGCGGCAACAGCTGAACGGATAAAACACTATGaa                                                                                                                                        >  5:206958/1‑98 (MQ=255)
                                         aTTGCGATTTCTTCAAAAACAACTGCCGTTTCCAAAACACGAAGCGCGATATGGTGGAAAGCATGGCGGCAACAGCTGAACGGATAAAACACTATGAACAGGTTCACGGGATAAAAGAAGTCGAATCTTTTTTAGATg                                                                                                >  4:358757/1‑138 (MQ=255)
                                         aTTGCGATTTCTTCAAAAACAACTGCCGTTTCCAAAAAACGAAGCGCGATATGGTGGAAAGCATGGCGGCAACAGCTGAACGGATAAAACACTATGCACAGGTTCACGGGATAAAAGAAGTCGAATCTTTTTTAGATg                                                                                                >  2:218815/1‑138 (MQ=255)
                                         aTTGCGATATCTTCAAAAACAACTGCCGTTTCCAAAACACGAAGCGCGATATGGTGGAAAGCATGGCGGCAACAGCTGAACGGATAAAACACTATGAACAGGTTCACGGGATAAAAGAAGTCGAATCTTTTTTAGATg                                                                                                >  1:231969/1‑138 (MQ=255)
                                           tGCGATTTCTTCAAAAACAACTGCCGTTTCCAAAACACGAAGCGCGATATGGTGGAAAGCATGGCGGCAACAGCTGAACGGATAAAACACTATGAACAGGTTCACGGGATAAAAGAAGTCGAATCtttttt                                                                                                     <  4:321232/131‑1 (MQ=255)
                                           tGCGATTTCTTCAAAAACAACTGCCGTTTCCAAAACACGAAGCGCGATATGGTGGAAAGCATGGCGGCAACAGCTGAACGGATAAAACACTATGAACAGGTTCACGGGATAAAAGAAGTCGAATCtttttt                                                                                                     >  3:321232/1‑131 (MQ=255)
                                           tGCGATTTCTTCAAAAACAACTGCCGTTTCCAAAACACGAAGCGCGATATGGTGGAAAGCATGGCGGCAACAGCTGAACGGATAAAACACTATGAACAGGTTCACGGGATAAAAGAAGTCGAATCTTTTTTAGATGCaa                                                                                             >  5:386647/1‑139 (MQ=255)
                                             cGATTTCATCAAAAACAACTGCCGTTTCCTAAACACGAAGCGAGATATGGTGGAAAGCATGGCGGCAACAGCTGAACGGATAAAACACTATTTACAGGTTCACGGGATAAAAGAAGTAGAATCTTTTTTAGATGCAAt                                                                                            <  6:137550/138‑1 (MQ=255)
                                               aTTTCTTCAAAAACAACTGCCGTTTCCAAAACACGAAGCGCGATATGGTGGAAAGCATGGCGGCAACAGCTGAACGGATAAAACACTATGAACAGGTTCACGGGATAAAAGAAGTCGAATCTTTTTTAGATGCAATCCt                                                                                         <  4:226233/139‑1 (MQ=255)
                                                  tcttcAAAATCAACTGCCGATTCAAAAACACGTACCGAGATAAGCTGGAAAGCATGGCTGCAAGAGCTGAACGGATAAAACACTATCTACTTGTTCACGGGATAAATGAAGTCGAATCTTTTTTAGTTGCAATCCttt                                                                                       <  6:264146/138‑1 (MQ=255)
                                                  tcttcAAAAACAACTGCCGTTTCCAAAACACGAAGCGCGATATGGTGGAAAGCATGGCGGCAACAGCTGAACGGATAAAACACTATGAACAGGTTCACGGGATAAAAGAAGTCGAATCTTTTTTAGATGCAATCCttt                                                                                       <  2:57012/138‑1 (MQ=255)
                                                  tcttcAAAAACAACTGCCGTTTCCAAAACACGAAGCGCGATATGGTGGAAAGCATGGCGGCAACAGCTGAACGGATAAAACACTATGAACAGGTTCACGGGATAAAAGAAGTCGAATCTTTTTTAGATGCAATCCttt                                                                                       <  1:36123/138‑1 (MQ=255)
                                                  tcttcAAAAACAACTGCCGTTTCCAAAACACGAAGCGCGATATGGTGGAAAGCATGGCGGCAACAGCTGAACGGATAAAACACTATGAACAGGTTCACGGGATAAAAGAAGTCGAATCTTTTTTAGATGCAATCCttt                                                                                       <  3:374136/138‑1 (MQ=255)
                                                  tcttcAAAAACAACTGCCGTTTCCAAAACACGAAGCGCGATATGGTGGAAAGCATGGCGGCAACAGCTGAACGGATAAAACACTATGAACAGGTTCACGGGATAAAAGAAGTCGAATCTTTTTTAGATGCAATCCttt                                                                                       <  2:131087/138‑1 (MQ=255)
                                                       aaaaaCAACTGCCTTTTCCAAAACACGTAGCGCTATATGGTTGAAAGCATGGCGGCAACAGCTGAACGGATAAAACACTATGAACAGGTTCACGGTATAAAAGAAGTCGAATCTTTTTTAGATGCAATCCTTTCGAtt                                                                                  <  7:229503/138‑1 (MQ=255)
                                                       aaaaaCAACTGCCGTTTCCAAAACACGAAGCGCGATATGGTGGAAAGCATGGCGGCAACAGCTGAACGGATAAAACACTATGAACAGGTTCACGGGATAAAAGAAGTCGAATCTTTTTTAGATGCAATCCTTTCGAtt                                                                                  <  1:24770/138‑1 (MQ=255)
                                                         aaaCAACTGCCTGTTCACAACCACGAAGCGCGATATGGTGGAAAGCATGGCGGCAACAGCTGAACGGATAAAACACTATTAACAGGTTCACGGGATAAAAGAAGTCGAATCTTTTTTAGATGCAATCCTTTCGATTCAg                                                                               <  4:227853/139‑1 (MQ=255)
                                                         aaaCAACTGCCGTTTCCAAAACACGAAGCGCGATATGGTGGAAAGCATGGCGGCAACAGCTGAACGGATAAAACACTATGAACAGGTTCACGGGATAAAAGAAGTCGAATCTTTTTTAGATGCAATCCTTTCGATTCAg                                                                               >  6:61777/1‑139 (MQ=255)
                                                           acaacTGCCGTCTCCAAAACACGAAGCGCGATATGGTGGAACGCAGGGCGGCAACAGCTGAACGGATAAAACACTATGAACAGGTTCACGGGATAAAAGAAGTCGAATCTTTTTAAGATGCAATCCTTTCGATTCAGGa                                                                             <  2:281877/139‑1 (MQ=255)
                                                                                        gATATGGTGGAAAGCATGGCGGCAACATCTGAACGGATAAAACACTATGAACAGGTTCACGGGATAAAAGAAGTCGAATCTTTTTTAGATGCAATCCTTTCGATTCAGGAGCATATCGATCCGTCGTTGGTCCGACCaa                                                <  3:169692/139‑1 (MQ=255)
                                                                                        gATATGGTGGAAAGCATGGCGGCAACAGCTGAACGGATAAAACACTATGAACAGGTTCACGGGATAAAAGAAGTCGAATCTTTTTTAGATGCAATCCTTTCGATTCAGGAGCATATCGATCCGTCGTTGGTCCGACCaa                                                <  1:196806/139‑1 (MQ=255)
                                                                                        gATATGGTGGAAAGCATGGCGGCAACAGCTGAACGGATAAAACACTATGAACAGGTTCACGGGATAAAAGAAGTCGAATCTTTTTTAGATGCAATCCTTTCGATTCAGGAGCATATCGATCCGTCGTTGGTCCGACCaa                                                <  3:186397/139‑1 (MQ=255)
                                                                                         atatGGTGGAAAGCATGGCGGCAACAGCTGAACGGATAAAACACTATGAACAGGTTCACGGGATAAaag                                                                                                                     >  3:446238/1‑69 (MQ=255)
                                                                                         atatGGTGGAAAGCATGGCGGCAACAGCTGAACGGATAAAACACTATGAACAGGTTCACGGGATAAaag                                                                                                                     <  4:446238/69‑1 (MQ=255)
                                                                                          tatGGTGGAAAGCATGGCGGCAACAGCTGAACGGATAAAACACTATGAACAGGTTCACTGGATAAAAGAAGGCGACTCTTTTTTAGATGGAATCCTTTGGATTCAGGAGCATATCGATCCGTCGTTGGGCCGACCAAAg                                              <  6:332893/139‑1 (MQ=255)
                                                                                          tatGGTGGAAAGCATGGCGGCAACAGCTGAACGGATAAAACACTATGAACAGGTTCACGGGATAAAAGAAGTCGAATCTTTTTTAGATGCAATCCTTTCGATTCAGGAGCATATCGATCCGTCGTTGGTCCGACCAAAg                                              <  5:61777/139‑1 (MQ=255)
                                                                                          tatGGTGGAAAGCATGGCGGCAACAGCTGAACGGATAAAACACTATGAACAGGTTCACGGGATAAAAGAAGTCGAATCTTTTTTAGATGCAATCCTTTCGATTCAGGAGCATATCGATCCGTCGTTGGTCCGACCAAAg                                              <  7:380139/139‑1 (MQ=255)
                                                                                              gttgAAAGCATGGCGGCAACAGCTGAACGGATAAAACACTCTGAACAGGTTCACGGGATAAAAGAAGTCGAATCTTTTTTAGATGCAATCCTTTCGATTCAGGAGCATATCGATCCGTCGTTGGTCCGACCAAAGCTc                                           <  5:277251/135‑1 (MQ=255)
                                                                                               tggAAAGCATGGCGGCAACAGCTGAACGGATAAAACACTATGAACAGGTTCa                                                                                                                                <  1:201161/52‑1 (MQ=255)
                                                                                               tggAAAGCATGGCGGCAACAGCTGAACGGATAAAACACTATGAACAGGTTCa                                                                                                                                >  2:201161/1‑52 (MQ=255)
                                                                                               tggAAAGCATGGCGGCAACAGCTGAACGGATAAAACACTATGAACAGGTTCACGGGATAAAAGAAGTCGAATCTTTTTTAGATGCAATCCTTTCGATTCAGGAGCATATCGATCCGTCGTTGGTCCGACCAAAGCTcc                                          >  8:130697/1‑138 (MQ=255)
                                                                                               tggAAAGCATGGCGGCAACAGCTGAACGGATAAAACACTATGAACAGGTTCACGGGATAAAAGAAGTCGAATCTTTTTTAGATGCAATCCTTTCGATTCAGGAGCATATCGATCCGTCGTTGGTCCGACCAAAGCTcc                                          >  8:136965/1‑138 (MQ=255)
                                                                                                  accgCATGGCGGCAACAGCTGAACGGATAAAACACTATGACCAGGTTCACGGGATAAACGAAGTCGAATCTTTTTTAGATGCAATCCTTTCGATTCAGGATCATATCGATCCGTCGTTGGTCCGACCAAAGCTCCAATg                                      <  4:164724/136‑1 (MQ=255)
                                                                                                  aaaGCATGGCGGCAACAGCTGAACGGATAAAACACTATGAACAGGTTCACGGGATAAAAGAAGTCGAATCTTTTTTAGATGCAATCCTTTCGATTCAGGAGCATATCGATCCGTCGTTGGTCCGACCAAAGCTCCAATg                                      >  7:270723/1‑139 (MQ=255)
                                                                                                  aaaGCATGGCGGCAACAGCTGAACGGATAAAACACTATGAACAGGTTCACGGGATAAAAGAAGTCGAATCTTTTTTAGATGCAATCCTTTCGATTCAGGAGCATATCGATCCGTCGTTGGTCCGACCAAAGCTCCAATg                                      <  7:130697/139‑1 (MQ=255)
                                                                                                  aaaGCATGGCGGCAACAGCTGAACGGATAAAACACTATGAACAGGTTCACGGGATAAAAGAAGTCGAATCTTTTTTAGATGCAATCCTTTCGATTCAGGAGCATATCGATCCGTCGTTGGTCCGACCAAAGCTCCAATg                                      <  7:136965/139‑1 (MQ=255)
                                                                                                    aGCATGGCGGCAACAGCTGAACGGATAAAACACTATGAACAGGTTCACGGGATAAAAGAAGTCGAATCTTTTTTAGATGCAATCCTTTCGATTCAGGAGCATATCGATCCGTCGTTGGTCCGACCAAAGCTCCAATg                                      >  7:291241/1‑137 (MQ=255)
                                                                                                    aGCATGGCGGCAACAGCTGAACGGATAAAACACTATGAACAGGTTCACGGGATAAAAGAAGTCGAATCTTTTTTAGATGCAATCCTTTCGATTCAGGAGCATATCGATCCGTCGTTGGTCCGACCAAAGCTCCAATg                                      >  4:432901/1‑137 (MQ=255)
                                                                                                     gCATGGCGGCAACAGCTGCACGGATAAAACACTATGAACAGGTTCACGGGATAAAAGAAGTCGAATCTTTTTTAGATGCAATCCTTTCGATTCAGGAGCATATCGATCCGTCGTTGGTCCGACCAAAGCTCCAATGGAg                                   <  4:240950/139‑1 (MQ=255)
                                                                                                     gCATGCCGGCAACAGATGAACGGATAACACACTCTGAACAGTGTCACTGGATAAAAGACGTCGAATCTTTTTTAGATGCAATCCTTTCGATTCAGGAGCATATCGATCCGTCGTTGGTCCGACCAAAGCTCCAATGGAg                                   <  2:415265/139‑1 (MQ=255)
                                                                                                      cATGGCGGCAACAGCTGAACGGATAAAACACTATGAACAGGTTCACGGGATAAAAGAAGTCGAATCTTTTTTAGATGCAAt                                                                                            >  8:252343/1‑81 (MQ=255)
                                                                                                      cATGGCGGCAACAGCTGAACGGATAAAACACTATGAACAGGTTCACGGGATAAAAGAAGTCGAATCTTTTTTAGATGCAAt                                                                                            <  7:252343/81‑1 (MQ=255)
                                                                                                      cATGGCGGCAACAGCTGAACGGATAAAACACTATGAACAGGTTCACGGGATAAAAGAAGTCGAATCTTTTTTAGATGCAATCCttt                                                                                       >  2:423641/1‑86 (MQ=255)
                                                                                                      cATGGCGGCAACAGCTGAACGGATAAAACACTATGAACAGGTTCACGGGATAAAAGAAGTCGAATCTTTTTTAGATGCAATCCttt                                                                                       <  1:423641/86‑1 (MQ=255)
                                                                                                        tGGCGGCAACAGCTGAACGGATAAAACACTATGAACAGGTTCACGGGATAAAAGAAGTCGAATCTTTTTTAGATGCAATCCTTTCGATTCAGGa                                                                             >  5:69544/1‑94 (MQ=255)
                                                                                                        tGGCGGCAACAGCTGAACGGATAAAACACTATGAACAGGTTCACGGGATAAAAGAAGTCGAATCTTTTTTAGATGCAATCCTTTCGATTCAGGa                                                                             <  6:69544/94‑1 (MQ=255)
                                                                                                         ggcggcAACAGCTGAACGGATAAAACACTATGAACAGGTTCACGGGATAAAAGAAGTCGAATCTTTTTTAGATGCAATCCTTTCGATTc                                                                                 >  1:52745/1‑89 (MQ=255)
                                                                                                         ggcggcAACAGCTGAACGGATAAAACACTATGAACAGGTTCACGGGATAAAAGAAGTCGAATCTTTTTTAGATGCAATCCTTTCGATTc                                                                                 <  2:52745/89‑1 (MQ=255)
                                                                                                         ggcggcAACAGCTGAACGGATAAAACACTATGAACAGGTTCACGGGATAAAAGAAGTCGAATCTTTTTTAGATGCAATCCTTTCGATTCAGGAGCATATCGATCCGTCGTTGGTCCGACCAAAGCTCCAATGGAGtgtg                               >  5:392217/1‑139 (MQ=255)
                                                                                                           cggcAACAGCTGTACGGATAAAACACAATGCCCAGGTTGACGGCAGAAAAGAAGCAGAATCTTTTTTAGATGCAATCCTTTCGATTCAGGAGCATATCGATCCGTCGTTGGTCCGACCAAAGCTCCAATGGAGTGTGga                             <  4:426763/139‑1 (MQ=255)
                                                                                                              cAACAGCTGAACGGATAAAACACTATGAACAGGTTCACGGGATAAAAGAAGTCGAATCTTTTTTAGATGCAATCCTTTCGATTCAGGAGCATATCGATc                                                                  <  2:45375/99‑1 (MQ=255)
                                                                                                              cAACAGCTGAACGGATAAAACACTATGAACAGGTTCACGGGATAAAAGAAGTCGAATCTTTTTTAGATGCAATCCTTTCGATTCAGGAGCATATCGATc                                                                  >  1:45375/1‑99 (MQ=255)
                                                                                                              cAACAGCTGAACGGATAAAACACTATGAACAGGTTCACGGGATAAAAGAAGGCGAATCTTTTTTAGATGCAATCCTTTCGATTCAGGAGCATATCGATCCGTCGTTGGTCCGACCAAAGCTCCAATGGAGTGTGgatga                          <  3:221373/139‑1 (MQ=255)
                                                                                                                  aGCTGAACTGATAAAACACTATGAACAGGTTCACGGGATAAAAGAAGTCGAATCTTTTTTAGATGCAATCCTTTCGATTCAGGAGCATATCGATCCGTCGTTGGTCCGACCAAAGCTCCAATGGAGTGTGGATGATGAg                      <  1:94468/139‑1 (MQ=255)
                                                                                                                  aGCTGAACGGATAAAACACTATGAACAGGTTCACGGGATAAAAGAAGTCGAATCTTTTTTAGATGCAATCCTTTCGATTCAGGAGCATATCGATCCGTCGTTGGTCCGACCAAAGCTCCAATGGAGTGTGGATGATGAg                      <  8:43789/139‑1 (MQ=255)
                                                                                                                       aaCGGATAAAACACTATGAACAGGTTCACGGGATAAAAGAAGTCGAATCTTTTTTAGATGCAATCCTTTCGATTCATGAGCAGATCGATCCGACGTTGGTCCGACCACAGCTCCCATGGAGAGTGTATGATGAGGaa                   >  4:46290/1‑137 (MQ=255)
                                                                                                                       aaCGGATAAAACACTATGAACAGGTTCACGGGATAAAAGAAGTCGAATCTTTTTTAGATGCAATCCTTTCGATTCAGGAGCATATCGATCCGTCGTTGGTCCGACCAAAGCTCCAATGGAGTGTGGATGATGAGGAaga                 <  8:270723/139‑1 (MQ=255)
                                                                                                                        ccgtaaaaaaCACTATGAACCGGTTCACGGGAAAAAAGAAGTCGAATGTTTTTTAGCTGCAATCCTTTCGATTCAGGAGCATATCGATCCGTCGTTGGTCCGACCAAAGCTCCAATGGAGTGTGGATGATGAGGAagag                <  4:390928/133‑1 (MQ=255)
                                                                                                                        cagtaTAAAACACTATTAAATGGTTCTCGGTATAAAATAAGTCGAATCTTTTTTATATGAAATCCTTTCGATTCAGGAGCATATCGATCCGTCGTTGGTCCGACCAAAGCTCCAATGGAGTGTGGATGATGAGGAaga                 <  4:332707/134‑1 (MQ=255)
                                                                                                                         cGGATAAAACACTATGAACAGGTTCACGGGATAAAAGAAGTCGAATCTTTTTTAGATGCAATCCTTTCGATTCAGGAGCATATCGATCCGTCGTTGGTCCGACCAAAGCTCCAATGGAGTGTGGATGATGAGGAagag                >  5:148480/1‑138 (MQ=255)
                                                                                                                         cGGATAAAACACTATGAACAGGTTCACGGGATAAAAGAAGTCGAATCTTTTTTAGATGCAATCCTTTCGATTCAGGAGCATATCGATCCGTCGTTGGTCCGACCAAAGCTCCAATGGAGTGTGGATGATGAGGAAgagg               >  3:379824/1‑139 (MQ=255)
                                                                                                                            aTAAAACACTATGAACAGGTTCACGGGATAAAAGAAGTCGAATCTTTTTTAGATGCAATCCTTTCGATTCAGGAGCATATCGATCCGTCGTTGGTCCGACCAAAGCTCCAATGGAGTGTGgatga                          <  6:108441/125‑1 (MQ=255)
                                                                                                                            aTAAAACACTATGAACAGGTTCACGGGATAAAAGAAGTCGAATCTTTTTTAGATGCAATCCTTTCGATTCAGGAGCATATCGATCCGTCGTTGGTCCGACCAAAGCTCCAATGGAGTGTGgatga                          >  5:108441/1‑125 (MQ=255)
                                                                                                                            aTAAAACACTATGAACAGGTTCACGGGATAAAAGAAGTCGAATCTTTTTTAGATGCAATCCTTTCGATTCAGGAGCATATCGATCCGTCGTTGGTCCGACCAAAGCTCCAATGGAGTGTGGATGATGAGGAag                  >  8:66562/1‑133 (MQ=255)
                                                                                                                            aTAAAACACTATGAACAGGTTCACGGGATAAAAGAAGTCGAATCTTTTTTAGATGCAATCCTTTCGATTCAGGAGCATATCGATCCGTCGTTGGTCCGACCAAAGCTCCAATGGAGTGTGGATGATGAGGAag                  <  7:66562/133‑1 (MQ=255)
                                                                                                                             tAAAACACTATTAACCGGCTCACGGTATAAAAGAAGTCGAATCTTTTTTAGATGCAATCCTTTCGATTCAGGAGCATATCGATCCGTCGTTGGTCCGACCAAAGCTCCAATGGAGTGTGGATGATGAGGAAGAGGAGAt           <  6:160819/139‑1 (MQ=255)
                                                                                                                                    cTATGACCAGGTTCACGGGATAAAAGCAGTCGAATCTTTTTTAGATGCAATCCTTTCGATTCAGGAGCATATCGATCCGTCGGGGGTCCGACCAAAGCTCCAATGGTGTGTGGATGATGAGGAAGAGGAGATTGAagag    <  7:302957/139‑1 (MQ=255)
                                                                                                                                     tATGAACAGGTTCACGGGATAAAAGAAGTCGAATCTTTTTTAGATGCAATCCTTTCGATTCAGGCGCATATCGATCCGTCGTTGGTCCGACCAAAGCTCCAATGGAGTGTGGATGATGAGGAAGAGGAGATTGAAGCgg   >  8:227021/1‑139 (MQ=255)
                                                                                                                                     tATGAACAGGTTCACGGGATAAAAGAAGTCGAATCTTTTTTAGATGCAATCCTTTCGATTCAGGAGCATATCGATCCTTCGTTGGTCCGACCAAAGCTCCAATGGAGTGTGGATGATGAGGAAGAGGAGATTGAAGAgg   >  5:5416/1‑139 (MQ=255)
                                                                                                                                      attaaCAGGATCACGGGAAAAAAGAAGTCGAATCATTTTAAGATGCAATCCTATCCATACAGTAGCATATCGCTCCGTCGTTCGTCCGACCAAAGCTACAATGGAGTGTGGATTAGGAGGAAGAGGAGATTGAAGACGc  >  8:23230/4‑139 (MQ=255)
                                                                                                                                      aTGTACAGGTTCACGGGATAAAAGAAGTCGAATCTTTTTTAGATGCAATCCTTTCGATTCAGGAGCATATCGATCCGTCGTTGGTCCGACCAAAGCTCCAATGGAGTGTGGATGATGAGGAAGAGGAGAtt          >  4:54948/1‑131 (MQ=255)
                                                                                                                                      aTGTACAGGTTCACGGGATAAAAGAAGTCGAATCTTTTTTAGATGCAATCCTTTCGATTCAGGAGCATATCGATCCGTCGTTGGTCCGACCAAAGCTCCAATGGAGTGTGGATGATGAGGAAGAGGAGAtt          <  3:54948/131‑1 (MQ=255)
                                                                                                                                      aTGTACAGGTTCACGGGATAAAAGAAGTCGAATCTTTTTTAGATGCAATCCTTTCGATTCAGGAGCATATCGATCCGTCGTTGGTCCGACCAAAGCTCCAATGGAGTGTGGATGATGAGGAAGAGGAGAtt          <  3:60959/131‑1 (MQ=255)
                                                                                                                                      aTGTACAGGTTCACGGGATAAAAGAAGTCGAATCTTTTTTAGATGCAATCCTTTCGATTCAGGAGCATATCGATCCGTCGTTGGTCCGACCAAAGCTCCAATGGAGTGTGGATGATGAGGAAGAGGAGAtt          >  4:60959/1‑131 (MQ=255)
                                                                                                                                      aTGAACAGGTTCACGGGATAAAAGAAGTCGAATCTTTTTTAGATGCAATCCTTTCGATTCAGGAGCATATCGATCCGTCGTTGGTCCGACCaa                                                <  8:106107/93‑1 (MQ=255)
                                                                                                                                      aTGAACAGGTTCACGGGATAAAAGAAGTCGAATCTTTTTTAGATGCAATCCTTTCGATTCAGGAGCATATCGATCCGTCGTTGGTCCGACCaa                                                >  7:106107/1‑93 (MQ=255)
                                                                                                                                      aTGAACAGGTTCACGGGATAAAAGAAGTCGAATCTTTTTTAGATGCAATCCTTTCGATTCAGGAGCATATCGATCCGTCGTTGGTCCGACCaa                                                >  5:233209/1‑93 (MQ=255)
                                                                                                                                      aTGAACAGGTTCACGGGATAAAAGAAGTCGAATCTTTTTTAGATGCAATCCTTTCGATTCAGGAGCATATCGATCCGTCGTTGGTCCGACCaa                                                <  6:233209/93‑1 (MQ=255)
                                                                                                                                      aTGAACAGGTTCACGGGATAAAAGAAGTCGAATCTTTTTTAGATGCAATCCTTTCGATTCAGGAGCATATCGATCCGTCGTTGGTCCGACCAAAGCTCCAATGGAg                                   <  7:305321/106‑1 (MQ=255)
                                                                                                                                      aTGAACAGGTTCACGGGATAAAAGAAGTCGAATCTTTTTTAGATGCAATCCTTTCGATTCAGGAGCATATCGATCCGTCGTTGGTCCGACCAAAGCTCCAATGGAg                                   >  8:305321/1‑106 (MQ=255)
                                                                                                                                      aTGAACAGGTTCACGGGATAAAAGAAGTCGAATCTTTTTTAGATGCAATCCTTTCGATTCAGGAGCATATCGATCCGTCGTTGGTCCGACCAAAGCTCCAATGGAGTTTGGATGATGAGGAAGAGGAGATTGAAGAgg   >  1:232206/1‑138 (MQ=255)
                                                                                                                                      aTGAACAGGTTCACGGGATAAAAGAAGTCGAATCTTTTTTAGATGCAATCCTTTCGATTCAGGAGCATATCGATCCGTCGTTGGTCCGACCAAAGCTCCAATGGAGTGTGGATGATGAGGAAGAGGAGATTGTagctgt  >  3:332707/1‑135 (MQ=255)
                                                                                                                                      aTGAACAGGTTCACGGGATAAAAGAAGTCGAATCTTTTTTAGATGCAATCCTTTCGATTCAGGAGCATATCGATCCGTCGTTGGTCCGACCAAAGCTCCAATGGAGTGTGGATGATGAGGAAGAGGAGATTGAAGAgg   >  6:316259/1‑138 (MQ=255)
                                                                                                                                         |                                                                                                                                       
ACTGATCCAAAATAAATTGATTGTCGCACACGTTCTTGCTCATTGCGATTTCTTCAAAAACAACTGCCGTTTCCAAAACACGAAGCGCGATATGGTGGAAAGCATGGCGGCAACAGCTGAACGGATAAAACACTATGAACAGGTTCACGGGATAAAAGAAGTCGAATCTTTTTTAGATGCAATCCTTTCGATTCAGGAGCATATCGATCCGTCGTTGGTCCGACCAAAGCTCCAATGGAGTGTGGATGATGAGGAAGAGGAGATTGAAGAGGC  >  CP006881/978962‑979234

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: