Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,757,734 T→G 6.2% W264G (TGG→GGG)  PP_1563 → terminase large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,757,7340TG 6.2% 242.4 / 4.1 78W264G (TGG→GGG) PP_1563terminase large subunit
Reads supporting (aligned to +/- strand):  ref base T (32/41);  new base G (2/3);  total (34/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGACACGATGGAGACCGGGATGGGGGCCCGCGAGCAACCCATCATGCTGGTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGC‑GAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTGGGATATCCGTCGCCGAGGACTTCCTGCTTGCCCGGCAGCGTGATGCGATGCAGTCGGCGA  >  NC_002947/1757586‑1757872
                                                                                                                                                     |                                                                                                                                          
cGACACGATGGAGACCGGGATGGGGGCCCGCGAGCAACCCATTATGCTGGTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGC‑GAATGCTGGAGGGAGTGATCGACCGGCCGGACCTAt                                                                                                                                            <  1:115685/149‑1 (MQ=255)
 gACACGATGGAGACCGGGATGGGGGCCCGCGAGCAACCCATCATGCTGGTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGC‑GAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATg                                                                                                                                           <  6:287786/149‑1 (MQ=255)
  acacGATGGAGACCGGGATGGGGGCCCGCGAGCAACCCATCATGCTGGTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGC‑GAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATgg                                                                                                                                          <  6:29377/149‑1 (MQ=255)
  acacGATGGAGACCGGGATGGGGGCCCGCGAGCAACCCATCATGCTGGTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGC‑GAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATgg                                                                                                                                          >  1:71145/1‑149 (MQ=255)
   cacGATGGAGACCGGGATGGGGGCCCGCGAGCAACCCATCATGCTGGTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGC‑GAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATggc                                                                                                                                         >  2:72194/1‑149 (MQ=255)
     cGATGGAGACCGGGATGGGGGCCCGCGAGCAACCCATCATGCTGGTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGC‑GAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATggccg                                                                                                                                       <  1:54991/149‑1 (MQ=255)
      gATGGAGACCGGGATGGGGGCCCGCGAGCAACCCATCATGCTGGTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGC‑GAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATggccgg                                                                                                                                      <  3:341302/149‑1 (MQ=255)
      gATGGAGACCGGGATGGGGGCCCGCGAGCAACCCATCATGCTGGTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGC‑GAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATggccgg                                                                                                                                      >  1:336571/1‑149 (MQ=255)
        tGGAGACCGGGATGGGGGCCCGCGAGCAACCCATCATGCTGGTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGC‑GAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATggccggcc                                                                                                                                    >  5:39306/1‑149 (MQ=255)
         ggAGACCGGGATGGGGGCCCGCGAGCAACCCATCATGCTGGTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGC‑GAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATggccggcc                                                                                                                                    <  2:71145/148‑1 (MQ=255)
          gagaCCGGGATGGGGGCCCGCGAGCAACCCATCATGCTGGTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGC‑GAATGCTGGAGGGTGTGATCGACCGGCCGGCCATAGgtccggccgg                                                                                                                                  <  6:215120/149‑3 (MQ=255)
          gagaCCGGGATGGGGGCCCGCGAGCAACCCATCATGCTGGTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGC‑GAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATggccggccg                                                                                                                                   <  6:39306/148‑2 (MQ=255)
              ccGGGATGGGGGCCCGCGAGCAACCCATCATGCTGGTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGC‑GAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTac                                                                                                                              <  1:72194/149‑1 (MQ=255)
               cGGGATGGGGGCCCGCGAGCAACCCATCATGCTGGTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGC‑GAATGCTGGAGGGTGTGATCGACCGGCCGGCCATAGgtccggccggtc                                                                                                                                <  5:171491/146‑5 (MQ=255)
               cGGGATGGGGGCCCGCGAGCAACCCATCATGCTGGTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGC‑GAATGCTGGAGGGTGTGATCGACCGGCCGGCCATAGgtccggccggtc                                                                                                                                >  6:171491/1‑142 (MQ=255)
               cGGGATGGGGGCCCGCGAGCAACCCATCATGCTGGTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGC‑GAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTaca                                                                                                                             <  6:138779/149‑1 (MQ=255)
               cGGGATGGGGGCCCGCGAGCAACCCATCATGCTGGTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGC‑GAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTaca                                                                                                                             <  5:299639/149‑1 (MQ=255)
               cGGGATGGGGGCCCGCGAGCAACCCATCATGCTGGTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGC‑GAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTac                                                                                                                              >  6:299639/1‑148 (MQ=255)
                  gATGGGGGCCCGCGAGCAACCCATCATGCTGGTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGC‑GAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATggccggc                                                                                                                                     >  3:116457/1‑138 (MQ=255)
                  gATGGGGGCCCGCGAGCAACCCATCATGCTGGTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGC‑GAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATggccggc                                                                                                                                     <  4:116457/138‑1 (MQ=255)
                      ggggCCCGCGAGCAACCCATCATGCTGGTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGC‑GAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTACACTATCGa                                                                                                                      <  3:328279/149‑1 (MQ=255)
                       gggCCCGCGAGCAACCCATCATGCTGGTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGC‑GAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTACACTATCGAc                                                                                                                     <  3:193819/149‑1 (MQ=255)
                           ccGCGAGCAACCCATCATGCTGGTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGC‑GAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATCGCCGGCCCTCTACACGATCGCCCatc                                                                                                                 >  3:244439/1‑147 (MQ=255)
                            cgcgAGCAACCCATCATGCTGGTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGC‑GAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAggg                                                                                                                >  4:329942/1‑149 (MQ=255)
                                 gCAACCCATCATGCTGGTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGC‑GAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGc                                                                                                               <  7:43559/145‑1 (MQ=255)
                                 gCAACCCATCATGCTGGTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGC‑GAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGc                                                                                                               >  8:43559/1‑145 (MQ=255)
                                    aCCCATCATGCTGGGGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGC‑GAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGAtt                                                                                                        <  8:307236/149‑1 (MQ=255)
                                     cccATCATGCTGGTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGC‑GAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTg                                                                                                       >  4:180231/1‑149 (MQ=255)
                                      ccATCATGCTGGTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGC‑GAATGCTGGAGGGTGTAATCGACCGGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTgg                                                                                                      <  2:144888/149‑1 (MQ=255)
                                             gCTGGTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGC‑GAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGc                                                                                                >  3:255534/1‑148 (MQ=255)
                                               tGGTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGC‑GAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATggccggcc                                                                                                                                    <  7:104523/110‑1 (MQ=255)
                                               tGGTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGC‑GAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATggccggcc                                                                                                                                    >  8:104523/1‑110 (MQ=255)
                                                 gTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGC‑GAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGaaa                                                                                            <  2:336571/148‑1 (MQ=255)
                                                 gTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGC‑GAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGaa                                                                                             <  3:329942/147‑1 (MQ=255)
                                                        ccacGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGC‑GAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATggccggcc                                                                                                                                    <  6:328292/101‑1 (MQ=255)
                                                        ccacGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGC‑GAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATggccggcc                                                                                                                                    >  5:328292/1‑101 (MQ=255)
                                                        ccacGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGC‑GAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTg                                                                                    >  6:157694/1‑149 (MQ=255)
                                                              cggGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGC‑GAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAg                                                                              <  3:180231/149‑1 (MQ=255)
                                                                   tCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGCTGAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCg                                                                           >  1:252354/1‑148 (MQ=255)
                                                                   tCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGC‑GAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTACACTATCGACCa                                                                                                                   >  5:322298/1‑107 (MQ=255)
                                                                   tCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGC‑GAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTACACTATCGACCa                                                                                                                   <  6:322298/107‑1 (MQ=255)
                                                                     aaGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGC‑GAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGaa                                                                         <  5:157694/147‑1 (MQ=255)
                                                                       gTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGC‑GAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTACACTATCGACCa                                                                                                                   >  5:49183/1‑103 (MQ=255)
                                                                       gTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGC‑GAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTACACTATCGACCa                                                                                                                   <  6:49183/103‑1 (MQ=255)
                                                                          ttGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGC‑GAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCaaa                                                                   <  2:252354/148‑1 (MQ=255)
                                                                          ttGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGC‑GAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCaaa                                                                   <  4:255534/148‑1 (MQ=255)
                                                                             gcggcCCCTGCCACCAGTTGATCCGCGATTCCGAGC‑GAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACttt                                                               >  2:21093/1‑149 (MQ=255)
                                                                                gcCCCTGCCACCAGTTGATCCGCGATTCCGAGC‑GAATGCTGTAGGGTTTGATCGACCGGCCGGACCTATGGCCGGCCATCTACAGTATCGACCAGGGCGAAGATTGGACCAGCGTAATCTCAATGCGAAAGGGGAATCAAAACTTTggg                                                            >  4:108879/1‑149 (MQ=255)
                                                                                    cTGCCACCAGTTGATCCGCGATTCCGAGC‑GAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTGGGatat                                                        >  5:276525/1‑149 (MQ=255)
                                                                                        caccaGTTGATCCGCGATTCCGAGC‑GAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTg                                                              >  6:120361/1‑139 (MQ=255)
                                                                                        caccaGTTGATCCGCGATTCCGAGC‑GAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTg                                                              <  5:120361/139‑1 (MQ=255)
                                                                                         accaGTTGATCCGCGATTCCGAGC‑GAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTGGGATATCCGTc                                                   >  6:45989/1‑149 (MQ=255)
                                                                                            aGTTGATCCGCGATTCCGAGC‑GAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTg                                                              >  5:351117/1‑135 (MQ=255)
                                                                                            aGTTGATCCGCGATTCCGAGC‑GAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTg                                                              <  6:351117/135‑1 (MQ=255)
                                                                                             gTTGATCCGCGATTCCGAGC‑GAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTGGGATATCCGTCGcc                                                >  8:68145/1‑148 (MQ=255)
                                                                                                gATCCGCGATTCCGAGC‑GAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTGGGATATCCGTCGCCGAg                                             <  8:135606/148‑1 (MQ=255)
                                                                                                   ccGCGATTCCGAGC‑GAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTGCACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTg                                                              <  8:224858/128‑1 (MQ=255)
                                                                                                   ccGCGATTCCGAGC‑GAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTGCACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTg                                                              >  7:224858/1‑128 (MQ=255)
                                                                                                    cgcgATTCCGAGC‑GAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTgg                                                             >  8:295414/1‑128 (MQ=255)
                                                                                                    cgcgATTCCGAGC‑GAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTgg                                                             <  7:295414/128‑1 (MQ=255)
                                                                                                     gcgATTCCGAGC‑GAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTGGGATATCCGTCACCGAGGACTTc                                       <  2:294059/149‑1 (MQ=255)
                                                                                                            cGAGC‑GAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTGGGATATCCGTCGCCGAGGACTTCCTGCtt                                 <  6:276525/148‑1 (MQ=255)
                                                                                                             gAGC‑GAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTGGGATATcc                                                      <  7:121947/126‑1 (MQ=255)
                                                                                                             gAGC‑GAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTGGGATATcc                                                      >  8:121947/1‑126 (MQ=255)
                                                                                                              aGC‑GAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTGGGATATCCGTCGCCGAGGACTTCCTGCTTGcc                              >  8:116977/1‑149 (MQ=255)
                                                                                                               gC‑GAATGCTGGAGGGTTTGATCGACCGGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTGGGATATCCGTCGCCGAGGACTTCCTGCTTGccc                             <  2:124423/149‑1 (MQ=255)
                                                                                                                   aaTGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTGGGATATCCGTCGCCGAGGACTTCCTGCTTGCCCGGc                          <  5:45989/149‑1 (MQ=255)
                                                                                                                             ggTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTGGGATATCCGTCGCCGAGGACTTCCTGCTTGCCCGGCAGCGTgatgc                >  4:268025/1‑149 (MQ=255)
                                                                                                                               tgagATCGACCGGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCACCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTGGGATATCCGTCGCCGAGGACTTCCTGCTTGCCCGGCAGCGTgatgcga              <  4:295157/146‑1 (MQ=255)
                                                                                                                                 tgATCGACCGGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTGGGATATCCGTCGCCGAGGACTTCCTGCTTGCCCGACAGCGTgatgcgat             <  5:274374/148‑1 (MQ=255)
                                                                                                                                  gATCGACCGGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTGGGATATCCGTCGCCGAGGACTTCCTGCTTGCCCGGCAGCGTgatgcgatg            >  7:24898/1‑148 (MQ=255)
                                                                                                                                    tCGACCGGCCGGCCCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTGGGATATCCGTCGCCGAGGACTTCCTGCTTGCCCGGCAGCGTgatgcgat             >  5:352274/1‑145 (MQ=255)
                                                                                                                                    tCGACCGGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTGGGATATCCGTCGCCGAGGACTTCCTGCTTGCCCGGCAGCGTgatgcgat             <  6:352274/145‑1 (MQ=255)
                                                                                                                                         cgggcggACCTTTGTCCGGCCCTCTACACTATCGACGTGGGCGTGGATTGTAACCTGGTAATCGCACTGCGCAAGGGGATTCGTTACTTTGGGATATTCGTCGCCGAGGATTTCCTGCTTGCGGGGGTGGGTTATGCGAGGCATTCGGc    <  4:244439/149‑1 (MQ=255)
                                                                                                                                          ggccggccATAGGTCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTGGGATATCCGTCGCCGAGGACTTCCTg                                    <  2:303464/116‑1 (MQ=255)
                                                                                                                                          ggccggccATAGGTCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTGGGATATCCGTCGCCGAGGACTTCCTg                                    >  1:303464/1‑116 (MQ=255)
                                                                                                                                           gccggACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTGGGATATCCGTCGCCGAGGACTTCCTGCTTGCCCGGCAGCGTGATGCGATGCAGTCGGc    <  5:53385/147‑1 (MQ=255)
                                                                                                                                           gccggACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTGGGATATCCGTCGCCGAGGACTTCCTGCTTGCCCGGCAGCGTGATGCGATGCAGTCGGc    >  6:53385/1‑147 (MQ=255)
                                                                                                                                           gccggACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTGGGATATCCGTCGCCGAGGACTTCCTGCTTGCCCGGCAGCGTGATGCGATGCAGTCGGCGa  <  7:116977/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                          
CGACACGATGGAGACCGGGATGGGGGCCCGCGAGCAACCCATCATGCTGGTGATCACCACGGCGGGCTCAAGTATTGGCGGCCCCTGCCACCAGTTGATCCGCGATTCCGAGC‑GAATGCTGGAGGGTGTGATCGACCGGCCGGACCTATGGCCGGCCCTCTACACTATCGACCAGGGCGACGATTGGACCAGCGAAATCGCACTGCGCAAGGCGAATCCAAACTTTGGGATATCCGTCGCCGAGGACTTCCTGCTTGCCCGGCAGCGTGATGCGATGCAGTCGGCGA  >  NC_002947/1757586‑1757872

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