Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,208,214 T→C 5.7% I64T (ATC→ACC)  ttdT → L‑tartrate/succinate antiporte

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,208,2140TC 5.7% 260.8 / 5.1 89I64T (ATC→ACC) ttdTL‑tartrate/succinate antiporte
Reads supporting (aligned to +/- strand):  ref base T (38/45);  new base C (3/2);  total (42/47)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.15e-01

GCGGTCATCGCCATTATTGCTCTACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTAATTCGGTTATCTGGCTGATTTTCGCCGCCTTTATGTTTGGCA  >  NC_000913/3208066‑3208363
                                                                                                                                                    |                                                                                                                                                     
gCGGTCATCGCGATTATTGCTCTACTTCCGGTTCCAGCAGGGCTCGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGAACGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCcatc                                                                                                                                                      >  2:303603/1‑150 (MQ=255)
gCGGTCATCGCCATTATTGCTCTACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCcatc                                                                                                                                                      <  1:1010144/150‑1 (MQ=255)
   gTCATCGCCATTATTGCTCTACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCatcatc                                                                                                                                                   >  2:1287755/1‑150 (MQ=255)
    tCATCGCCATTATTGCTCTACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTTGTGGCGATGATGGATATACCCACCATCg                                                                                                                                                  <  1:1389150/150‑1 (MQ=255)
      aTCGCCATTATTGCTCTACTTCAGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGAGCGCTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGATGGATATACCCACCATCGcc                                                                                                                                                <  1:828755/150‑1 (MQ=255)
       tCGCCATTATTGCTCTACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCa                                                                                                                                               <  1:337867/150‑1 (MQ=255)
        cGCCATTATTGCTCTACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCAt                                                                                                                                              <  2:1395064/150‑1 (MQ=255)
        cGCCATTATTGCTCTACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCACCAt                                                                                                                                              <  2:476350/150‑1 (MQ=255)
         gCCATTATTGCTCTACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATa                                                                                                                                             >  2:642804/1‑150 (MQ=255)
         gCCATTATTGCTCTACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGATTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGAGGGTATATCCATCATCGCCATa                                                                                                                                             >  2:647392/1‑150 (MQ=255)
           cATTATTGCTCTACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATAct                                                                                                                                           <  2:75869/150‑1 (MQ=255)
            attattGCTCTACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCGCTCCTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATActc                                                                                                                                          <  1:309692/150‑1 (MQ=255)
                ttGCTCTACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACtctctc                                                                                                                                      <  2:928450/150‑1 (MQ=255)
                   ctctACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCt                                                                                                                                   <  2:826072/150‑1 (MQ=255)
                   ctctACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCt                                                                                                                                   >  1:672716/1‑150 (MQ=255)
                             ggTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGtt                                                                                                                         >  1:152227/1‑150 (MQ=255)
                              gTTCCCGCAGGGCTGGAGAACCATTCCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGtta                                                                                                                        >  1:496791/1‑149 (MQ=255)
                                tCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAg                                                                                                                      >  2:791684/1‑150 (MQ=255)
                                 cccGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTAATGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCCATGGTGGGTATATACATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGc                                                                                                                     >  1:747222/1‑150 (MQ=255)
                                  ccGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGcc                                                                                                                    >  2:947022/1‑150 (MQ=255)
                                  ccGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGATGTTCAGcc                                                                                                                    >  2:1308344/1‑150 (MQ=255)
                                  ccGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGATGGACTGATCCTCTATCCCGTGCCGGGTGCCGTGTTGGCGATGGTGGGTATATACATAATCTCCATACTCTCTCACTGGATGCTGTTAAGcc                                                                                                                    >  2:68410/1‑150 (MQ=255)
                                   cGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGCACCTTGGACTGATCCTAGAACCCCTGCCGCGTGCCGTGGTGGCGATGGTGGGTATAGCCAAAATCGCCAAACTCACACCCTGGCAGCGGTTCAGccc                                                                                                                   >  1:1558786/1‑150 (MQ=255)
                                   cGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGGGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGccc                                                                                                                   >  1:807734/1‑150 (MQ=255)
                                     cAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCgg                                                                                                                 <  1:570036/150‑1 (MQ=255)
                                     cAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGGGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCgg                                                                                                                 <  2:931357/150‑1 (MQ=255)
                                        ggCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGagc                                                                                                              <  1:1345622/150‑1 (MQ=255)
                                         gCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTCGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCCGCCCGGagct                                                                                                             >  1:221091/1‑149 (MQ=255)
                                           tGGAGAACCATACCTGGCTCTACTTTGCAGTTATTACTGGCGTGATCGTTGGACTGATAATCGAACCCGTGCCCGGTGCCGTGGTGGCGATGGGGGGTATATCCATCATAGCAATTCTCTCACCCTGGCTGATTATCAGAACCGagcagc                                                                                                           >  2:639175/1‑150 (MQ=255)
                                            ggAGAACCATACCTGGCTCTACTTTGCCGTTTTTCCTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAgct                                                                                                          <  1:290533/150‑1 (MQ=255)
                                            ggAGAACCATACCTGGCGCTACTTAGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAgct                                                                                                          <  2:1131100/150‑1 (MQ=255)
                                              agaACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAgctcg                                                                                                        <  1:1708133/150‑1 (MQ=255)
                                               gaACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAgctcgc                                                                                                       >  2:1715089/1‑150 (MQ=255)
                                                    aTACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGc                                                                                                  >  2:1519927/1‑150 (MQ=255)
                                                    aTACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGc                                                                                                  >  2:1221778/1‑150 (MQ=255)
                                                     tACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGcc                                                                                                 >  1:28898/1‑150 (MQ=255)
                                                        cTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCCTACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAgg                                                                                              >  2:929874/1‑150 (MQ=255)
                                                           gCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCtt                                                                                           >  1:1040885/1‑150 (MQ=255)
                                                           gCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCAGTGGTGGCGATGGTGGGTATATACATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCtt                                                                                           >  1:1055087/1‑150 (MQ=255)
                                                            ctctACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCttt                                                                                          >  2:640875/1‑150 (MQ=255)
                                                            ctctACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCttt                                                                                          >  1:1180787/1‑150 (MQ=255)
                                                               tACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTTCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTaaa                                                                                       <  1:870303/150‑1 (MQ=255)
                                                                aCTTTGCCGTTTTTACTGGCGTGACCGCTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAAt                                                                                      <  2:1465559/150‑1 (MQ=255)
                                                                 cTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAAtt                                                                                     >  1:1088625/1‑150 (MQ=255)
                                                                 cTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAAtt                                                                                     >  1:1566858/1‑150 (MQ=255)
                                                                   ttGCCGGTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCa                                                                                   <  1:1646923/150‑1 (MQ=255)
                                                                    tGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCAc                                                                                  <  1:947022/150‑1 (MQ=255)
                                                                      ccGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTg                                                                                <  2:56318/150‑1 (MQ=255)
                                                                      ccGTTTTTACTGGCGGGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTg                                                                                <  1:243537/150‑1 (MQ=255)
                                                                             tACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTTTATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATc                                                                         <  2:218005/150‑1 (MQ=255)
                                                                             tACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATc                                                                         >  2:173140/1‑150 (MQ=255)
                                                                               cTGGCGTGATCGTTTGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATccc                                                                       <  2:398492/150‑1 (MQ=255)
                                                                                   cGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCtctc                                                                   <  2:450301/150‑1 (MQ=255)
                                                                                    gTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCg                                                                  >  1:69470/1‑150 (MQ=255)
                                                                                    gTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGATCAGCCCGGAGCAGCTCGCTCATCCAGGCTTTAAATTCACTGCAAAATCCCTCTCg                                                                  <  1:303603/150‑1 (MQ=255)
                                                                                      gATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTg                                                                <  2:1476603/150‑1 (MQ=255)
                                                                                        tCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTggg                                                              <  2:1088625/150‑1 (MQ=255)
                                                                                           ttGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCg                                                           <  2:672716/150‑1 (MQ=255)
                                                                                                cTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTcc                                                      <  2:425158/150‑1 (MQ=255)
                                                                                                cTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCAGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTcc                                                      >  1:853676/1‑150 (MQ=255)
                                                                                                  gATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCgg                                                    >  1:783131/1‑150 (MQ=255)
                                                                                                  gATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCgg                                                    >  1:1272985/1‑150 (MQ=255)
                                                                                                      cTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGtttt                                                <  1:1287755/150‑1 (MQ=255)
                                                                                                        cGAACCCGTGCCGTGTGCCTTGGTGGCGATGGTGGGTATATCAATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGTCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGtttttg                                              <  1:68410/150‑2 (MQ=255)
                                                                                                        cGAACCCATGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTc                                              <  1:944928/150‑1 (MQ=255)
                                                                                                         gAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAATATCCCTCTCGTGGGCCGTTTCCGGTTTTTCt                                             <  1:647392/150‑1 (MQ=255)
                                                                                                         gAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCt                                             <  1:173140/150‑1 (MQ=255)
                                                                                                          aaCCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTa                                            <  2:13112/150‑1 (MQ=255)
                                                                                                          aaCCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTa                                            >  2:474137/1‑150 (MQ=255)
                                                                                                          aaCCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTa                                            <  2:1180787/150‑1 (MQ=255)
                                                                                                          aaCCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTa                                            <  2:358611/150‑1 (MQ=255)
                                                                                                           aCCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTaa                                           <  1:1191545/150‑1 (MQ=255)
                                                                                                              cGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTAATTc                                        <  2:982317/150‑1 (MQ=255)
                                                                                                                tGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTAATTCgg                                      <  2:1566858/150‑1 (MQ=255)
                                                                                                                        gCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTACTTCGGTAATCTgg                              >  2:1142645/1‑150 (MQ=255)
                                                                                                                         ccGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTAATTCGGTTATCTGGc                             >  1:992913/1‑150 (MQ=255)
                                                                                                                           gtggtgGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTAATTCGGTTATCTGGCTg                           >  2:1407455/1‑150 (MQ=255)
                                                                                                                            gtatggCGATGATGGATATACCCACCATCGCCATACTCTCTCCCTGGATGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGGTTTAAATTCACTGCAAAATCCCTCACGTGGGCCGATTCCGGTGTTTCTAATTCGGTTATCTGGCTGa                          >  1:1107589/4‑150 (MQ=255)
                                                                                                                               tggCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTAATTCGGTTATCTGGCTGAttt                       <  1:1715089/150‑1 (MQ=255)
                                                                                                                               tggCGATGATGGATATACCCACCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCCGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTAATTCGGTTATCTGGCTGAttt                       >  2:177635/1‑150 (MQ=255)
                                                                                                                               tggCGATGATGGATATACCCACCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTAATTCGGTTATCTGGCTGAttt                       >  2:928860/1‑150 (MQ=255)
                                                                                                                                  cGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTAATTCGGTTATCTGGCTGATTTTcg                    >  2:57763/1‑150 (MQ=255)
                                                                                                                                  cGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTAATTCGGTTATCTGGCTGATTTTcg                    <  1:1519927/150‑1 (MQ=255)
                                                                                                                                  cGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTAATTCGGTTATCTGGCTGATTTTcg                    >  1:1321846/1‑150 (MQ=255)
                                                                                                                                  cGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTAATTCGGTTATCTGGCTGATTTTcg                    <  1:1221778/150‑1 (MQ=255)
                                                                                                                                    aTGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCATCCAGGCTTTAAATTCACTGCAAAATCCCTGTCGTGGGCCGTTTCCGGTTTTTCTAATTCGGGTATCTGGCTGATTTTcgcc                  <  2:152227/150‑1 (MQ=255)
                                                                                                                                      ggtggGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGTGCAGCTCGCTCAGCCAGGCTTTATATTCACTTCAAAATCCCTCGCGTGGGCGGTGTCCGGTTTGTCTAATTCGGTTATCTGTCTGATTTTcgccgc                >  1:311865/1‑150 (MQ=255)
                                                                                                                                       gtggGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTAATTCGGTTATCTGGCTGATTTTCgccgcc               <  1:1131293/150‑1 (MQ=255)
                                                                                                                                          ggTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAATTCCCTCTCGTGGGCCGTTTCCGGTTTTTCTAATTCGGTTATCTGGCTGATTTTCGCCGCCttt            <  1:219968/150‑1 (MQ=255)
                                                                                                                                          ggTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTAATTCGGTTATCTGGCTGATTTTCGCCGCCttt            >  1:1365182/1‑150 (MQ=255)
                                                                                                                                           gTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTAATTCGGTTATCTGGCTGATTTTCGCCGCCTTTa           <  2:783131/150‑1 (MQ=255)
                                                                                                                                              tatCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTAATTCGGTTATCTGGCTGATTTTCGCCGCCTTTATGt        >  1:346615/1‑150 (MQ=255)
                                                                                                                                                  catcatCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTAATTCGGTTATCTGGCTGATTTTCGCCGCCTTTATGTTTgg    <  1:1475890/150‑1 (MQ=255)
                                                                                                                                                  catcatCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTCGGCCGTTTCCGGTTTTTCTAATTCGGTTATCTGGCTGATTTTCGCCGCCTTTATGTTTgg    >  1:1024833/1‑150 (MQ=255)
                                                                                                                                                    tcatcGCCATACTCGCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTGTTCTAATTCGGTTATCTGGCTGATTTTCGCCGCCTTTATGTTTGGca  <  2:1628136/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCGGTCATCGCCATTATTGCTCTACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTAATTCGGTTATCTGGCTGATTTTCGCCGCCTTTATGTTTGGCA  >  NC_000913/3208066‑3208363

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