Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,639,373 G→C 6.4% V115L (GTG→CTG)  PP_3206 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,639,3730GC 6.4% 298.0 / 9.3 93V115L (GTG→CTG) PP_3206hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (40/47);  new base C (2/4);  total (42/51)
Fisher's exact test for biased strand distribution p-value = 6.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

TGCTAGCTGACGGTGTCGATGTGGTCTTCCATCTCGTCAGCGTGCCTGGCGGCGCTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCATGGTGGAGATCGCGTTGCAGGATCTTGCCCGCCGTGG  >  NC_002947/3639227‑3639518
                                                                                                                                                  |                                                                                                                                                 
tGCTAGCTGACGGTGTCGATGTGGTCTTCCATCTCGTCAGCGTGCCTGGCGGCGCTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGc                                                                                                                                                <  1:1794024/150‑1 (MQ=255)
tGCTAGCTGACGGTGTCGATGTGGTCTTCCATCTCGTCAGCGTGCCTGGCGGCGCTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGc                                                                                                                                                >  1:14387/1‑150 (MQ=255)
  cTAGCTGACGGTGTCGATGTGGTCTTCCATCTCGTCAGCGTGCCTGGCGGCGCTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTg                                                                                                                                              <  1:1031027/150‑1 (MQ=255)
   tAGCTGACGGTGTCGATGTGGTCTTCCATCTCGTCAGCGTGCCTGGCGGCGCTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGCTGCTgg                                                                                                                                             >  1:166063/1‑150 (MQ=255)
      cTGACGGTGTCGATGTGGTCTTCCATCTCGTCAGCGTGCCTGGTGCCGCGTCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGgtgt                                                                                                                                          <  1:205717/150‑1 (MQ=255)
      cTGACGGTGTCGATGTGGTCTTCCATCTCGTCAGCGTGCCTGGCGGCGCTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGgtgt                                                                                                                                          <  1:246830/150‑1 (MQ=255)
      cTGACGGTGTCGATGTGGTCTTCCATCTCGTCAGCGTGCCTGGCGGCGCTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGgt                                                                                                                                            >  1:1131572/1‑148 (MQ=255)
      cTGACGGTGTCGATGTGGTCTTCCATCTCGTCAGCGTGCCTGGCGGCGCTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGgt                                                                                                                                            <  2:1131572/148‑1 (MQ=255)
        gACGGTGTCGATGTGGTCTTCCATCTCGTCAGCGTGCCTGGCGGCGCTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTAc                                                                                                                                        <  1:1504400/150‑1 (MQ=255)
        gACGGTGTCGATGTGGTCTTCCATCTCGTCAGCGTGCCTGGCGGCGCTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTAc                                                                                                                                        <  2:1700767/150‑1 (MQ=255)
               tCGATGTGGTCTTCCATCTCGTCAGCGTGCCTGGCGGCGCTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCcagca                                                                                                                                 <  1:740698/150‑1 (MQ=255)
                 gATGTGGTCTTCCATCTCGTCAGCGTGCCTGGCGGCGCTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTAc                                                                                                                                        <  1:1016336/141‑1 (MQ=255)
                 gATGTGGTCTTCCATCTCGTCAGCGTGCCTGGCGGCGCTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTAc                                                                                                                                        >  2:1016336/1‑141 (MQ=255)
                 gATGTGGTCTTCCATCTCGTCAGCGTGCCTGGCGGCGCTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGCTGCTGGCGTACACcagcacc                                                                                                                               <  1:1432223/150‑3 (MQ=255)
                  aTGTGGTCTTCCATCTCGTCAGCGTGCCTGGCGGCGCTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAGCAGCg                                                                                                                              >  2:1525860/1‑150 (MQ=255)
                   tgtgGTCTTCCATCTCGTCAGCGTGCCTGGCGGCGCTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGCTGCTGGCGTACACCAGCACCgg                                                                                                                             <  2:166063/150‑2 (MQ=255)
                   tgtgGTCTTCCATCTCGTCAGCGTGCCTGGCGGCGCTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGCTGCTGGCGTACACCAGCACCgg                                                                                                                             <  2:1108875/150‑2 (MQ=255)
                       gTCTTCCATCTCGTCAGCGTGCCTGGCGGCGCTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGTCAGCAGCGTcgcg                                                                                                                         >  2:1230931/1‑150 (MQ=255)
                       gTCTTCCATCTCGTCAGCGTGCCTGGCGGCGCTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAGCAGCGTcgcg                                                                                                                         >  2:1076008/1‑150 (MQ=255)
                        tCTTCCATCTCGTCAGCGTGCCTGGCGGCGCTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAGCAGCGTCGCgg                                                                                                                        >  2:1155827/1‑150 (MQ=255)
                         cTTCCATCTCGTCAGCGTGCCTGGCGGCGCTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTTTTCAACCCGTTGAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAGCAGCGTCGCGgt                                                                                                                       <  2:1606676/150‑1 (MQ=255)
                         cTTCCATCTCGTCAGCGTGCCTGGCGGCGCTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAGCAGCGTCGCGgt                                                                                                                       <  2:1643895/150‑1 (MQ=255)
                          ttCCATCTCGTCAGCGTGCCTGGCGGCGCTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGc                                                                                                                                      >  2:1792232/1‑134 (MQ=255)
                          ttCCATCTCGTCAGCGTGCCTGGCGGCGCTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGc                                                                                                                                      <  1:1792232/134‑1 (MQ=255)
                                 tcGTCAGCGTGCCTGGCGGCGCTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAGCAGCGTCGCGGTGTAcggcg                                                                                                               <  2:1429193/150‑1 (MQ=255)
                                   gTCAGCGTGCCTGGCGGCGCTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCATCTGTCCACCGGTGCTGGTGTACGCCAGCAGCGTCGCGGTGTACggcggc                                                                                                             >  2:999937/1‑150 (MQ=255)
                                    tCAGCGTGCCTGGCGGCGCTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAGCAGCGTCGCGGTGTACGGCGGCa                                                                                                            >  2:534046/1‑150 (MQ=255)
                                        cGTGCCTGGCGGCGCTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAGCAGCGTCGCGGTGTACGGCGGCAAGtt                                                                                                        <  1:1204947/150‑1 (MQ=255)
                                          tGCCTGGCGGCGCTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCAGTGCTGGTGTACGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGc                                                                                                      <  1:728137/150‑1 (MQ=255)
                                              tGGCGGCGCTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGAc                                                                                                  >  2:1475123/1‑150 (MQ=255)
                                                gcggcgCTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGAcgc                                                                                                <  1:1306631/150‑1 (MQ=255)
                                                gcggcgCTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGAcgc                                                                                                <  1:197337/150‑1 (MQ=255)
                                                   gcgcTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCa                                                                                             <  2:1236083/150‑1 (MQ=255)
                                                    cgcTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCAt                                                                                            <  2:1236531/150‑1 (MQ=255)
                                                       tgacGAGATGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGa                                                                                         <  2:1672943/147‑1 (MQ=255)
                                                           gagaCGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGTGCATGGACGaa                                                                                     <  2:117679/150‑1 (MQ=255)
                                                           gagaCGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGaa                                                                                     <  1:844283/150‑1 (MQ=255)
                                                                gCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGgcca                                                                                >  1:1206084/1‑150 (MQ=255)
                                                                gCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGgcca                                                                                >  1:868234/1‑150 (MQ=255)
                                                                   gTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGgccagcc                                                                             >  1:1255507/1‑150 (MQ=255)
                                                                    tACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGgccagcca                                                                            <  1:1053550/150‑1 (MQ=255)
                                                                           tgggTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTcgc                                                                     <  2:141210/149‑1 (MQ=255)
                                                                           ggggTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTcgc                                                                     >  2:645535/1‑150 (MQ=255)
                                                                           ggggTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAACGTcgc                                                                     <  1:1525860/150‑1 (MQ=255)
                                                                               ttACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAGCAGCGTCGCGgtg                                                                                                                      <  1:643336/97‑1 (MQ=255)
                                                                               ttACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAGCAGCGTCGCGgtg                                                                                                                      >  2:643336/1‑97 (MQ=255)
                                                                               ttACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTcgccgcc                                                                 >  2:1251080/1‑149 (MQ=255)
                                                                               ttACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCa                                                                 <  2:868234/150‑1 (MQ=255)
                                                                                tACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAg                                                                <  1:1369046/150‑1 (MQ=255)
                                                                                 aCCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAGCAGCGTCGCGgt                                                                                                                       >  1:1350224/1‑94 (MQ=255)
                                                                                 aCCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAGCAGCGTCGCGgt                                                                                                                       <  2:1350224/94‑1 (MQ=255)
                                                                                 aCCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGc                                                               <  2:1255507/150‑1 (MQ=255)
                                                                                  ccAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAGCAGCGTCGCGGTGCTAGAAACGATGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCt                                                              >  2:528267/1‑150 (MQ=255)
                                                                                   cAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTg                                                             >  1:1394887/1‑150 (MQ=255)
                                                                                   cAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAGCAGCGGCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTg                                                             >  2:59502/1‑150 (MQ=255)
                                                                                      gTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTg                                                                                                                                              >  1:1080080/1‑66 (MQ=255)
                                                                                      gTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTg                                                                                                                                              <  2:1080080/66‑1 (MQ=255)
                                                                                       tAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCt                                                         <  2:14387/150‑1 (MQ=255)
                                                                                            cTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATgcg                                                    >  2:391133/1‑150 (MQ=255)
                                                                                            cTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATgcg                                                    >  1:160541/1‑150 (MQ=255)
                                                                                            cTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGCTGCTGGCGTAcac                                                                                                                                      <  1:1631334/68‑3 (MQ=255)
                                                                                            cTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGCTGCTGGCGTAcac                                                                                                                                      >  2:1631334/1‑66 (MQ=255)
                                                                                              gcaagcaaGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATgcggc                                                  <  2:1238838/150‑1 (MQ=255)
                                                                                              gcaagcaaGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATgcggc                                                  <  1:1132069/150‑1 (MQ=255)
                                                                                                 agcaagcCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGca                                               >  1:904728/1‑150 (MQ=255)
                                                                                                 agcaagcCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGca                                               <  1:1475123/150‑1 (MQ=255)
                                                                                                 agcaagcCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCATATGCTGCGca                                               >  1:701213/1‑150 (MQ=255)
                                                                                                     agcCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACaaa                                           <  2:487522/150‑1 (MQ=255)
                                                                                                      gcCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAAc                                          <  1:756104/150‑1 (MQ=255)
                                                                                                          ggAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCAt                                      <  1:1081597/150‑1 (MQ=255)
                                                                                                             gctgctCGACCACCTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCAtggt                                   <  1:391133/150‑1 (MQ=255)
                                                                                                             gctgctCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCAtggt                                   >  1:742008/1‑150 (MQ=255)
                                                                                                              ctgctCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCAtggtg                                  <  1:1155827/150‑1 (MQ=255)
                                                                                                              ctgctCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCAtggtg                                  <  1:1094186/150‑1 (MQ=255)
                                                                                                                gctCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCATGGTGga                                >  1:1033188/1‑150 (MQ=255)
                                                                                                                gctCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCATGGTGga                                <  1:1655735/150‑1 (MQ=255)
                                                                                                                 ctCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCATGGTGgag                               >  2:1169958/1‑150 (MQ=255)
                                                                                                                    gACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCATGGTGGAGATc                            >  2:603035/1‑150 (MQ=255)
                                                                                                                     aCCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCATGGTGGAGATcg                           >  1:909803/1‑150 (MQ=255)
                                                                                                                       cAACTGCGTAATCCTCTCTGTCCACCGGTGCTGGTGTACGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCATGGTGGAGATcgcg                         >  1:105609/1‑150 (MQ=255)
                                                                                                                       cAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCCTGGTGGAGATcgcg                         >  2:595110/1‑150 (MQ=255)
                                                                                                                        aaCTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGgcca                                                                                <  2:1689201/94‑1 (MQ=255)
                                                                                                                        aaCTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGgcca                                                                                >  1:1689201/1‑94 (MQ=255)
                                                                                                                          cTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCATGGTGGAGATCGCGTTg                      >  2:279535/1‑150 (MQ=255)
                                                                                                                          cTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCAATCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCATGGTGGAGATCGCGTTg                      <  2:1112614/150‑1 (MQ=255)
                                                                                                                            gCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCATGGTGGAGATCGCGTTGCa                    >  1:168861/1‑150 (MQ=255)
                                                                                                                               tAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCATGGTGGAGATCGCGTTGCAGGa                 <  1:645535/150‑1 (MQ=255)
                                                                                                                                 aTCCCCTCTGTCCACCGGTGCTGGTGTACGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCATGGTGGAGATCGCGTTGCAGGATc               >  2:224512/1‑150 (MQ=255)
                                                                                                                                 aTCCCCTCTGTCCACCGGTGCTGGTGTACGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCATGGTGGAGATCGCGTTGCAGGATc               <  2:672362/150‑1 (MQ=255)
                                                                                                                                    cccTCTGTCCACCGGTGCTGGTGTACGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAgg                                                                                   <  2:312297/79‑1 (MQ=255)
                                                                                                                                    cccTCTGTCCACCGGTGCTGGTGTACGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAgg                                                                                   >  1:312297/1‑79 (MQ=255)
                                                                                                                                           tCCACCGGTGCTGGTGTACGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCATGGTGGAGATCGCGTTGCAGGATCTTGCccgccg     >  2:1865436/1‑150 (MQ=255)
                                                                                                                                              aCCGGTGCTGGTGTACGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCATGGTGGAGATCGCGTTGCAGGATCTTGCCCGCCGTgg  >  1:685861/1‑150 (MQ=255)
                                                                                                                                                  gTGCTGGTGTACGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCATGGTGGAGATcgc                          >  1:1764870/1‑122 (MQ=255)
                                                                                                                                                  gTGCTGGTGTACGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCATGGTGGAGATcgc                          <  2:1764870/122‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TGCTAGCTGACGGTGTCGATGTGGTCTTCCATCTCGTCAGCGTGCCTGGCGGCGCTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTACGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCATGGTGGAGATCGCGTTGCAGGATCTTGCCCGCCGTGG  >  NC_002947/3639227‑3639518

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