Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,826,830 G→C 8.0% W81C (TGG→TGC PP_2479 → cytochrome c family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,826,8300GC 8.0% 228.7 / 9.5 75W81C (TGG→TGCPP_2479cytochrome c family protein
Reads supporting (aligned to +/- strand):  ref base G (34/35);  new base C (3/3);  total (37/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.32e-01

TACCTGACGCGCGCTGCCAACTGTGTGGCTTGCCATACCGTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGACATCCTCGCGATCAAGGGTTATCTCGACAGCCTGCAACCG  >  NC_002947/2826684‑2826977
                                                                                                                                                   |                                                                                                                                                     
tACCTGACGCGCGCTGCCAACTGTGTGGCTTGCCATACCGTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGAc                                                                                                                                                     <  2:294951/149‑1 (MQ=255)
   cTGACGCGCGCTGCCAACTGTGTGGCTTGCCATACCGTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGa                                                                                                                                                  >  7:216001/1‑149 (MQ=255)
      aCGCGCGCTGCCAACTGTGTGGCTTGCCATACCGTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTCCATGCACCGATgt                                                                                                                                               <  4:132035/149‑2 (MQ=255)
        gcgcgcTGCCAACTGTGTGGCTTGCCATACCGTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGAc                                                                                                                                              >  2:17830/1‑148 (MQ=255)
         cgcgcTGCCAACTGTGTGGCTTGCCATACCGTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGa                                                                                                                                            <  2:296422/149‑1 (MQ=255)
              tGCCAACTGTGTGGCTTGCCATACCGTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCg                                                                                                                                       >  1:206056/1‑149 (MQ=255)
                ccAACTGTGTGGCTTGCCATACCGTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTc                                                                                                                                        <  7:229900/146‑1 (MQ=255)
                ccAACTGTGTGGCTTGCCATACCGTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTc                                                                                                                                        >  8:229900/1‑146 (MQ=255)
                 cAACTGTGTGGCTTGCCATACCGTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCa                                                                                                                                    <  1:17830/149‑1 (MQ=255)
                    cTGTGTGGCTTGCCATACCGTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTg                                                                                                                                 <  8:32400/149‑1 (MQ=255)
                       tgtgGCTTGCCATACCGTGCCTGGCGGTCAACCCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTgcg                                                                                                                               >  1:46856/1‑149 (MQ=255)
                       tgtgGCTTGCCATACCGTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTgcgc                                                                                                                              >  1:206021/1‑149 (MQ=255)
                         tgGCTTGCCATACCGTGCCTGGCGGTCAACGCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACTTTACCCCTGACCCGCATACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTg                                                                                                                            <  4:234648/149‑1 (MQ=255)
                           gCTTGCCATACCGTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGt                                                                                                                                      <  5:55551/137‑1 (MQ=255)
                           gCTTGCCATACCGTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGt                                                                                                                                      >  6:55551/1‑137 (MQ=255)
                                ccATACCGTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCg                                                                                                                     <  2:229008/149‑1 (MQ=255)
                                 cATACCGTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCgg                                                                                                                    >  2:336834/1‑149 (MQ=255)
                                 cATACCGTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGATGCAGACCgg                                                                                                                    >  8:258872/1‑149 (MQ=255)
                                       gTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGAc                                                                                                                                                     <  7:299238/110‑1 (MQ=255)
                                       gTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCACGGACCCGCAAACCGGCATCGGTGCATGGAc                                                                                                                                                     >  8:299238/1‑110 (MQ=255)
                                                     aCCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGc                                                                                               <  2:46856/149‑1 (MQ=255)
                                                        ccTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGAc                                                                                                       >  5:167696/1‑140 (MQ=255)
                                                        ccTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGAc                                                                                                       <  6:167696/140‑1 (MQ=255)
                                                          ttCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCac                                                                                            >  7:43457/1‑147 (MQ=255)
                                                               cGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTAcc                                                                                     <  3:123317/149‑1 (MQ=255)
                                                                      gtTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACtactac                                                                                      <  7:55788/141‑1 (MQ=255)
                                                                      gtTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACtactac                                                                                      >  8:55788/1‑141 (MQ=255)
                                                                       ttGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGttt                                                                             >  5:90646/1‑149 (MQ=255)
                                                                          aGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTccc                                                                          <  2:206021/149‑1 (MQ=255)
                                                                          aGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTcc                                                                           <  6:90646/148‑1 (MQ=255)
                                                                           gTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTAcc                                                                                     >  3:261608/1‑137 (MQ=255)
                                                                           gTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTAcc                                                                                     <  4:261608/137‑1 (MQ=255)
                                                                           gTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTccc                                                                          >  5:286610/1‑148 (MQ=255)
                                                                            ttCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTg                                                                                                               <  6:32294/112‑1 (MQ=255)
                                                                              cAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCACGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTg                                                                                                               >  5:325637/1‑110 (MQ=255)
                                                                              cAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCACGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTg                                                                                                               <  6:325637/110‑1 (MQ=255)
                                                                               aaGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATAc                                                                     <  8:216001/149‑1 (MQ=255)
                                                                                 gCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTcc                                                                           >  8:306845/1‑141 (MQ=255)
                                                                                 gCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTcc                                                                           <  7:306845/141‑1 (MQ=255)
                                                                                             cTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTcc                                                                           >  4:302404/1‑129 (MQ=255)
                                                                                                 cTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTc                                                   >  1:220891/1‑149 (MQ=255)
                                                                                                           ccAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCAACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTcgc                                                 <  1:248487/142‑1 (MQ=255)
                                                                                                           ccAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTcgc                                                 >  2:248487/1‑141 (MQ=255)
                                                                                                           ccAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATgc                                          >  8:116848/1‑147 (MQ=255)
                                                                                                            cAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGt                                                                                                                                  <  4:238892/61‑1 (MQ=255)
                                                                                                            cAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGt                                                                                                                                  >  3:238892/1‑61 (MQ=255)
                                                                                                            cAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGAc                                         >  3:281099/1‑148 (MQ=255)
                                                                                                                     cccGGACCCGCAAACCGGCATCGGTGGATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGACATCCTCGCGa                               >  1:186110/1‑149 (MQ=255)
                                                                                                                     cccGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGACATCCTCGCGa                               >  4:157988/1‑149 (MQ=255)
                                                                                                                       cGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGACATCCTCGCGATc                             <  4:281099/149‑1 (MQ=255)
                                                                                                                           cccGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCTGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGACATCCTCGCGATCAAgg                         <  7:258872/149‑1 (MQ=255)
                                                                                                                           cccGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGACATCCTCGCGATCAAgg                         <  8:43457/149‑1 (MQ=255)
                                                                                                                           cccGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGACATCCTCGCGATCAAgg                         <  1:336834/149‑1 (MQ=255)
                                                                                                                            ccGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGACATCCTCGCGATCAAgg                         >  5:342638/1‑148 (MQ=255)
                                                                                                                            ccGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGACATCCTCGCGATCAAgg                         <  6:342638/148‑1 (MQ=255)
                                                                                                                              gCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGAc                                         >  7:32129/1‑130 (MQ=255)
                                                                                                                              gCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGAc                                         <  8:32129/130‑1 (MQ=255)
                                                                                                                              gCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGAc                                         >  5:236297/1‑130 (MQ=255)
                                                                                                                              gCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGAc                                         <  6:236297/130‑1 (MQ=255)
                                                                                                                              gCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGACATCCTCGCGATCAAGGGtt                      <  2:220891/149‑1 (MQ=255)
                                                                                                                               cAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGACATCCTCGCGATCAAGGGTTa                     <  6:286610/149‑1 (MQ=255)
                                                                                                                                 aaCCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGACATc                                      >  7:164426/1‑130 (MQ=255)
                                                                                                                                 aaCCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGACATc                                      <  8:164426/130‑1 (MQ=255)
                                                                                                                                   tcgtcATCGGTCCATGCACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGACATCCTCGCGATCAAGGGTTAtctc                 >  6:348769/5‑149 (MQ=255)
                                                                                                                                    cGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAATCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGACATCCTCGCGATCAAGGGTTAtctc                 <  1:109397/148‑1 (MQ=255)
                                                                                                                                    cGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGACATCCTCGCGATCAAGGGTTAtctc                 >  2:109397/1‑148 (MQ=255)
                                                                                                                                     gtcATCGGTCCATGCACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCaa                                                         >  8:10806/3‑107 (MQ=255)
                                                                                                                                     gtcATCGGTCCATGCACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCaa                                                         <  7:10806/105‑1 (MQ=255)
                                                                                                                                      tcATCGGTCCATGCACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTc                                                   <  2:8374/111‑1 (MQ=255)
                                                                                                                                      tcATCGGTCCATGCACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTc                                                   >  1:8374/2‑112 (MQ=255)
                                                                                                                                      gCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGACATCCTCGCGATCAAGGGTTATCTCGAc              <  7:116848/149‑1 (MQ=255)
                                                                                                                                       cATCGGTGCATGGACCCATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGACATCCTCGCGATCAAGGGTTATCTCGACa             <  8:59923/149‑1 (MQ=255)
                                                                                                                                            gTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCATGGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGACATCCTCGCGATCAAGGGTTATCTCGACAGCCt         >  1:100321/1‑149 (MQ=255)
                                                                                                                                            gTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGACATCCTCGCGATCAAGGGTTATCTCGACAGCCTg        >  5:95254/1‑149 (MQ=255)
                                                                                                                                            gTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGACATCCTCGCGATCAAGGGTTATCTCGACAGCCTg        >  2:236387/1‑149 (MQ=255)
                                                                                                                                                  ggACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGACATCCTCGCGATCAAGGGTTATCTCGACAGCCTGCAACCg  >  7:119551/1‑149 (MQ=255)
                                                                                                                                                   gACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGACATCCTCGCGATCAAGGGTTATCTCGACAGCCt         <  7:255085/141‑1 (MQ=255)
                                                                                                                                                   gACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGACATCCTCGCGATCAAGGGTTATCTCGACAGCCt         >  8:255085/1‑141 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TACCTGACGCGCGCTGCCAACTGTGTGGCTTGCCATACCGTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGACATCCTCGCGATCAAGGGTTATCTCGACAGCCTGCAACCG  >  NC_002947/2826684‑2826977

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: