Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 978,341 C→A 6.0% L109M (CTG→ATG) ‡ iscR → DNA‑binding transcriptional regulator IscR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947978,3410CA 6.0% 213.0 / 2.3 67L109M (CTG→ATG) ‡iscRDNA‑binding transcriptional regulator IscR
Reads supporting (aligned to +/- strand):  ref base C (34/29);  new base A (2/2);  total (36/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.53e-01

GTCCAGGCGGTGGCTACCAGCTGTCGCGGGGCATGGAAACCATCCAGGTGGCCCAGGTTATCGACGCGGTCAACGAATCGGTCGACGCCACCCGCTGCCAGGGGCTCGGTGACTGCCATGCCGGTGATACCTGCCTGACCCACCACCTGTGGTGCGACCTCAGCCAGCAGATCCATGAATTCCTCAGCGGCATAAGCCTGGCCGACCTCGTCAAGCGCCGTGATGTGCAGGAAGTCGCCCAGCGCCAGGACCTGCGCCGTGTCGCAGGACGAACCGCCCAGCTGGACAAGATTG  >  NC_002947/978195‑978488
                                                                                                                                                  |                                                                                                                                                   
gTCCAGGCGGTGGCTACCAGCTGTCGCGGGGCATGGAAACCATCCAGGTGGCCCAGGTTATCGACGCGGTCAACGAATCGGTCGACGCCACCCGCTGCCAGGGGCTCGGTGACTGCCATGCCGGTGATACCTGCCTGACCCACCACCtgt                                                                                                                                                  <  2:1090737/150‑1 (MQ=255)
 tCCAGGCGGTGGCTCCCAGCTGTCGCGGGGCATGGAAACGATCCAGGTGGCCCAGGTTATCGACGCGGTCAACGAATCGGGCGGCGCCACCCGCTGCCAGGGGCTCGGTGACTGCCATGCTGGTGATACCTGCCTGACCCACCACCtgtg                                                                                                                                                 <  2:600767/150‑1 (MQ=255)
 tCCAGGCGGTGGCTACCAGCTGTCGCGGGGCATGGAAACCATCCAGGTGGCCCAGGTTATCGACGCGGTCAACGAATCGGTCGACGCCACCCGCTGCCAGGGGCTCGGTGACTGCCATGCCGGTGATACCTGCCTGACCCACCACCtgtg                                                                                                                                                 <  1:535109/150‑1 (MQ=255)
 cccAGGCGGTGGCTACCAGCTGTCGCGGGGCATGGAAACCATCCAGGTGGCCCAGGTTATCGACGCGGTCAACGAATCGGTCGACGCCACCCGCTGCCAGGGGCTCGGTGACTGCCATGCCGGTGATACCTGCCTGACCCaccacccgcg                                                                                                                                                 >  2:729281/2‑146 (MQ=255)
  ccAGGCGGTGGCTACCAGCTGTCGCGGGGCATGGAAACCATCCAGGTGGCCCAGGTTATCGACGCGGTCAACGAATCGGTCGACGCCACCCGCTGCCAGGGGCTCGGTGACTGCCATGCCGGTGATACCTGCCTGACCCACCACCTgtgg                                                                                                                                                <  1:843952/150‑1 (MQ=255)
   cAGGCGGTGGCTACCAGCTGTCGCGGGGCATGGAAACCATCCAGGTGGCCCAGGTTATCGACGCGGTCAACGAATCGGTCGACGCCACCCGCTGCCAGGGGCTCGGTGACTGCCATGCCGGTGATACCTGCCTGACCCACCACCTgtggt                                                                                                                                               <  1:1448830/150‑1 (MQ=255)
    aGGCGGTGGCTACCAGCTGTCGCGGGGCATGGAAACCATCCAGGTGGCCCAGGTTATCGACGCGGTCAACGAATCGGTCGACGCCACCCGCTGCCAGGGGCTCGGTGACTGCCATGCCGGTGATACCTGCCTGACCCACCACCTgtggtg                                                                                                                                              >  1:1069892/1‑150 (MQ=255)
     ggCGGTGGCTACCAGCTGTCGCGGGGCATGGAAACCATCCAGGTGGCCCAGGTTATCGACGCGGTCAACGAATCGGTCGACGCCACCCGCTGCCAGGGGCTCGGTGACTGCCATGCCGGTGATACCTGCCTGACCCACCACCTGTGGTGc                                                                                                                                             <  1:1101070/150‑1 (MQ=255)
      gCGGTGGCTACCAGCTGTCGCGGGGCATGGAAACCATCCAGGTGGCCCAGGTTATCGACGCGGTCAACGAATCGGTCGACGCCACCCGCTGCCAGGGGCTCGGTGACTGCCATGCCGGTGATACCTGCCTGACCCACCACCTGTGGTGCg                                                                                                                                            >  1:1233245/1‑150 (MQ=255)
                  aGCTGTCGCGGGGCATGGAAACCATCCAGGTGGCCCAGGTTATCGACGCGGTCAACGAATCGGTCGACGCCACCCGCTGCCAGGGGCTCGGTGACTGCCATGCCGGTGATACCTGCCTGACCCACCACCTGTGGTGCGACCTcagccagc                                                                                                                                <  2:1616811/150‑1 (MQ=255)
                          cgGGGCATGGAAACCATCCAGGTGGCCCAGGTTATCGACGCGGTCAACGAATCGGTCGACGCCACCCGCTGCCAGGGGCTCGGTGACTGCCATGCCGGTGATACCTGCCTGACCCACCACCTGTGGTGCGACCTCAGCCAGCAGATCCAt                                                                                                                        <  2:500630/150‑1 (MQ=255)
                            gggCATGGAAACCATCCAGGTGGCCCAGGTTATCGACGCGGTCAACGAATCGGTCGACGCCACCCGCTGCCAGGGGCTCGGTGACTGCCATGCCGGTGATACCTGCCTGACCCACCACCTGTGGTGCGACCTCAGCCAGCAGATCCATGa                                                                                                                      >  1:952672/1‑150 (MQ=255)
                             ggCATGGAAACCATCCAGGTGGCCCAGGTTATCGACGCGGTCAACGAATCGGTCGACGCCACCCGCTGCCAGGGGCTCGGTGACTGCCATGCCGGTGATACCTGCCTGACCCACCACCTGTGGTGCGACCTCAGCCAGCAGATCCATGaa                                                                                                                     <  1:955548/150‑1 (MQ=255)
                                aTGGAAACCATCCAGGTGGCCCAGGTTATCGACGCGGTCAACGAATCGGTCGACGCCACCCGCTGCCAGGGGCTCGGTGACTGCCATGCCGGTGATACCTGCCTGACCCACCACCTGTGGTGCGACCTCAGCCAGCAGATCCATGAATTc                                                                                                                  <  2:1069892/150‑1 (MQ=255)
                                 tGGAAACCATCCAGGTGGCCCAGGTTATCGACGCGGTCAACGAATCGGTCGACGCCACCCGCTGCCAGGGGCTCGGTGACTGCCATGCCGGTGATACCTGCCTGACCCACCACCTGTGGTGCGACCTCAGCCAGCAGATCCATGAATTcc                                                                                                                 >  1:1526106/1‑150 (MQ=255)
                                  ggAAACCATCCAGGTGGCCCAGGTTATCGACGCGGTCAACGAATCGGTCGACGCCACCCGCTGCCAGGGGCTCGGTGACTGCCATGCCGGTGATACCTGCCTGACCCACCACCTGTGGTGCGACCTCAGCCAGCAGATCCATGAATTCCt                                                                                                                >  2:1064594/1‑150 (MQ=255)
                                              ggTGGCCCAGGTTATCGACGCGGTCAACGAATCGGTCGACGCCACCCGCTGCCAGGGGCTCGGTGACTGCCATGCCGGTGATACCTGCCTGACCCACCACCTgtgg                                                                                                                                                >  2:1522051/1‑106 (MQ=255)
                                              ggTGGCCCAGGTTATCGACGCGGTCAACGAATCGGTCGACGCCACCCGCTGCCAGGGGCTCGGTGACTGCCATGCCGGTGATACCTGCCTGACCCACCACCTgtgg                                                                                                                                                <  1:1522051/106‑1 (MQ=255)
                                               gTGGCCCAGGTTATCGACGCGGTCAACGAATCGGTCGACGCCACCCGCTGCCAGGGGCTCGGTGACTGCCATGCCGGTGATACCTGCCTGACCCACCACCTGTGGTGCGACCTCAGCCAGCAGATCCATGAATTCCTCAGCGGCATAAGc                                                                                                   <  1:1230383/150‑1 (MQ=255)
                                                              gACGCGGTCAACGAATCGGTCGACGCCACCCGCTGCCAGGGGCTCGGTGACTGCCATGCCGGTGATACCTGCCTGACCCACCACCTGTGGTGCGACCTCAGCCAGCAGATCCATGAATTCCTCAGCGGCATAAGCCTGGCCGACCTCGTc                                                                                    <  2:1233245/150‑1 (MQ=255)
                                                              gACGCGGTCAACGAATCGGTCGACGCCACCCGCTGCCAGGGGCTCGGTGACTGCCATGCCGGTGATACCTGCCTGACCCACCACCTGTGGTGCGACCTCAGCCAGCAGATCCATCAATTCCTCAGCGGCATAAGCCTGGCCGACCTCGTc                                                                                    >  2:1335883/1‑150 (MQ=255)
                                                                  cgGTCAACGAATCGGTCGACGCCACCCGCTGCCAGGGGCTCGGTGACTGCCATGCCGGTGATACCTGCCTGACCCACCACCTGTGGTGCGACCTCAGCCAGCAGATCCATGAATTCCTCAGCGGCATAAGCCTGGCCGACCTCGTCAAgc                                                                                <  1:1877916/150‑1 (MQ=255)
                                                                    gTCAACGAATCGGTCGACGCCACCCGCTGCCAGGGGCTCGGTGACTGCCATGCCGGTGATACCTGCCTGACCCACCACCTGTGGTGCGACCTCAGCCAGCAGATCCATGAATTCCTCAGCGGCATAAGCCTGGCCGACCTCGTCAAgcgc                                                                              >  1:751520/1‑150 (MQ=255)
                                                                     tCAACGAATCGGTCGACGCCACCCGCTGCCAGGGGCTCGGTGACTGCCATGCCGGTGATACCTGCCTGACCCACCACCTGTGGTGCGACCTCAGCCAGCAGATCCATGAATTCCTCAGCGGCATAAGCCTGGCCGACCTCGTCAAGCGcc                                                                             >  1:275139/1‑150 (MQ=255)
                                                                      cAACGAATCGGTCGACGCCACCCGCTGCCAGGGGCTCGGTGACTGCCATGCCGGTGATACCTGCCTGACCCACCACCTGTGGTGCGACCTCAGCCAGCAGATCCATGAATTCCTCAGCGGCATAAGCCTGGCCGACCTCGTCAAGCGCCg                                                                            >  1:331792/1‑150 (MQ=255)
                                                                       caCGAGTCGGTCGCCGCCACCGGCTGCCAGGGGCTCGGTGACTGCCCTGCCGGTGATACCTGGCTGACCCACCACCTGTGGTGCGACCTCAGCCAGCAGATCCATGAATGCCTCAGCGGCATAAGCCTGGCCGACCTCGTCAAGCGCCGt                                                                           <  2:301004/149‑1 (MQ=255)
                                                                       aaCGAATCGGTCGACGCCACCCGCTGCCAGGGGCTCGGTGACTGCCATGCCGGTTATACCTGCCTGACCCACCACCTCTGGTGCGACCTCAGCCAGCACATCCATCAATTCCTCAGCGGCATAAGCCTGGCCGACCTCGTCAAGCGCCGt                                                                           >  2:740392/1‑150 (MQ=255)
                                                                       aaCGAATCGGTCGACGCCACCCGCTGCCAGGGGCTCGGTGACTGCCATGCCGGTGATACCTGCCTGACCCACCACCTGTGGTGCGACCTCAGCCAGCAGATCCATGAATTCCTCAGCGGCATAAGCCTGGCCGACCTCGTCAAGCGCCGt                                                                           >  2:92215/1‑150 (MQ=255)
                                                                          gAATCGGTCGACGCCACCCGCTGCCAGGGGCTCGGTGACTGCCATGCCGGTGATACCTGCCTGACCCACCACCTGTGGTGCGACCTCAGCCAGCAGATCCATGAATTCCTCAGCGGCATAAGCCTGGCCGACCTCGTCAAGCGCCGTGAt                                                                        <  1:678613/150‑1 (MQ=255)
                                                                               ggTCGACGCCACCCGCTGCCAGGGGCTCGGTGACTGCCATGCCGGTGATACCTGCCTGACCCACCACCTGTGGTGCGACCTCAGCCAGCAGATCCATGAATTCCTCAGCGGCATAAGCCTGGCCGACCTCGTCAAGCGCCGTGATGTGCa                                                                   >  2:1677264/1‑150 (MQ=255)
                                                                                gTCGACGCCACCCGCTGCCAGGGGCTCGGTGACTGCCATGCCGGTGATACCTGCCTGACCCACCACCTGTGGTGCGACCTCAGCCAGCAGATCCATGAATTCCTCAGCGGCATAAGCCTGGCCGACCTCGTCAAGCGCCGTGATGTGCAg                                                                  <  2:290138/150‑1 (MQ=255)
                                                                                gTCGACGCCACCCGCTGCCAGGGGCTCGGTGACTGCCATGCCGGTGATACCTGCCTGACCCACCACCTGTGGTGCGACCTCAGCCAGCAGATCCATGAATTCCTCAGCGGCATAAGCCTGGCCGACCTCCTCAAGCGCCGTGGTGTTCAg                                                                  <  2:1368256/150‑1 (MQ=255)
                                                                                gTCGACGCCACCCGCTGCCAGGGGCTCGGTGACTGCCATGCCGGTGATACCTGCCTGACCCACCACAGgtggtgggt                                                                                                                                           >  1:246861/1‑74 (MQ=255)
                                                                                gTCGACGCCACCCGCTGCCAGGGGCTCGGTGACTGCCATGCCGGTGATACCTGCCTGACCCACCACAGgtggtgggt                                                                                                                                           <  2:246861/77‑4 (MQ=255)
                                                                                 tCGACGCCACCCGCTGCCAGGGGCTCGGTGACTGCCATGCCGGTGATACCTGCCTGACCCACCACCTGTGGTGCGACCTCAGCCAGCAGATCCATGAATTCCTCAGCGGCATAAGCCTGGCCGACCTc                                                                                       >  2:1727842/1‑128 (MQ=255)
                                                                                 tCGACGCCACCCGCTGCCAGGGGCTCGGTGACTGCCATGCCGGTGATACCTGCCTGACCCACCACCTGTGGTGCGACCTCAGCCAGCAGATCCATGAATTCCTCAGCGGCATAAGCCTGGCCGACCTc                                                                                       <  1:1727842/128‑1 (MQ=255)
                                                                                 tCGACGCCACCCGCTGCCAGGGGCTCGGTGACTGCCATGCCGGTGATACCTGCCTGACCCACCACCTGTGGTGCGACCTCAGCCAGCAGATCCATGAATTCCTCAGCGGCATAAGCCTGGCCGACCTCGTCAAGCGCCGTGATGTGCAgg                                                                 >  2:1884080/1‑150 (MQ=255)
                                                                                 gcGACGCCACCCGCTGCCAGGGGCTCGCTGTCTGCCATGCCGGTGATGCCTGCCTGACCCACCACCTGTGGTGCGACCTCAGCCAGCAGATCCATGAATTCCTCAGCGGCATAAACCTGGCCGACCTCGTCAAGCGCCGTGATGTGCAgg                                                                 <  2:1016693/149‑1 (MQ=255)
                                                                                   gACGCCACCCGCTGCCAGGGGCTCGGTGACTGCCATGCCGGTGATACCTGCCTGACCCACCACCtgtg                                                                                                                                                 >  2:1700271/1‑68 (MQ=255)
                                                                                   gACGCCACCCGCTGCCAGGGGCTCGGTGACTGCCATGCCGGTGATACCTGCCTGACCCACCACCtgtg                                                                                                                                                 <  1:1700271/68‑1 (MQ=255)
                                                                                      gCCACCCGCTGCCAGGGGCTCGGTGACTGCCATGCCGGTGATACCTGCCTGACCCACCACAGgtggtgg                                                                                                                                             <  1:1187263/69‑2 (MQ=255)
                                                                                      gCCACCCGCTGCCAGGGGCTCGGTGACTGCCATGCCGGTGATACCTGCCTGACCCACCACAGgtggtgg                                                                                                                                             >  2:1187263/1‑68 (MQ=255)
                                                                                      gCCACCCGCTGCCAGGGGCTCGGTGACTGCCATGCCGATGATACCTGCCTGACCCACCACCTGTGGTGCGACCTCAGCCAGCAGATCCATGAATTCCTCAGCGGCATAAGCCTGGCCGACCTCGTCAAGCGCCGTGATGTGCAGGAAGTc                                                            >  1:1187787/1‑150 (MQ=255)
                                                                                        cACCCGCTGCCAGGGGCTCGGTGACTGCCATGCCGGTGATACCTGCCTGACCCACCACCTGTGGTGCGACCTCAGCCAGCAGATCCATGAATTCCTCAGCGGCATAAGCCTGGCCGACCTCGTCAAGCGCCGTGATGTGCAGGAAGTCGc                                                          >  1:1677602/1‑150 (MQ=255)
                                                                                          cccGCTGCCAGGGGCTCGGTGACTGCCATGCCGGTGATACCTGCCTGACCCACCACCTGTGGTGCGACCTCAGCCAGCAGATCCATGAATTCCTCAGCGGCATAAGCCTGGCCGACCTCGTCAAGCGCCGTGATGTGCAGGAAGTCGccc                                                        >  1:196917/1‑150 (MQ=255)
                                                                                               tGCCAGGGGCTCGGTGACTGCCATGCCGGTGATACCTGCCTGACCCACCACCTGTGGGGCGCCCTCAGCCAGCAGGTCCATGAATTCCTCAGCGGCATAAGCCTGGCCGACCTCGTCAAGCGCCGTGCTGTACGGGAAGTCGCCCAgcgc                                                   >  1:1038248/1‑150 (MQ=255)
                                                                                                 ccAGGGGCTCGGTGACTGCCATGCCGGTGATACCTGCCTGACCCACCACCTGTGGTGCGACCTCAGCCAGCAGATCCATGAATTCCTCAGCGGCATAAGCCTGGCCGACCTCGTCAAGCGCCGTGATGTGCAGGAAGTCGc                                                          <  1:136683/141‑1 (MQ=255)
                                                                                                 ccAGGGGCTCGGTGACTGCCATGCCGGTGATACCTGCCTGACCCACCACCTGTGGTGCGACCTCAGCCAGCAGATCCATGAATTCCTCAGCGGCATAAGCCTGGCCGACCTCGTCAAGCGCCGTGATGTGCAGGAAGTCGc                                                          >  2:136683/1‑141 (MQ=255)
                                                                                                    ggggCTCGGTGACTGCCATGCCGGTGATACCTGCCTGACCCACCACCTGTGGTGCGACCTCAGCCAGCAGATCCATGAATTCCTCAGCGGCATAAGCCTGGCCGACCTCGTCAAGCGCCGTGATGTGCAGGAAGTCGCCCAGCGCCAGGa                                              <  1:1653699/150‑1 (MQ=255)
                                                                                                      ggCTCGGTGACTGCCATGCCGGTGATACCTGCCTGACCCACCACCTGTGGTGCGACCTCAGCCAGCAGATCCATGAATTCCTCAGCGGCATAAGCCTGGCCGACCTCGTCAAGCGCCGTGATGTGCAGGAAGTCGCCCAGCGCCAGGAcc                                            <  1:1677264/150‑1 (MQ=255)
                                                                                                             tGACTGCCATGCCGGTGATACCTGCCTGACCCACCACCTGTGGTGCGACCTCAGCCAGCAGATCCATGAATTCCTCAGCGGCATAAGCCTGGCCGACCTCGTCAAGCGCCGTGATGTGCAGGAAGTCGCCCAGCGCCAGGACCTGCGCCg                                     >  1:1494220/1‑150 (MQ=255)
                                                                                                              gACTGCCATGCCGGTGATACCTGCCTGACCCACCACCTGTGGTGCGACCTCAGCCAGCAGATCCATGAATTCCTCAGCGGCATAAGCCTGGCCGACCTCGTCAAGCGCCGTGATGTGCAGGAAGTCGCCCAGCGCCAGGACCTGCGCCgt                                    <  1:1587849/150‑1 (MQ=255)
                                                                                                              gACTGCCATGCCGGTGATACCTGCCTGACCCACCACCTGTGGTGCGACCTCAGCCAGCAGATCCATGAATTCCTCAGCGGCATAAGCCTGGCCGACCTCGTCAAGCGCCGTGATGTGCAGGAAGTCGCCCAGCGCCAGGACCTGCGCCgt                                    >  1:1560084/1‑150 (MQ=255)
                                                                                                                 tGCCATGCCGGTGATACCTGCCTGACCCACCACCTGTGGTGCGACCTCAGCCAGCAGATCCATGAATTCCTCAGCGGCATAAGCCTGGCCGACCTCGTCAAGCGCCGTGATGTGCAGGAAGTCGCCCAGCGCCAGGACCTGCGCCGTGTc                                 >  1:350960/1‑150 (MQ=255)
                                                                                                                   ccATGCCGGTGATACCTGCCTGACCCACCACCTGTGGTGCGACCTCAGCCAGCAGATCCATGAATTCCTCAGCGGCATAAGCCTGGCCGACCTCGTCAAGCGCCGTGATGTGCAGGAAGTCGTCCAGCGCCAGGACCTGCGCCGTGTCGc                               <  2:1835822/150‑1 (MQ=255)
                                                                                                                    cATGCCGGTGATACCTGCCTGACCCACCACCTGTGGTGCGACCTCAGCCAGCAGATCCATGAATTCCTCAGCGGCATAAGCCTGGCCGACCTCGTCAAGCGCCGTGATGTGCAGGAAGTCGCCCAGCGCCAGGACCTGCGCCGTGTCGCa                              >  1:1778376/1‑150 (MQ=255)
                                                                                                                         cGGTGATACCTGCCTGACCCACCACCTGTGGTGCGACCTCAGCCAGCAGATCCATGAATTCCTCAGCGGCATAAGCCTGGCCGACCTCGTCAAGCGCCGTGATGTGCAGGAAGTCGCCCAGCGCCAGGACCTGCGCCGTGTCGCAGGACg                         >  2:940636/1‑150 (MQ=255)
                                                                                                                           gTGATACCTGCCTGACCCACCACCTGTGGTGCGACCTCAGCCAGCAGATCCATGAATTCCTCAGCGGCATAAGCCTGGCCGACCTCGTCAAGCGCCGTGATGTGCAGGAAg                                                              >  2:592655/1‑111 (MQ=255)
                                                                                                                           gTGATACCTGCCTGACCCACCACCTGTGGTGCGACCTCAGCCAGCAGATCCATGAATTCCTCAGCGGCATAAGCCTGGCCGACCTCGTCAAGCGCCGTGATGTGCAGGAAg                                                              <  1:592655/111‑1 (MQ=255)
                                                                                                                                aCCTGCCTGACCCACCACCTGTGGTGCGACCTCAGCCAGCAGATCCATGAATTCCTCAGCGGCATAAGCCTGGCCGACCTCGTCAAGCGCCGTGATGTGCAGGAAGTCGCCCAGCGCCAGGACCTGCGCCGTGTCGCAGGACGAACCGcc                  >  1:1449575/1‑150 (MQ=255)
                                                                                                                                  ctgcctgACCCACCACCTGTGGTGCGACCTCAGCCAGCAGATCCATGAATTCCTCAGCGGCATAAGCCTGGCCGACCTCGTCAAGCGCCGTGATGTGCAGGAAGTCGCCCAGCGCCAGGACCTGCGCCGTGTCGCAGGACGAACCGCCCa                >  1:263133/1‑150 (MQ=255)
                                                                                                                                       tgACCCACCACCTGTGGTGCGACCTCAGCCAGCAGATCCATGAATTCCTCAGCGGCATAAGCCTGGCCGACCTCGTCAAGCGCCGTGATGTGCAGGAAGTCGCCCAGCGCCAGGACCTGCGCCGTGTCGCAGGACGAACCGCCCAGCTgg           >  1:1863191/1‑150 (MQ=255)
                                                                                                                                        gACCCACCACCTGTGGTGCGACCTCAGCCAGCAGATCCATGAATTCCTCAGCGGCATAAGCCTGGCCGACCTCGTCAAGCGCCGTGATGTGCAGGAAGTCGCCCAGCGCCAGGACCTGCGCCGTGTCGCAGGACGAACCGCCCAGCTGGa          >  1:867186/1‑150 (MQ=255)
                                                                                                                                           ccaccaCCTGTGGTGCGACCTCAGCCAGCCGATCCATGAATTCCTCAGCGGCATAAGCCTGGCCGACCTCGTCAAGCGCCGTGATGTGCAGGAAAGCGCCCAGCGCCAGGACCTGCGCCGTGTTGAAGGACGAACCGCCCCGCTGAACaa       >  1:218420/1‑150 (MQ=255)
                                                                                                                                           ccaccaCCTGTGGTGCGACCTCAGCCAGCAGATCCATGAATTCCTCAGCGGCATAAGCCTGGCCGACCTCGTCAAGCGCCGTGAt                                                                        <  1:914475/85‑1 (MQ=255)
                                                                                                                                           ccaccaCCTGTGGTGCGACCTCAGCCAGCAGATCCATGAATTCCTCAGCGGCATAAGCCTGGCCGACCTCGTCAAGCGCCGTGAt                                                                        >  2:914475/1‑85 (MQ=255)
                                                                                                                                           ccaccaCCTGTGGTGCGACCTCAGCCAGCAGATCCATGAATTCCTCAGCGGCATAAGCCTGGCCGACCTCGTCAAGCGCCGTGATGTGCAGGAAGTCGCCCAGCGCCAGGACCTGCGCCGTGTCGCAGGACGAACCGCCCAGCTGGACaa       <  1:214718/150‑1 (MQ=255)
                                                                                                                                            caccacCTGTGGTGCGACCTCAGCCAGCAGATCCATGAATTCCTCAGCGGCATAAGCCTGGCCGACCTCGTCAAGCGCCGTGATGTGCAGGAAGTCGCCCAGCGCCAGGACCTGCGCCGTGTCGCAGGACGAACCGCCCAGCTGGACAAg      >  2:824768/1‑150 (MQ=255)
                                                                                                                                             accaccTGTGGTGCGACCTCAGCCAGCAGATCCATGAATTCCTCAGCGGCATAAGCCTGGCCGACCTCGTCAAGCGCCGTGATGTGCAGGAAGTCGCCCAGCGCCAGGACCTGCGCCGTGTCGCAGGACGAACCGCCCAGCTGGACAAGa     <  1:824768/150‑1 (MQ=255)
                                                                                                                                               caccTGTGGTGCGACCTCAGCCAGCAGATCCATGAATTCCTCAGCGGCATAAGCCTGGCCGACCTCGTCAAGCGCCGTGATGTGCAGGAAGTCGCCCAGCGCCAGGACCTGCGCCGTGTCGCAGGACGAACCGCCCAGCTGGACAAGAtt   <  2:1449575/150‑1 (MQ=255)
                                                                                                                                                accTGTGGTGCGACCTCAGCCAGCAGATCCATGAATTCCTCAGCGGCATAAGCCTGGCCGACCTCGTCAAGCGCCGTGATGTGCAGGAAGTCGCCCAGCGCCCGGACCTGCGCCGTGTCGCAGGACGAACCGCCCAGCTGGACAAGATTg  <  1:1884080/150‑1 (MQ=255)
                                                                                                                                                accTGTGGTGCGACCTCAGCCAGCAGATCCATGAATTCCTCAGCGGCATAAGCCTGGCCGACCTCGTCAAGCGCCGTGATGTGCAGGAAGTCGCCCAGCGCCAGGACCTGCGCCGTGTCGCAGGACGAACCGCCCAGCTGGACAAGATTg  >  2:1876716/1‑150 (MQ=255)
                                                                                                                                                accTGTGGTGCGACCTCAGCCAGCAGATCCATGAATTCCTCAGCGGCATAAGCCTGGCCGACCTCGTCAAGCGCCGTGATGTGCAGAAAGTCGCCCAGCGCAAGGACCTCCGCCGTGTCGCAGGCGGAACCGCCCAGCTGGACAAGATTg  >  1:676369/1‑150 (MQ=255)
                                                                                                                                                  cTGTGGTGCGACCTCAGCCAGCAGATCCATGAATTCCTCAGCGGCATAAGCCTGGCCGACCTCGTCAAGCGCCGTGATGTGCAGGAAGTCGCCCAGCGCCAGGACCTGCGCCGTGTCGCAGGACGAACCGCCCAGCTGGACAAGAtt   <  2:1147158/147‑1 (MQ=255)
                                                                                                                                                  cTGTGGTGCGACCTCAGCCAGCAGATCCATGAATTCCTCAGCGGCATAAGCCTGGCCGACCTCGTCAAGCGCCGTGATGTGCAGGAAGTCGCCCAGCGCCAGGACCTGCGCCGTGTCGCAGGACGAACCGCCCAGCTGGACAAGAtt   >  1:1147158/1‑147 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GTCCAGGCGGTGGCTACCAGCTGTCGCGGGGCATGGAAACCATCCAGGTGGCCCAGGTTATCGACGCGGTCAACGAATCGGTCGACGCCACCCGCTGCCAGGGGCTCGGTGACTGCCATGCCGGTGATACCTGCCTGACCCACCACCTGTGGTGCGACCTCAGCCAGCAGATCCATGAATTCCTCAGCGGCATAAGCCTGGCCGACCTCGTCAAGCGCCGTGATGTGCAGGAAGTCGCCCAGCGCCAGGACCTGCGCCGTGTCGCAGGACGAACCGCCCAGCTGGACAAGATTG  >  NC_002947/978195‑978488

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