Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,826,868 G→C 5.7% A80P (GCA→CCA)  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,8680GC 5.7% 278.0 / 7.2 88A80P (GCA→CCA) PP_2479cytochrome c family protein
Reads supporting (aligned to +/- strand):  ref base G (40/43);  new base C (3/2);  total (43/45)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.78e-01

GGGCGAATACCTGACGCGCGCTGCCAACTGTGTGGCTTGCCATACCGTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGACATCCTCGCGATCAAGGGTTATCTCGACA  >  NC_002947/2826720‑2827009
                                                                                                                                                    |                                                                                                                                             
ggacGAATACCTGACGCGCGCTGCCAACTGTGTGGCTTGCCATACCGTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGc                                                                                                                                              >  2:1678015/4‑150 (MQ=255)
 ggCGAATACCTGACGCGCGCTGCCAACTGTGTGGCTTGCCATACCGTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCa                                                                                                                                             >  1:1798344/1‑150 (MQ=255)
      aTACCTGACGCGCGCTGCCAACTGTGTGGCTTGCCATACCGTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGAc                                                                                                                                        <  2:1266366/150‑1 (MQ=255)
          cTGACGCGCGCTGCCAACTGTGTGGCTTGCCATACCGTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGAt                                                                                                                                    >  2:2021850/1‑150 (MQ=255)
          cTGACGCGCGCTGCCAACTGTGTGGCTTGCCATACCGTGCCCGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGAt                                                                                                                                    >  2:721385/1‑150 (MQ=255)
           tGACGCGCGCTGCCAACTGTGTGGCTTGCCATACCGTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATg                                                                                                                                   >  1:2137566/1‑150 (MQ=255)
              cgcgcgCTGCCAACTGTGTGGCTTGCCATACCGTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTCCATGCACCGATGCCg                                                                                                                                <  2:233122/150‑1 (MQ=255)
               gcgcgcTGCCAACTGTGTGGCTTGCCATACCGTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGa                                                                                                                               <  2:2137566/150‑1 (MQ=255)
                cgcgcTGCCAACTGTGTGGCTTGCCATACCGTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAg                                                                                                                              <  2:874247/150‑1 (MQ=255)
                     tGCCAACTGTGTGGCTTGCCATACCGTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGt                                                                                                                         <  2:1798344/150‑1 (MQ=255)
                       ccAACTGTGTGGCTTGCCATACCGTGCCTGGCGGTCAACCCTTCGCCGGTGGAGTTGAGTTCCCGTTGCCGCTTGGTTCCCTTTGTTCTCCCAACGTGACCCCGTACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCa                                                                                                                       <  2:1478157/150‑1 (MQ=255)
                       ccAACTGTGTGGCTTGCCATACCGTGCCTGGCGGTCAACCCTTCACCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCa                                                                                                                       <  1:742570/150‑1 (MQ=255)
                       ccAACTGTGTGGCTTGCCATACCGTGCCCGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCa                                                                                                                       <  1:721385/150‑1 (MQ=255)
                         aaCTGTGTGGCTTGCCATACCGTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTCCATGCAcc                                                                                                                                       >  1:1309739/1‑132 (MQ=255)
                         aaCTGTGTGGCTTGCCATACCGTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTCCATGCAcc                                                                                                                                       <  2:1309739/132‑1 (MQ=255)
                          aCTGTGTGGCTTGCCATACCGTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTg                                                                                                                    <  2:950211/150‑1 (MQ=255)
                             gtgtgGCTTGCCATACTTTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTgcgc                                                                                                                 <  1:2043098/150‑1 (MQ=255)
                              tgtgGCTTGCCATACCGTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCt                                                                                                                >  2:212000/1‑150 (MQ=255)
                                tgGCTTGCCATACCGTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGa                                                                                                                               >  1:864028/1‑133 (MQ=255)
                                tgGCTTGCCATACCGTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGa                                                                                                                               <  2:864028/133‑1 (MQ=255)
                                  gCTTGCCATACCGTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAg                                                                                                            >  1:1215560/1‑150 (MQ=255)
                                   cTTGCCCTACCGTGCCTGCACGTCAACCCTGCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTTTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGa                                                                                                           <  1:42089/150‑1 (MQ=255)
                                    ttGCCATAGCGTGCCTGGCGGTCAACCATTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCCGACCGCTGTTTTCTCCCAACATTACCAGGGACCCGCAAACCGGCATAGTTGCATGGACCGATGACGACTTTGTCCGTGCGCTTCCGAc                                                                                                          >  2:232736/1‑150 (MQ=255)
                                    ttGCCATACCGTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGAc                                                                                                          <  2:1215560/150‑1 (MQ=255)
                                      gCCATACCGTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCg                                                                                                        <  1:2271443/150‑1 (MQ=255)
                                        cATACCGTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCAGCATCGGTGCATGGTCCTATTACAAGGTCGTCAGTGAGCTCCacaaaggc                                                                                                      >  2:999961/1‑143 (MQ=255)
                                         aTACCGTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCa                                                                                                     >  2:447450/1‑150 (MQ=255)
                                           aCCGTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCAtt                                                                                                   <  1:999961/150‑1 (MQ=255)
                                            ccGTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTg                                                                                                  <  2:1195034/150‑1 (MQ=255)
                                                 catGGCGGTCAACCCTTAGCCGGTGGTGTAGAGTTCAAGCTGCCGTTCGGCTCGGTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGAACGGCATTGGCCGt                                                                                             <  2:905141/148‑1 (MQ=255)
                                                  cTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGc                                                                                                <  1:1435170/146‑1 (MQ=255)
                                                  cTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGc                                                                                                >  2:1435170/1‑146 (MQ=255)
                                                  cTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTg                                                                                            >  1:681411/1‑150 (MQ=255)
                                                   tGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGa                                                                                           >  1:87902/1‑150 (MQ=255)
                                                     gCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCAt                                                                                                                                            >  1:1449358/1‑99 (MQ=255)
                                                     gCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCAt                                                                                                                                            <  2:1449358/99‑1 (MQ=255)
                                                        gTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATGACCCGGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCa                                                                                      >  1:2123296/1‑150 (MQ=255)
                                                            aCCCTGCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAGCATTACCCCGGATCCGCGGAGCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAgga                                                                                  <  2:536485/150‑3 (MQ=255)
                                                             cccTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCac                                                                                 >  2:2210311/1‑150 (MQ=255)
                                                              ccTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCact                                                                                >  1:367354/1‑150 (MQ=255)
                                                              ccTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCact                                                                                >  1:1677842/1‑150 (MQ=255)
                                                               cTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCacta                                                                               <  2:1737594/150‑1 (MQ=255)
                                                               cTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCacta                                                                               >  2:1238456/1‑150 (MQ=255)
                                                                   gCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACtactac                                                                           <  1:212000/150‑1 (MQ=255)
                                                                            gtTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGttt                                                                  <  1:1678015/150‑1 (MQ=255)
                                                                             ttGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCTCGCCCACGCAAACCGGCATCGGGGCATGGACCGAGGACTAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTc                                                                 <  2:1476959/150‑1 (MQ=255)
                                                                             ttGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTc                                                                 >  2:2160000/1‑150 (MQ=255)
                                                                              tGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTcc                                                                <  1:2019069/150‑1 (MQ=255)
                                                                               gAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCaa                                                                                     <  2:60116/128‑1 (MQ=255)
                                                                               gAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCaa                                                                                     >  1:60116/1‑128 (MQ=255)
                                                                                        cTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATAcctcc                                                      <  2:87902/150‑1 (MQ=255)
                                                                                            cGTTCGGCTCGCTGTTTTCTCCCAACCTTACTCGGCACCCGCAAACCGGCATCGGTGCATGGACCGATGACGCGTTCGTCAGTTCGCTTCTGACCGGCATTGGCCGTGACGGTAAGCACTACTACCCGGCGTTTCCCTATACCTCCcaca                                                  <  2:143293/150‑1 (MQ=255)
                                                                                                  gCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCaaaa                                            >  1:758361/1‑150 (MQ=255)
                                                                                                  gCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCaaaa                                            >  1:2027198/1‑150 (MQ=255)
                                                                                                   cTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAATCCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCa                                                                                  >  1:1948865/1‑111 (MQ=255)
                                                                                                   cTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCa                                                                                  <  2:1948865/111‑1 (MQ=255)
                                                                                                      gCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTc                                        <  1:1238456/150‑1 (MQ=255)
                                                                                                             tCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGAt                                 <  1:588987/150‑1 (MQ=255)
                                                                                                                cccAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCa                                                                                      >  2:603809/1‑94 (MQ=255)
                                                                                                                cccAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCa                                                                                      <  1:603809/94‑1 (MQ=255)
                                                                                                                cccAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAAt                                           <  2:775605/137‑1 (MQ=255)
                                                                                                                cccAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAAt                                           >  1:775605/1‑137 (MQ=255)
                                                                                                                 ccAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGACa                             <  1:2160000/150‑1 (MQ=255)
                                                                                                                  cAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAg                                                                                    <  2:259458/94‑1 (MQ=255)
                                                                                                                  cAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAg                                                                                    >  1:259458/1‑94 (MQ=255)
                                                                                                                  cAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCactac                                                                              <  2:1242953/100‑1 (MQ=255)
                                                                                                                  cAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCactac                                                                              >  1:1242953/1‑100 (MQ=255)
                                                                                                                   aaCATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCa                                                                                  <  1:868224/95‑1 (MQ=255)
                                                                                                                   aaCATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCa                                                                                  >  2:868224/1‑95 (MQ=255)
                                                                                                                       ttACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGACATCCTcg                       >  1:1340980/1‑150 (MQ=255)
                                                                                                                         aCCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGACATCCTcgcg                     >  2:2002518/1‑150 (MQ=255)
                                                                                                                          ccccGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGACATCCTCGCGa                    <  1:2002518/150‑1 (MQ=255)
                                                                                                                           cccGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGACATCCTCGCGAt                   >  1:2112880/1‑150 (MQ=255)
                                                                                                                           cccGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGACATCCTCGCGAt                   >  2:2380005/1‑150 (MQ=255)
                                                                                                                            ccGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCaa                                                                                     <  1:1114007/83‑1 (MQ=255)
                                                                                                                            ccGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCaa                                                                                     >  2:1114007/1‑83 (MQ=255)
                                                                                                                            ccGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGACATCCTCGCGATc                  >  1:1840139/1‑150 (MQ=255)
                                                                                                                              ggACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTgg                                                                                                 <  2:1251730/69‑1 (MQ=255)
                                                                                                                              ggACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTgg                                                                                                 >  1:1251730/1‑69 (MQ=255)
                                                                                                                               gACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGACATCCTCGCGATCAAg               <  1:1089434/150‑1 (MQ=255)
                                                                                                                                  ccGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGACATCCTCGCGATCAAGGGt            >  2:1135197/1‑150 (MQ=255)
                                                                                                                                       aaCCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGACATCCTCGCGATCAAGGGTTAtct       <  1:1674521/150‑1 (MQ=255)
                                                                                                                                        aCCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGc                                                                                   >  2:557417/1‑73 (MQ=255)
                                                                                                                                        aCCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGc                                                                                   <  1:557417/73‑1 (MQ=255)
                                                                                                                                         tcgtcATCGGTCCATGCACCGATGCCGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGACATCCTCGCGATCAAGGGTTATCTCg     >  2:1396074/5‑150 (MQ=255)
                                                                                                                                            tcATCGGTCCATGCACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGACATCCTCGCGATCAAGGGTTATCTCGACa  >  2:1855881/2‑150 (MQ=255)
                                                                                                                                              aTCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGACATCCTCGCGATCAAGGGTTAtc        <  2:437589/142‑1 (MQ=255)
                                                                                                                                              aTCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGACATCCTCGCGATCAAGGGTTAtc        >  1:437589/1‑142 (MQ=255)
                                                                                                                                                    gCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTcc                                                                <  2:1467019/80‑1 (MQ=255)
                                                                                                                                                    gCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTcc                                                                >  1:1467019/1‑80 (MQ=255)
                                                                                                                                                    |                                                                                                                                             
GGGCGAATACCTGACGCGCGCTGCCAACTGTGTGGCTTGCCATACCGTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGACATCCTCGCGATCAAGGGTTATCTCGACA  >  NC_002947/2826720‑2827009

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