Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,416,082 Δ1 bp 5.7% intergenic (+24/+95) PP_3901 → / ← PP_3902 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,416,0820C. 5.7% 222.0 / 6.7 70intergenic (+24/+95)PP_3901/PP_3902hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (32/34);  new base . (2/2);  total (34/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.77e-01

GGGGAATCCCAGGAAGACTGGGGCGCGCTCGTCGAGGCCATGGAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGATAGACACCGGCCAGCAACAAAGGCAGCGAAA  >  NC_002947/4415936‑4416229
                                                                                                                                                  |                                                                                                                                                     
ggggAATCCCAGGAAGACTGGGGCGCGCTCGTCGAGGCCATGGAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCt                                                                                                                                                     >  6:34010/1‑149 (MQ=255)
 gggAATCCCAGGAAGACTGGGGCGCGCTCGTCGAGGCCATGGAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTc                                                                                                                                                    >  6:11161/1‑149 (MQ=255)
        ccAGGAAGACTGGGGCGCGCTCGTCGAGGCCATGGAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCggg                                                                                                                                            >  8:254472/1‑149 (MQ=255)
            gAAGACTGGGGCGCGCTCGTCGAGGCCATGGAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCttt                                                                                                                                       <  5:78572/149‑1 (MQ=255)
            gAAGACTGGGGCGCGCTCGTCGAGGCCATGGAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCt                                                                                                                                         <  4:96470/147‑1 (MQ=255)
            gAAGACTGGGGCGCGCTCGTCGAGGCCATGGAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCt                                                                                                                                         >  3:96470/1‑147 (MQ=255)
               gACTGGGGCGCGCTCGTCGAGGCCATGGAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCTCGAGGCGGGcaa                                                                                                                                         >  7:16934/1‑143 (MQ=255)
               gACTGGGGCGCGCTCGTCGAGGCCATGGAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCTCGAGGCGGGcaa                                                                                                                                         <  8:16934/145‑3 (MQ=255)
                   ggggCGCGCTCGTCGAGGCCATGGAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTgcg                                                                                                                                >  7:50427/1‑149 (MQ=255)
                    gggCGCGCTCGTCGAGGCCATGGAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGAGGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTgc                                                                                                                                 <  2:4726/147‑1 (MQ=255)
                        gcgcTCGTCGAGGGCATGGAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATa                                                                                                                           <  1:12048/149‑1 (MQ=255)
                                 gAGGCCATGGAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGTGGGCTTTTTTTGCGCCATAGAAAAGTGa                                                                                                                  <  5:34010/149‑1 (MQ=255)
                                  aGGCCATGGAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGaa                                                                                                                 <  2:217829/149‑1 (MQ=255)
                                      cATGGAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTAt                                                                                                             >  2:10156/1‑149 (MQ=255)
                                        tGGAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATgg                                                                                                           >  3:48321/1‑149 (MQ=255)
                                            aGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTac                                                                                                       >  3:171162/1‑149 (MQ=255)
                                              gCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTacca                                                                                                     >  2:79661/1‑149 (MQ=255)
                                                aCCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTaccacc                                                                                                   <  4:171162/149‑1 (MQ=255)
                                                    aGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATa                                                                                                                           <  2:107923/121‑1 (MQ=255)
                                                    aGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATa                                                                                                                           >  1:107923/1‑121 (MQ=255)
                                                     gggCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTgc                                                                                                                                 >  1:207254/1‑114 (MQ=255)
                                                     gggCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTgc                                                                                                                                 <  2:207254/114‑1 (MQ=255)
                                                       gCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTATAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACAt                                                                                            <  1:261327/149‑1 (MQ=255)
                                                         gTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACAtt                                                                                           >  8:24378/1‑148 (MQ=255)
                                                            aGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTAt                                                                                        <  7:42044/148‑1 (MQ=255)
                                                             gCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTgcg                                                                                                                                >  6:99639/1‑107 (MQ=255)
                                                             gCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTgcg                                                                                                                                <  5:99639/107‑1 (MQ=255)
                                                             gCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGa                                                                                      <  4:48321/149‑1 (MQ=255)
                                                                 ggCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTgcgc                                                                                                                               >  6:118403/1‑104 (MQ=255)
                                                                 ggCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTgcgc                                                                                                                               <  5:118403/104‑1 (MQ=255)
                                                                         aaGATGATGGCGGCGTTCTGGTGACCTGGGAAAAGCCGCTAGAATACTGAATTCGAAAAACCCACCTTCCCCGCCTCTAGCGGGCTTTTTTTGCGCCATAGAATAGTGAATTATGGGTACCACCATAACATTTATGAGCATACCCAtt                                                                           <  2:259182/148‑1 (MQ=255)
                                                                         aaGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCAtt                                                                           <  1:79661/148‑1 (MQ=255)
                                                                         aaGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCAtt                                                                           <  7:254472/148‑1 (MQ=255)
                                                                          aGATGATGGCGGCGTTCGTGTGCCCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGCAAAGTGAATTATGGGTACCACCAGAACATTTATGAGCAGACCCATTgc                                                                         >  6:310890/1‑148 (MQ=255)
                                                                             tgatgGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAAGGCATTGaca                                                                       >  1:106182/1‑148 (MQ=255)
                                                                                    ggcgTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACCGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCAtt                                                                           <  5:146475/138‑1 (MQ=255)
                                                                                    ggcgTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACCGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCAtt                                                                           >  6:146475/1‑138 (MQ=255)
                                                                                      cgTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGa                                                             <  8:50427/149‑1 (MQ=255)
                                                                                      cgTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGa                                                             <  6:273196/149‑1 (MQ=255)
                                                                                           gtgtGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATg                                                                                                            >  4:212982/1‑97 (MQ=255)
                                                                                           gtgtGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATg                                                                                                            <  3:212982/97‑1 (MQ=255)
                                                                                                  cTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATACCATTTATGAGCAAACCCATTGACACAAAATATGAGcatactcatac                                                 >  8:95399/1‑149 (MQ=255)
                                                                                                  cTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGcatactcatag                                                 >  2:119061/1‑148 (MQ=255)
                                                                                                    gggAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGcatactca                                                    >  3:110662/1‑144 (MQ=255)
                                                                                                    gggAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGcatactca                                                    <  4:110662/144‑1 (MQ=255)
                                                                                                      gAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCGTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGc                                             >  1:215901/1‑149 (MQ=255)
                                                                                                      gAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGc                                             <  5:11161/149‑1 (MQ=255)
                                                                                                               cTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACaa                                                                    >  4:34020/1‑117 (MQ=255)
                                                                                                               cTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACaa                                                                    <  3:34020/117‑1 (MQ=255)
                                                                                                                  gaagaCTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCgg                                  <  2:215901/148‑1 (MQ=255)
                                                                                                                    agaCTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGATa                               <  1:10156/149‑1 (MQ=255)
                                                                                                                     gaCTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGATAg                              <  6:238177/149‑1 (MQ=255)
                                                                                                                      aCTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGATAGa                             >  4:256294/1‑149 (MQ=255)
                                                                                                                           tttCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGATAGACACCg                        <  5:310890/149‑1 (MQ=255)
                                                                                                                            ttCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGATAGACACCgg                       >  7:152663/1‑149 (MQ=255)
                                                                                                                               gACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGTCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATTAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGATAGACACCGGc                      >  5:296298/1‑147 (MQ=255)
                                                                                                                                aCAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAatat                                                               >  3:125589/1‑105 (MQ=255)
                                                                                                                                aCAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAatat                                                               <  4:125589/105‑1 (MQ=255)
                                                                                                                                 cAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGATAGACACCGGCCAg                   >  8:26391/1‑148 (MQ=255)
                                                                                                                                  aaaCCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGATAGACACCGGCCAGc                  >  2:94001/1‑148 (MQ=255)
                                                                                                                                     cccACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGATAGACACCGGCCAGcaac               >  1:81657/1‑148 (MQ=255)
                                                                                                                                       cACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGATAGACACCGGCCAGCAACaaa            <  3:256294/149‑1 (MQ=255)
                                                                                                                                       aaaaaagCCCGCTCGAGGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGATAGACACCGGCCAGcaaca              >  7:10001/7‑147 (MQ=255)
                                                                                                                                          aaagCCCGCTCGAGGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGATAGACACCGGCCAGCAACAAAGGc         <  8:10001/146‑1 (MQ=255)
                                                                                                                                               ccGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAAc                                                                                <  5:60556/73‑1 (MQ=255)
                                                                                                                                               ccGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAAc                                                                                >  6:60556/1‑73 (MQ=255)
                                                                                                                                               ccGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGc                                             >  7:29574/1‑108 (MQ=255)
                                                                                                                                               ccGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGc                                             <  8:29574/108‑1 (MQ=255)
                                                                                                                                                 gCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGATAGACACCGGCCAGCAACAAAGGCAGCGaaa  <  2:106182/149‑1 (MQ=255)
                                                                                                                                                 gCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGATAGACACCGGCCAGCAACAAAGGCAGCGaaa  <  7:26391/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GGGGAATCCCAGGAAGACTGGGGCGCGCTCGTCGAGGCCATGGAAGGCACCGAGGGCGTAAGCCTGGCCTTTGAAGATGATGGCGGCGTTCGTGTGACCTGGGAAAAGCCGCTAGAAGACTGATTTCGACAAACCCACCTTGCCCGCCTCGAGCGGGCTTTTTTTGCGCCATAGAAAAGTGAATTATGGGTACCACCATAACATTTATGAGCAAACCCATTGACACAAAATATGAGCATACTCATACTCGCTCCATCACCGGATAGACACCGGCCAGCAACAAAGGCAGCGAAA  >  NC_002947/4415936‑4416229

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