Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,826,873 G→C 5.7% 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,8730GC 5.7% 277.9 / 5.9 87W81C (TGG→TGCPP_2479cytochrome c family protein
Reads supporting (aligned to +/- strand):  ref base G (38/44);  new base C (3/2);  total (41/46)
Fisher's exact test for biased strand distribution p-value = 6.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.17e-01

ATACCTGACGCGCGCTGCCAACTGTGTGGCTTGCCATACCGTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGACATCCTCGCGATCAAGGGTTATCTCGACAGCCTGCAACCGGT  >  NC_002947/2826726‑2827022
                                                                                                                                                   |                                                                                                                                                     
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                                                                                                                                            <  2:864028/133‑1 (MQ=255)
                          tgGCTTGCCATACCGTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGa                                                                                                                                            >  1:864028/1‑133 (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                                                                                                             >  2:1435170/1‑146 (MQ=255)
                                            cTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGc                                                                                                             <  1:1435170/146‑1 (MQ=255)
                                            cTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTg                                                                                                         >  1:681411/1‑150 (MQ=255)
                                             tGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGa                                                                                                        >  1:87902/1‑150 (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:2027198/1‑150 (MQ=255)
                                                                                            gCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCaaaa                                                         >  1:758361/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                                                                                                   <  1:603809/94‑1 (MQ=255)
                                                                                                          cccAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCa                                                                                                   >  2:603809/1‑94 (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                                                                                           >  1:1242953/1‑100 (MQ=255)
                                                                                                            cAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCactac                                                                                           <  2:1242953/100‑1 (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                                                                                                              >  1:1251730/1‑69 (MQ=255)
                                                                                                                        ggACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTgg                                                                                                              <  2:1251730/69‑1 (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                                                                             >  1:1467019/1‑80 (MQ=255)
                                                                                                                                              gCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTcc                                                                             <  2:1467019/80‑1 (MQ=255)
                                                                                                                                                aTGGACCGATGCCGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGACATCCTCGCGATCAAGGGTTATCTCGACAGCCTGCAAcc     <  2:833425/150‑1 (MQ=255)
                                                                                                                                                   gACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGACATCCTCGCGATCAAGGGTTATCTCGACAGCCTGCAACCGGt  <  2:758361/150‑1 (MQ=255)
                                                                                                                                                   gACCGATGACGAGTTCGTCAGGGCGCTGGAGACCGGCATCGGCCGAGAGGGCAAGCACTACTACCCGGCGTATCCCTATACCTCCTAGAGCAAAATGTCGCGCGATGACCTCCTGGCGGACAAGGGTTATCTCGACAGTCTGCAACCggg  >  2:917305/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
ATACCTGACGCGCGCTGCCAACTGTGTGGCTTGCCATACCGTGCCTGGCGGTCAACCCTTCGCCGGTGGTGTTGAGTTCAAGCTGCCGTTCGGCTCGCTGTTTTCGCCCAACATTACCCCGGACCCGCAAACCGGCATCGGTGCATGGACCGATGACGAGTTCGTCAGTGCGCTGCAGACCGGCATTGGCCGTGACGGCAAGCACTACTACCCGGCGTTTCCCTATACCTCCTACAGCAAAATGTCGCGCGATGACATCCTCGCGATCAAGGGTTATCTCGACAGCCTGCAACCGGT  >  NC_002947/2826726‑2827022

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: