Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 4,031,607 T→C 5.5% T25A (ACG→GCG)  U712_20820 ← Protein jag

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068814,031,6070TC 5.5% 221.9 / 3.4 73T25A (ACG→GCG) U712_20820Protein jag
Reads supporting (aligned to +/- strand):  ref base T (28/40);  new base C (2/2);  total (30/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

TCGTTTGAAGATATAACTTCGCTTGCTGAATCGGATCAATCTTTTCAACCAGCTTTACGATCGCAGGCTTTTTACCGAAAATACCGAGAAACCCTTTGTTTCCCTCTTCAATAACAGTAATTTCTACCTTATCCTTCGTTAATCCCAACTCTTGCAATCCGGACTGTACTGCTTCATCGACATTTCGCCCTGCAGCAGTCACATTCCTCACTTTTTCTTTCCTCCGGCTTTTTGCGGCTCAGGATTTTTTTTAATATCAGGTCCTTTAATGAGG  >  CP006881/4031469‑4031742
                                                                                                                                          |                                                                                                                                       
tCGTTTGAAGATATAACTTCGCTTGCTGAATCGGATCAATCTTTTCAACCAGCTTTACGATCGCAGGCTTTTTACCGAAAATACCGAGAAACCCTTTGTTTCCCTCTTCAATAACAGTAATTTCTACCTTATCCTTCGt                                                                                                                                         <  7:201267/139‑1 (MQ=255)
          atatAACTTCGCTTGCTGAATCGGATCAATCTTTTCAACCAGCTTTACGATCGCAGGCTTTTTACCGAAAATACCGAGAAACCCTTTGTTTCCCTCTTCAATAACAGTAATTTCTACCTTATCCTTCGTTAATCCCAAc                                                                                                                               <  3:225392/139‑1 (MQ=255)
               aCTTCGCTTGCTGAATCGGATCAATCTTTTCAACCAGCTTTACGATCGCAGGCTTTTTACCGAAAATACCGAGAAACCCTTTGTTTCCCTCTTCAATAACAGTAATTTCTACCTTATCCTTCGTTAATCCCAACTCtt                                                                                                                           <  4:60638/138‑1 (MQ=255)
               aCTTCGCTTGCTGAATCGGATCAATCTTTTCAACCAGCTTTACGATCGCAGGCTTTTTACCGAAAATACCGAGAAACCCTTTGTTTCCCTCTTCAATAACAGTAATTTCTACCTTATCCTTCGTTAATCCCAACTCtt                                                                                                                           <  3:272588/138‑1 (MQ=255)
               aCTTCGCTTGCTGAATCGGATCAATCTTTTCAACCAGCTTTACGATCGCAGGCTTTTTACCGAAAATACCGAGAAACCCTTTGTTTCCCTCTTCAATAACAGTAATTTCTACCTTATCCTTCGTTAATCCCAACTCtt                                                                                                                           <  2:291437/138‑1 (MQ=255)
               aCTTCGCTTGCTGAATCGGATCAATCTTTTCAACCAGCTTTACGATCGCAGGCTTTTTACCGAAAATACCGAGAAACCCTTTGTTTCCCTCTTCAATAACAGTAATTTCTACCTTATCCTTCGTTAATCCCAACTCTTg                                                                                                                          <  6:353979/139‑1 (MQ=255)
                   cGCTTGCTGAATCGGATCAATCTTTTCAACCAGCTTTACGATCGCAGGCTTTTTACCGAAAATTCCGAGAAACCCTTTGTTTCCCTCTTCAATAACAGTAATTTCTACCTTATCCTTCGTTAATCCCATCTCTTGCAAt                                                                                                                      <  4:44138/139‑1 (MQ=255)
                   cGCTTGCTGAATCGGATCAATCTTTTCAACCAGCTTTACGATCGCAGGCTTTTTACCGAAAATACCGAGAAACCCTTTGTTTCCCTCTTCAATAACAGTAATTTCTACCTTATCCTTCGTTAATCCCAACTCTTGCAAt                                                                                                                      >  8:103646/1‑139 (MQ=255)
                   cGCTTGCTGAATCGGATCAATCTTTTCAACCAGCTTTACGATCGCAGGCTTTTTACCGAAAATACCGAGAAACCCTTTGTTTCCCTCTTCAATAACAGTAATTTCTACCTTATCCTTCGCTAATCCCAACTCTTGCaa                                                                                                                       >  1:215247/1‑138 (MQ=255)
                   cGCTTGCTGAATCGGATCAATCTTTTCAACCAGCTTTACGATCGCAGGCTTTTTACCGAAAATACCGAGAAACCCTTTGTTTCCCTCTTCAATAACAGTAATTTCTACCTTATCCTTCGCTAATCCCAACTCTTGCAAt                                                                                                                      >  1:210024/1‑139 (MQ=255)
                    gCTTGCTGAATCGGATCAATCTTTTCAACCAGCTTTACGATCGCAGGCTTTTTACCGAAAATACCGAGAAACCCTTTGTTTCCCTCTTCAATAACAGTAATTTCTACCTTATCCTTCGTTAATCCCAACTCTTGCAAt                                                                                                                      <  5:345912/138‑1 (MQ=255)
                    gCTTGCTGAATCGGATCAATCTTTTCAACCAGCTTTACGATCGCAGGCTTTTTACCGAAAATACCGAGAAACCCTTTGTTTCCCTCTTCAATAACAGTAATTTCTACCTTATCCTTCGTTAATCCCAACTCTTGCAATc                                                                                                                     >  5:332096/1‑139 (MQ=255)
                     cTTGCTGAATCGGATCAATCTTTTCAACCAGCTTTACGATCGCAGGCTTTTTACCGAAAATACCGAGAAACCCTTTGTTTCCCTCTTCAATAACAGTAATTTCTACCTTATCCTTCGTTAATCCCAACTCTTGCAATcc                                                                                                                    >  1:150172/1‑139 (MQ=255)
                     cTTGCTGAATCGGATCAATCTTTTCAACCAGCTTTACGATCGCAGGCTTTTTACCGAAAATACCGAGAAACCCTTTGTTTCCCTCTTCAATAACAGTAATTTCTACCTTATCCTTCGTTAATCCCAACTCTTGCAATcc                                                                                                                    >  5:349829/1‑139 (MQ=255)
                        gatGAATCGGATCAATCTTTTCAACCAGCTTTACGATCGCAGGCTTTTTACCGAAAATACCGAGAAACCCTTTGTTTCCCTCTTCAATAACAGTAATTTCTACCTTATCCTTCGTTAATCCCAACTCTTGCAATCCGGa                                                                                                                 >  2:201711/3‑139 (MQ=255)
                        gCTGAATCGGATCAATCTTTTCAACCAGCTTTACGATCGCAGGCTTTTTACCGAAAATACCGAGAAACCCTTTGTTTCCCTCTTCAATAACAGTAATTTCTACCTTATCCTTCGTTAATCCCAACTCTTGCAATCCg                                                                                                                   <  8:351123/137‑1 (MQ=255)
                         cTGAATCGGATCAATCTTTTCAACCAGCTTTACGATCGCAGGCTTTTTACCGAAAATACCGAGAAACCCTTTGTTTCCCTCTTCAATAACAGTAATTTCTACCTTATCCTTCGTTAATCCCAACTCTTGCAATCCGGa                                                                                                                 <  8:55600/138‑1 (MQ=255)
                         cTGAATCGGATCAATCTTTTCAACCAGCTTTACGATCGCAGGCTTTTTACCGAAAATACCGAGAAACCCTTTGTTTCCCTCTTCAATAACAGTAATTTCTACCTTATCCTTCGTTAATCCCAACTCTTGCAATCCGGa                                                                                                                 <  1:249390/138‑1 (MQ=255)
                         cTGAATCGGATCAATCTTTTCAACCAGCTTTACGATCGCAGGCTTTTTACCGAAAATACCGAGAAACCCTTTGTTTCCCTCTTCAATAACAGTAATTTCTACCTTATCCTTCGTTAATCCCAACTCTTGCAATCCGGa                                                                                                                 <  5:28048/138‑1 (MQ=255)
                         cTGAATCGGATCAATCTTTTCAACCAGCTTTACGATCGCAGGCTTTTTACCGAAAATACCGAGAAACCCTTTGTTTCCCTCTTCAATAACAGTAATTTCTACCTTATCCTTCGTTAATCCCAACTCTTGCAATCCGGa                                                                                                                 <  4:210920/138‑1 (MQ=255)
                         cTGAATCGGATCAATCTTTTCAACCAGCTTTACGATCGCAGGCTTTTTACCGAAAATACCGAGAAACCCTTTGTTTCCCTCTTCAATAACAGTAATTTCTACCTTATCCTTCGTTAATCCCAACTCTTGCAATCCGGa                                                                                                                 <  3:73015/138‑1 (MQ=255)
                         cTGAATCGGATCAATCTTTTCAACCAGCTTTACGATCGCAGGCTTTTTACCGAAAATACCGAGAAACCCTTTGTTTCCCTCTTCAATAACAGTAATTTCTACCTTATCCTTCGTTAATCCCAACTCTTGCAATCCGGa                                                                                                                 <  6:332096/138‑1 (MQ=255)
                         cTGAATCGGATCAATCTTTTCAACCAGCTTTACGATCGCAGGCTTTTTACCGAAAATACCGAGAAACCCTTTGTTTCCCTCTTCAATAACAGTAATTTCTACCTTATCCTTCGTTAATCCCAACTCTTGCAATCCGGa                                                                                                                 <  3:275599/138‑1 (MQ=255)
                         cTGAATCGGATCAATCTTTTCAACCAGCTTTACGATCGCAGGCTTTTTACCGAAAATACCGAGAAACCCTTTGTTTCCCTCTTCAATAACAGTAATTTCTACCTTATCCTTCGTTAATCCCAACTCTTGCAATCCGGa                                                                                                                 <  8:49527/138‑1 (MQ=255)
                         cTGAATCGGATCAATCTTTTCAACCAGCTTTACGATCGCAGGCTTTTTACCGAAAATACCGAGAAACCCTTTGTTTCCCTCTTCAATAACAGTAATTTCTACCTTATCCTTCGTTAATCCCAACTCTTGCAATCCGGa                                                                                                                 <  6:372867/138‑1 (MQ=255)
                         cTGAATCGGATCAATCTTTTCAACCAGCTTTACGATCGCAGGCTTTTTACCGAAAATACCGAGAAACCCTTTGTTTCCCTCTTCAATAACAGTAATTTCTACCTTATCCTTCGTTAATCCCAACTCTTGCAATCCGGa                                                                                                                 <  7:203417/138‑1 (MQ=255)
                         cTGAATCGGATCAATCTTTTCAACCAGCTTTACGATCGCAGGCTTTTTACCGAAAATACCGAGAAACCCTTTGTTTCCCTCTTCAATAACAGTAATTTCTACCTTATCCTTCGTTAATCCCAACTCTTGCAATCCGGa                                                                                                                 <  7:297603/138‑1 (MQ=255)
                         cTGAATCGGATCAATCTTTTCAACCAGCTTTACGATCGCAGGCTTTTTACCGAAAATACCGAGAAACCCTTTGTTTCCCTCTTCAATAACAGTAATTTCTACCTTATCCTTCGTTAATCCCAACTCTTGCAATCCGGa                                                                                                                 <  2:150172/138‑1 (MQ=255)
                         cTGAATCGGATCAATCTTTTCAACCAGCTTTACGATCGCAGGCTTTTTACCGAAAATACCGAGAAACCCTTTGTTTCCCTCTTCAATAACAGTAATTTCTACCTTATCCTTCGGTAATCCCAACTCTTGCAATCCGGa                                                                                                                 <  6:236011/138‑1 (MQ=255)
                             aTCGGATCAATCTTTTCAACCAGCTTTACGATCGCAGGCTTTTTACCGAAAATACCGAGACACCCTTTGTTTCCCTCTTCAATAACAGTAATTTCTACCTTATCCTTCGTTAATCCCAACTCTTGCAATCCGGACTGTa                                                                                                            <  4:272994/139‑1 (MQ=255)
                                      aTCTTTTCAACCAGCTTTACGATCGCAGGCTTTTTACCGAAAATACCGAGAAACCCTTTGTTTCCCTCTTCAATAACAGTAATTTCTACCTTATCCTTCGTTAATCCCAACTCTTGCAATCCGGACTGTACTGCTTCAt                                                                                                   <  6:349829/139‑1 (MQ=255)
                                           ttCAACCAGCTTTACGATCGCAGGCTTTTTACCGAAAATACCGAGAAACCCTTTGTTTCCCTCTTCAATAACAGTAATTTCTACCTTATCCTTCGTTAATCCCAACTCTTGCAATCCGGACTGTACTGCTTCATc                                                                                                  <  7:135757/135‑1 (MQ=255)
                                            tCAACCAGCTTTACGATCGCAGGCTTTTTACCGAAAATACCGAGAAACCCTTTGTTTCCCTCTTCAATAACAGTAATTTCTACCTTATCCTTCGTTAATCCCAACTCTTGCAATCCGGACTGTACTGCTTCATCGACAt                                                                                             >  5:231333/1‑139 (MQ=255)
                                                  aGCTTACCGATCGCAGGCTTTTACCCAAAAAACCCGAGAAACCCTTTGTTACCATCTTCAAAAACAGTATTTTATACCTTATCCTCCTTAAAACCCACCTCTTCCAACCCGGAGTGTCCTGCTTCATCGACATTTCGcc                                                                                       <  6:75591/139‑1 (MQ=255)
                                                   gCTTTACGATCGCAGGCTTTTTACCGAAAATACCGAGAAACCCTTTGTTTCCCTCTTCAATAACAGTAATTTCTACCTTATCCTTCGTTAATCCCAACTCTTGCAATCCg                                                                                                                   >  8:3678/1‑110 (MQ=255)
                                                   gCTTTACGATCGCAGGCTTTTTACCGAAAATACCGAGAAACCCTTTGTTTCCCTCTTCAATAACAGTAATTTCTACCTTATCCTTCGTTAATCCCAACTCTTGCAATCCg                                                                                                                   <  7:3678/110‑1 (MQ=255)
                                                   gCTTTACGATCGCAGGCTTTTTACCGAAAATACCGAGAAACCCTTTGTTTCCCTCTTCAATAACAGTAATTTCTACCTTATCCTTCGTTAATCCCAACTCTTGCAATCCGGACTGTACTGCTTCATCGACATTTCGccc                                                                                      <  6:231333/139‑1 (MQ=255)
                                                         cGATCGCAGGCTTTTTACCGAAAATACCGAGAAACCCTTTGTTTCCCTCTTCAATAACAGTAATTTCTACCTTATCCTTCGTTAATCCCAACTCTTGCAATCCGGACTGTACTGCTTCATCGACATTTCGCCCTgcag                                                                                 <  5:185763/138‑1 (MQ=255)
                                                         cGATCGCAGGCGTTTTACCGACAATACCGAGAAACCCTTTGTTTCCCTCTTCAATAACAGTAATTTCTACCTTATCCTTCGTTAATCCCAACTCTTGCAATCCGGACTGTACTGCTTCATCGACATTTCGCCCTgcag                                                                                 <  3:139247/138‑1 (MQ=255)
                                                              gCAGGCTTTTTACCGAAAATACCGAGAAACCCTTTGTTTCCCTCTTCAATAACAGTAATTTCTACCTTATCCTTCGTTAATCCCAACTCTTGCAATCCGGACTGTACTGCTTCATCGCCATTTCGCCCGGCAGCAGTca                                                                           >  7:325736/1‑139 (MQ=255)
                                                              gCAGGCTTTATACCGAAAATACCGAGAAACCCTTTGTTTCCCTCTTCAATAACAGTAATTTCTACCTTATACTTCGTTAATCCCAACTCTAGCAATCAGTACTGAACTGATTCATCGCAATTTCTCCATGCAGAAGTca                                                                           >  4:37248/1‑139 (MQ=255)
                                                                aGGCTTTTTACCGAAAATACCGAGAAACCCTTTGTTTCCCTCTTCAATAACAGTAATTTCTACCTTATCCTTCGTTAATCCCAACTCTTGCAATCCGGACTGTACTGCTTCATCGACATTTCGCCCTGCAGCAGTcaca                                                                         <  2:243630/139‑1 (MQ=255)
                                                                        tACCGAAAATACCGAGAAACCCTTTGTTTCCCTCTTCAATAACAGTAATTTCTACCTTATCCTTCGTTAATCCCAACTCTTGCAATCCGGa                                                                                                                 <  6:248879/91‑1 (MQ=255)
                                                                        tACCGAAAATACCGAGAAACCCTTTGTTTCCCTCTTCAATAACAGTAATTTCTACCTTATCCTTCGTTAATCCCAACTCTTGCAATCCGGa                                                                                                                 >  5:248879/1‑91 (MQ=255)
                                                                        tACCGAAAATACCGAGAAACCCTTTGTTTCCCTCTTCAATAACAGTAATTTCTACCTTATCCTTCGTTAATCCCAACTCTTGCAATCCGGACTGTACTGCTTCATCGACATTTCGCCCTGCAGCAGTCACATTCCTCa                                                                  <  4:315424/138‑1 (MQ=255)
                                                                                                      ccTCTTCAATAACAGTAATTTCTACCTTATCCTTCGTTAATCCCAACTCTTGCAATCCGGACTGTACTGCTTCATCGACATTTCGCCCTgcag                                                                                 >  6:273196/1‑93 (MQ=255)
                                                                                                      ccTCTTCAATAACAGTAATTTCTACCTTATCCTTCGTTAATCCCAACTCTTGCAATCCGGACTGTACTGCTTCATCGACATTTCGCCCTgcag                                                                                 <  5:273196/93‑1 (MQ=255)
                                                                                                      ccTCTTCAATAACAGTAATTTCTACCTTATCCTTCGTTAATCCCAACTCTTGCAATCCGGACTGTACTGCTTCATCGACATTTCGCCCTGCAGCAGTCACATTCCTCACTTTTTCTTTCCTCCGGCTTTTTGCGGCt                                     >  1:324480/1‑137 (MQ=255)
                                                                                                      ccTCTTCAATAACAGTAATTTCTACCTTATCCTTCGTTAATCCCAACTCTTGCAATCCGGACTGTACTGCTTCATCGACATTTCGCCCTGCAGCAGTCACATTCCTCACTTTTTCTTTCCTCCGGCTTTTTGCGGCTCa                                   <  2:324480/139‑1 (MQ=255)
                                                                                                        tctTCAATAACAGTAATTTCTACCTTATCCTTCGTTAATCCCAACTCTTGCAATCCGGACTGTACTGCTTCATCGACATTTCGCCCTGCAGCAGTCACATTCCTCACTTTTTCTTTCCTCCGGCTTTTTGCGGCTCAgg                                 >  4:76736/1‑139 (MQ=255)
                                                                                                          ttCAATAACAGTAATTTCTACCTTATCCTTCGTTAATCCCAACTCTTGCAATCCGGACTGTACTGCTTCATCGACATTTCGCCCTGCAGCAGTCACATTCCTCACTTTTTCTTTCCTCCGGCTTTTTGCGGCTCAGGAt                               >  7:157046/1‑139 (MQ=255)
                                                                                                             aaTAACAGTAATTTCTACCTTATCCTTCGCTAATCCCAACTCTTGCAATCCGGACTGTACTGCTTCATCGACATTTCGCCCTGCAGCAGTCACATTCCTCACTTTTTCTTTCCTCCGGCTTTTTGCGGCTCAGGAttt                             <  2:210024/138‑1 (MQ=255)
                                                                                                             aaTAACAGTAATTTCTACCTTATCCTTCGCTAATCCCAACTCTTGCAATCCGGACTGTACTGCTTCATCGACATTTCGCCCTGCAGCAGTCACATTCCTCACTTTTTCTTTCCTCCGGCTTTTTGCGGCTCAGGAttt                             <  2:215247/138‑1 (MQ=255)
                                                                                                                    gttaTTTCTACCTTATCCTTCGTTAATCCCAAATCTTGCAATCCGGACTTTACTGCTTCATCGACATTTCGCCCTGCACCAGTCACATTCCTCACTTTTTCTTTCCTCAGGCTTTTTGCGGCTCAGCATTTTTTTTat                      >  8:182423/4‑137 (MQ=255)
                                                                                                                    gTAATTTCTACCTTATCCTTCGTTAATCCCAACTCTTGCAATCCGGACTGTACTGCTTCATCGACATTTCGCCCTGCAGCAGTCACATTCCTCACTTTTTCTTTCCTCCGGCTTTTTGCGGCTCAGGATTTTTTTTaa                      >  6:3231/1‑138 (MQ=255)
                                                                                                                    gTAATTTCTACCTTATCCTTCGTTAATCCCAACTCTTGCAATCCGGACTGTACTGCTTCATCGACATTTCGCCCAGCAGAAGTCACATTCCTCACTTTTTCTTTCCTCAGGCTTTTTGAGGCTCAGGAATTTTTTTaa                      >  4:46700/1‑138 (MQ=255)
                                                                                                                       aTTTCTACCTTATCCTTCGTTAATCCCAACTCTTGCAATCCGGACTGTACTGCTTCATCGACATTTCGCCCTGCAGCAGTCACATTCCTCACTTTTTCTTTCCTCCGGCTTTTTGCGGCTCAGGATTTTTTTTAATATc                  <  1:201711/139‑1 (MQ=255)
                                                                                                                          tCTACCTTATCCTTCGTTAATCCCAACTCTTGCAATCCGGACTGTACTGCTTCATCGACATTTCGCCCTGCAGCAGTCACATTCCTCACTTTTTCTTTCCTCCGGCTTTTTGCGGCTCAGGATTTTTTTTAATATCAg                <  5:340113/138‑1 (MQ=255)
                                                                                                                                ttATCCTTCGTTAATCCCAACTCTTGCAATCCGGACTGTACTGCTTCATCGACATTTCGCCCTGCAGCAGTCACATTCCTCACTTTTTCTTTCCTCCGGCTTTTTGCGGCTCAGGATTTTTTTTAATATCAGGTCCtt          >  4:308616/1‑138 (MQ=255)
                                                                                                                                 tATCCTTCGTTAATCCCAACTCTTGCAATCCGGACTGTACTGCTTCATCGACATTTCGCCCTGCAGCAGTCACATTCCTCACTTTTTCTTTCCTCCGGCTTTTTGCGGCTCAGGATTTTTTTTAATATCAGGTCCttt         >  3:150825/1‑138 (MQ=255)
                                                                                                                                 tATCCTTCGTTAATCCCAACTCTTGCAATCCGGACTGTACTGCTTCATCGACATTTCGCCCTGCAGCAGTCACATTCCTCACTTTTTCTTTCCTCCGGCTTTTTGCGGCTCAGGATTTTTTTTAATATCAGGTCCttt         >  7:145313/1‑138 (MQ=255)
                                                                                                                                  aTCCTTCGTTAATCCCAACTCTTGCAATCCGGACTGTACTGCTTCATCGACATTTCGCCCTGCAGCAGTCACATTCCTCa                                                                  >  4:99221/1‑80 (MQ=255)
                                                                                                                                  aTCCTTCGTTAATCCCAACTCTTGCAATCCGGACTGTACTGCTTCATCGACATTTCGCCCTGCAGCAGTCACATTCCTCa                                                                  <  3:99221/80‑1 (MQ=255)
                                                                                                                                   tCCTTCGTTAATCCCAACTCTTGCAATCCGGACTGTACTGCTTCATCGACATTTCGCCCTGcagcag                                                                              <  6:282330/67‑1 (MQ=255)
                                                                                                                                   tCCTTCGTTAATCCCAACTCTTGCAATCCGGACTGTACTGCTTCATCGACATTTCGCCCTGcagcag                                                                              >  5:282330/1‑67 (MQ=255)
                                                                                                                                   tCCTTCGTTAATCCCAACTCTTGCAATCCGGACTGTACTGCTTCATCGACATTTCGCCCTGCAGCAGTCACATTCCTCACTTTTTCTTTCCTCCGGCTTTTTGCGGCTCAGGATTTTTTTTAATATc                  <  5:6520/127‑1 (MQ=255)
                                                                                                                                   tCCTTCGTTAATCCCAACTCTTGCAATCCGGACTGTACTGCTTCATCGACATTTCGCCCTGCAGCAGTCACATTCCTCACTTTTTCTTTCCTCCGGCTTTTTGCGGCTCAGGATTTTTTTTAATATCAGGTCGTTTAAt      >  1:169948/1‑139 (MQ=255)
                                                                                                                                   tCCTTCGTTAATCCCAACTCTTGCAATCCGGACTGTACTGCTTCATCGACATTTCGCCCTGCAGCAGTCACATTCCTCACTTTTTCTTTCCTCCGGCTTTTTGCGGCTCAGGATTTTTTTTAATATCAGGTCCTTTAAt      >  4:230288/1‑139 (MQ=255)
                                                                                                                                   tCCTTCGTTAATCCCAACTCTTGCAATCCGGACTGTACTGCTTCATCGACATTTCGCCCTGCAGCAGTCACATTCCTCACTTTTTCTTTCCTCCGGCGTTTTGCGGCTCAGGATTTTTTTTAATATc                  >  6:6520/1‑127 (MQ=255)
                                                                                                                                        cGTTAATCCCAACTCTTGCAATCCGGACTGTACTGCTTCATCGACATTTCGCCCTGCAGCAGTCACATTc                                                                      <  6:198834/70‑1 (MQ=255)
                                                                                                                                        cGTTAATCCCAACTCTTGCAATCCGGACTGTACTGCTTCATCGACATTTCGCCCTGCAGCAGTCACATTc                                                                      >  5:198834/1‑70 (MQ=255)
                                                                                                                                        cGTTAATCCCAACTCTTGCAATCCGGACTGTACTGCTTCATCGACATTTCGCCCTGCAGCAGTCACATTCCTCACTTTTTCTTTCCTCCGGCTTTTTGCGGCTCAGGAttttttt                         >  2:176626/1‑115 (MQ=255)
                                                                                                                                        cGTTAATCCCAACTCTTGCAATCCGGACTGTACTGCTTCATCGACATTTCGCCCTGCAGCAGTCACATTCCTCACTTTTTCTTTCCTCCGGCTTTTTGCGGCTCAGGAttttttt                         <  1:176626/115‑1 (MQ=255)
                                                                                                                                        cGTTAATCCCAACTCTTGCAATCCGGACTGTACTGCTTCATCGACATTTCGCCCTGCAGCAGTCACATTCCTCACTTTTTCTTTCCTCCGGCTTTTTGCGGCTCAGGATTTTTTTTAATATCAGGTCCTTTAATGAgg  >  3:52753/1‑138 (MQ=255)
                                                                                                                                          |                                                                                                                                       
TCGTTTGAAGATATAACTTCGCTTGCTGAATCGGATCAATCTTTTCAACCAGCTTTACGATCGCAGGCTTTTTACCGAAAATACCGAGAAACCCTTTGTTTCCCTCTTCAATAACAGTAATTTCTACCTTATCCTTCGTTAATCCCAACTCTTGCAATCCGGACTGTACTGCTTCATCGACATTTCGCCCTGCAGCAGTCACATTCCTCACTTTTTCTTTCCTCCGGCTTTTTGCGGCTCAGGATTTTTTTTAATATCAGGTCCTTTAATGAGG  >  CP006881/4031469‑4031742

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: ‑