Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,143,980 C→A 8.5% G62V (GGC→GTC)  hybB ← putative hydrogenase 2 cytochrome b type component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,143,9800CA 8.5% 263.7 / 3.5 97G62V (GGC→GTC) hybBputative hydrogenase 2 cytochrome b type component
Reads supporting (aligned to +/- strand):  ref base C (51/36);  new base A (2/8);  total (53/44)
Fisher's exact test for biased strand distribution p-value = 3.97e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACCCACGTCGATAGTGATCGACAAGCCACCCAGTGAGTAACCAAACAGACTCGCCAACAGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGAACACCAGACGCTTCACAATCAGGAGCATACAGATGACGATTAACGGTCCAAAAATCATGACCGGTTTACTGATG  >  NC_000913/3143831‑3144129
                                                                                                                                                     |                                                                                                                                                     
taCCCACGTCGATAGTGATCGACAAGCCACCCAGTGAGTAACCAACCAGACTCGCCAACAGCGCCGCTCGCACCAGCGTATGGTATTGCCAACCGGTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCCAAGCCGGTGc                                                                                                                                                       <  2:1212229/149‑1 (MQ=255)
  cccACGTCGATAGTGATCGACAAGCCCCCCAGTGAGTAACCAAACAGACTCGCAAACAGCGCCGGACGCACCAGCGGAGGGTATTGCCCACGGTAAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCa                                                                                                                                                     <  2:1046070/150‑1 (MQ=255)
  cccACGTCGATAGTGATCGACAAGCCACCCAGTGAGTAACCAAACAGACTCGCCAACAGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCa                                                                                                                                                     <  2:493487/150‑1 (MQ=255)
   ccACGTCGATAGTGATCGACAAGCCACCCAGTGAGTAACCAAACAGACTCGCCAACAGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCat                                                                                                                                                    >  2:426714/1‑149 (MQ=255)
     aCGTCGATAGTGATCGACAAGCCACCCAGTGAGTAACCAAACAGACTCGCCAACAGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATc                                                                                                                                                  <  1:653880/150‑1 (MQ=255)
     aCGTCGATAGTGATCGACAAGCCACCCAGTGAGTAACCAAACAGACTCGCCAACAGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATc                                                                                                                                                  <  1:1228470/150‑1 (MQ=255)
      cGTCGATAGTGATCGACAAGCCACCCAGTGAGTAACCAAACAGACTCGCCAACAGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATca                                                                                                                                                 >  1:1501051/1‑150 (MQ=255)
         cGATAGTGATCGACAAGCGACCCAGGGAGTAACCAAACAGACTCGCCAACAGCGCCGGACGCACCAGCGGATGGTATTGCCCCCGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATcagca                                                                                                                                              <  2:140307/150‑1 (MQ=255)
           aTAGTGATCGACAAGCCACCCAGTGAGTAACCAAACAGACTCGCCAACAGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAgg                                                                                                                                            <  1:879832/150‑1 (MQ=255)
           aTAGTGATCGACAAGCCACCCAGTGAGTAACCAAACAGACTCGCCAACAGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAgg                                                                                                                                            <  1:1645058/150‑1 (MQ=255)
              gTGATCGACAAGCCCCCCAGTGAGTAACCAAACAGACTCGCCAACAGCGCCGGCCGCACCAGCGGATGGTATTGCCCCCGGTTAAAGACGACTACCGCCCACGCCAGCGCCCAGCCGCCACCGGCAAAGCCGGTGCCAAACAGCAGGTCa                                                                                                                                         >  1:985749/1‑150 (MQ=255)
              gTGATCGACAAGCCACCCAGTGAGTAACCAAACAGACTCGCCAACAGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCGCCGGCAAAGCCGGGGCCAATCGGCAGGTCa                                                                                                                                         >  2:1166776/1‑150 (MQ=255)
                   cGACAAGCCACCCAGTGAGTAACCAAACAGACTCGCCAACAGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGCCGTATACCGCCCACGCCCGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAAcgc                                                                                                                                    >  2:211605/1‑150 (MQ=255)
                   cGACAAGCCACCCAGTGAGTAACCAAACAGACTCGCCAACAGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAACcgc                                                                                                                                    <  2:1203576/150‑1 (MQ=255)
                    gACAAGCCACCCAGTGAGTAACCAAACAGACTCGCCAACAGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAAcgcg                                                                                                                                   >  1:1820590/1‑150 (MQ=255)
                    gACAAGCCACCCAGTGAGTAACCAAACAGACTCGCCAACAGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAAcgcg                                                                                                                                   <  2:422138/150‑1 (MQ=255)
                     aCAAGCCACCCAGTGAGTAACCAAACAGACTCGCCAACAGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAAcgcgc                                                                                                                                  >  1:569219/1‑149 (MQ=255)
                      cAAGCCACCCAGTGAGTAACCAAACAGACTCGCCAACAGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCCGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGag                                                                                                                                 >  1:348902/1‑149 (MQ=255)
                       aaGCCACCCAGTGAGTAACCAAACAGACTCGCCAACAGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATc                                                                                                                                >  2:830759/1‑150 (MQ=255)
                        aGCCACCCAGTGAGTAACCAAACAGACTCGCCAACAGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATcc                                                                                                                               <  2:1326853/150‑1 (MQ=255)
                        aGCCACCCAGTGAGTAACCAAACAGACTCGCCAACAGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTCCCAATCAGCAGGTCAAACGCGATcc                                                                                                                               >  2:784043/1‑150 (MQ=255)
                        aGCCACCCAGTGAGAAACCAAACAGACTCGCCAACAGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACAGCCAACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTAAAACGCGATcc                                                                                                                               >  1:1144560/1‑150 (MQ=255)
                         gCCACCCAGTGAGTAACCAAACAGACTCGCCAACAGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCCGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATccc                                                                                                                              >  2:538665/1‑149 (MQ=255)
                         gCCACCCAGTGAGTAACCAAACAGACTCGCCAACAGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCca                                                                                                                              >  2:518811/1‑150 (MQ=255)
                          ccACCCAGTGAGTAACCAAACAGACTCGCCAACAGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCcac                                                                                                                             <  2:963240/150‑1 (MQ=255)
                           cACCCAGTGAGTAACCAAACAGACTCGCCAACAGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCcacc                                                                                                                            >  2:1811064/1‑149 (MQ=255)
                               cAGTGAGTACCCACACAGACTCGCCAACAGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTACAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGcc                                                                                                                        <  2:1849424/150‑1 (MQ=255)
                               cAGTGAGTAACCAAACAGACTCGCCAACAGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAATCGCGATCCACACGcc                                                                                                                        <  1:856651/150‑1 (MQ=255)
                               cAGTGAGTAACCAAACAGACTCGCCAACAGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGcc                                                                                                                        >  2:211113/1‑150 (MQ=255)
                                 gTGAGTAACCAAACAGACTCGCCAACAGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGATAAAGACGTATACCGCCCACGCCGGCGCCCAGCCGCCACAGGCACAGCCGGTGCCAATCAGCAGGTCACACGCGATCCACACGcccc                                                                                                                      >  2:742597/1‑150 (MQ=255)
                                   gAGTAACCAAACAGACTCGCCAACAGCGCCGGACGCTCCAGCGGAGGGTATTGCCCACGGTTAAAGACGTATGCCGCCCACGCCAGCGCCCAGCCGCCCCAGGCAACGCCGGTGCCAATCCGCAGGGCAAACGCGATCCACACGCCCCAc                                                                                                                    >  2:1772064/1‑150 (MQ=255)
                                   gAGTAACCAAACAGACTCGCCAACAGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCCCACGCCCCAc                                                                                                                    >  2:1607158/1‑150 (MQ=255)
                                    aGTAACCAAACAGACTCGCCAACCGCGCCGGACGCACCAGCGGCTGGTATTGCCCCCGGTTAAAGCCGTAAACCGCCCACGCCAGCGCCCCGCCGCCCCCGGCAAAGCCGGTGCCAAACAGCCGGTCAACCGCGAGCCCCGCGCCCCCCg                                                                                                                   >  1:36239/1‑150 (MQ=255)
                                     gTAACCAAACAGACTCGCCAACAGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACgg                                                                                                                  >  2:239353/1‑150 (MQ=255)
                                      tAACCAAACAGCCTCGCCAACAGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACggg                                                                                                                 <  1:480244/150‑1 (MQ=255)
                                      tAACCAAACAGACTCGCCAACAGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACggg                                                                                                                 <  1:907584/150‑1 (MQ=255)
                                         ccAAACAGCCTCGCCAACAGCGCCGGACGCACAAGCGGATGGTCTTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGAAATGCTGGTTACATTCAGCTGGTCAAACGAGATCTACACGCCCCACGGGAAg                                                                                                              <  2:922805/150‑1 (MQ=255)
                                          cAAACAGACTCGCCAACAGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCCGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACCCGCCCCACGGGAAgc                                                                                                             >  1:1808923/1‑150 (MQ=255)
                                          cAAACAGACTCGCCAACAGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAgc                                                                                                             >  1:1816086/1‑150 (MQ=255)
                                          cAAACAGACTCGCCAACAGCGCCGGACCCACCAGCGGATGGTATTGCCCACGGTTAAAGACCTAAACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCCAACGCGATCCACCCGCCGCACGGGAAgc                                                                                                             >  2:1078296/1‑150 (MQ=255)
                                               aGACTCGCCAACAGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCACTCCGCAGGGCACCCGCGATCCCCCCGCCCCACGGGCCGacgccg                                                                                                        >  1:1013303/1‑150 (MQ=255)
                                                gACTCGCCAACAGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCCGCGCCCCGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGCTCCACACGCCCCCCGGGAAGCCGCCGt                                                                                                       >  2:1800956/1‑150 (MQ=255)
                                                  cTCGCCAACAGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTc                                                                                                     <  1:819360/150‑1 (MQ=255)
                                                   tCGCCAACCGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAATACGTTGACAATAAGCTGGTCAAACGAGATCAACACGCCCCACGGGAAGCCGCCGTTCa                                                                                                    <  2:740264/150‑1 (MQ=255)
                                                     gCCAACAGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCGGCAGGTCAAACGCGATCCACACGCCCCACGGGCAGCCGCCGTTCAgg                                                                                                  >  2:1709015/1‑150 (MQ=255)
                                                     gCCAACAGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAgg                                                                                                  <  1:830759/150‑1 (MQ=255)
                                                     cccAACAGCGCCGGACGCACCAGCGGATGGTTTTGCCCACGGTTAAACACGGATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAATGGCGGTGACTTTCAGCAGGTCAAACGCGATCAACACGCCCCACGGGAAGCCGCCGTTCAgg                                                                                                  <  2:136884/149‑1 (MQ=255)
                                                      ccAACAGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGTTTGTGACAATCAGCAGGTCAAACGAGATCTACACGCCCCACGGGAAGCCGCCGTTCAGGt                                                                                                 <  2:1089039/150‑1 (MQ=255)
                                                             cgcCGGACGCAGCAGCGGATGGTATTGCCCAAGGTTAAAGACGTTTACCGCCCACGCCAGCGCCCAGCCGCCACAGGCATAGTCGGTGCCAATAAGCAGGTCAAACGAGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGAc                                                                                          <  2:1338108/150‑1 (MQ=255)
                                                             ccccGGACGCACCAGCGGATGGGTTTGCCCACGGTTAACGACGTATACCGCCCACGCCACCGCCCAGCCGCCACAGGATTTTTTTTTTAATTTTAGGTGGTCAACCGAGATCTACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGAc                                                                                          <  2:1124224/148‑1 (MQ=255)
                                                              gcCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACa                                                                                         >  2:892416/1‑150 (MQ=255)
                                                              gcCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACa                                                                                         >  2:1095972/1‑150 (MQ=255)
                                                                  gACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCATATTTTTTTAAATTTAGGTGGTCAAACGAGATCTACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGc                                                                                     <  2:31733/150‑1 (MQ=255)
                                                                    cGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCAACAGGCCAAGCCGGTGCCAATCAGCTGGTCAAGCGGTCTCCCCACGCCCCCCGGGCAGCCGCCGTGCTGGTCAGCGCCagcggcc                                                                                   >  2:335897/1‑145 (MQ=255)
                                                                       aCCAGCGGATGGTATTGCCCACGGTTAACGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCATAGCCGGTGCCTATCAGCTGGTCATACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGa                                                                                <  2:903026/150‑1 (MQ=255)
                                                                         cAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCGACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGAGCCACACGCCCCACGGTAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGAcc                                                                              <  1:886433/150‑1 (MQ=255)
                                                                          aGCGGATGGTATTGCCCAAGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCATAGCTGGTTCCTTTCAGGTGGTCAAACGCGATCAACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCg                                                                             <  2:1358554/150‑1 (MQ=255)
                                                                           gCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGa                                                                            <  1:675513/150‑1 (MQ=255)
                                                                            cGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGaa                                                                           >  2:1624322/1‑150 (MQ=255)
                                                                             ggATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAATGTATTTTAATTATAGGTGGTCAAACGAGATCTACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGAac                                                                          <  2:1220016/150‑1 (MQ=255)
                                                                              gATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACCGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGAacc                                                                         >  1:1610304/1‑149 (MQ=255)
                                                                                tGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCATTGTTGTTTCAATTTAGGTGGTCAAACGAGATCTACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGAACAcc                                                                       <  2:1377474/150‑1 (MQ=255)
                                                                                tGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAACCGCGCTCCACACGCCCCCCGGGACGCCGCCGTTCAGGTCAGAGACAGAGCCCAGCCCGCACAcc                                                                       >  1:920404/1‑150 (MQ=255)
                                                                                tGGTATTGCCCACGGTTAAAGACGTAAACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGCTGACAATCAGCAGGTCAAACGCGCTACACACGCCCCACGGGACGCCTCCGGTCAGGACAGAGACAGAGCCCCGCCCGCAcccc                                                                       >  1:361563/1‑150 (MQ=255)
                                                                                  gTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCCCAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGCGCCCAGACCGAACACCAg                                                                     >  1:1234585/1‑150 (MQ=255)
                                                                                  gTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCAGAGGCAAAGACGGTGCCAATCAGCAGGGCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGAACACCAg                                                                     <  2:471585/150‑1 (MQ=255)
                                                                                     ttGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGACTGTTCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGAACACCAGACg                                                                  <  2:1435167/150‑1 (MQ=255)
                                                                                        cccACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGAACACCAGACGCtt                                                               >  2:2162/1‑150 (MQ=255)
                                                                                                tAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGACAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGAACACCAGACGCTTCACAATCa                                                       <  2:161709/150‑1 (MQ=255)
                                                                                                  aaGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCTGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGAACACCAGACGCTTCACAATCAgg                                                     <  2:380379/150‑1 (MQ=255)
                                                                                                     aCGTATACCGCCCCCGCCAGCGCCCCGCCGCCACAGGCAAAGACGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGAACACCAGACGCTTCACAATCAGGAGc                                                  <  2:101503/150‑1 (MQ=255)
                                                                                                       gTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGAACACCAGACGCTTCACAATCAGGAGCAt                                                <  1:1709628/150‑1 (MQ=255)
                                                                                                         ataCCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGAACACCAGACGCTTCACAATCAGGAGCATAc                                              <  2:1501051/150‑1 (MQ=255)
                                                                                                          taCCGCCCACGCCAGCGCCCAGCCGCCACCGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGCGCCCAGCCCGCACCCCAGCCGCTTCCCAAGCAGGAGCAGACa                                             >  2:1532310/1‑150 (MQ=255)
                                                                                                          taCCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACCCGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGAACACCCGACGCTTCACACTCAGGAGCATAcc                                             >  1:1400723/1‑149 (MQ=255)
                                                                                                             cGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGAACCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGAACACCAGCCGCTTCACAATCAGTAGCATACAGAt                                          >  1:1284715/1‑150 (MQ=255)
                                                                                                                ccACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACCCGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCCGACCGAACACCAGACGCTTCACAATCAGGAGCATACAGATGAc                                       >  1:810156/1‑150 (MQ=255)
                                                                                                                 cACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGAACACCAGACGCTTCACAATCAGGAGCATACAGATGACg                                      >  2:755639/1‑150 (MQ=255)
                                                                                                                       aGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGAACACCAGACGCTTCACAATCAGGAGCATACAGATGACGATTAAc                                <  1:211113/150‑1 (MQ=255)
                                                                                                                       aGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGAACACCAGACGCTTCACAATCAGGAGCATACAGATGACGATTAAc                                >  2:36237/1‑150 (MQ=255)
                                                                                                                         cgcCCAGCCGCCACAGGCAAAGCCGGTGACAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGAACACCAGACGCTTCACAATCAGGAGCATACAGATGACGATTAACgg                              >  2:1700942/1‑150 (MQ=255)
                                                                                                                            ccAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGAACACCAGACGCTTCACAATCAGGAGCATACAGATGACGATTAACGGTcc                           >  1:1548367/1‑150 (MQ=255)
                                                                                                                             cAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCCCACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGCCCGAGCGCCAGACGCTTCAGAATCAGGTGCATCCAGATGTCGATGAACGGTCCa                          >  2:1243796/1‑150 (MQ=255)
                                                                                                                             cAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGAACACCAGACGCTTCACAATCAGGAGCATACAGATGACGCTTAACGGTCCa                          >  2:1374657/1‑150 (MQ=255)
                                                                                                                             cAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGAACACCAGACGCTTCACAATCAGGAGCATACAGATGACGATTAACGGTCCa                          >  1:111480/1‑150 (MQ=255)
                                                                                                                               gccgccACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGAACACCAGACGCTTCACAATCAGGAGCATACAGATGACGATTAACGGTCCaaa                        >  2:204497/1‑150 (MQ=255)
                                                                                                                                 cgccACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGAACACCAGACGCTTCACAATCAGGAGCATACAGATGACGATTAACGGTCCaaaaa                      <  2:1820590/150‑1 (MQ=255)
                                                                                                                                  gccACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGACCACCAGACGCGTCCCAATCAGGAGCATAGAGATGACGATTAACGGTCCAAAAAt                     >  2:325737/1‑150 (MQ=255)
                                                                                                                                  gccACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGAACACCAGACGCTTCCCAATCAGGAGCATACAGATGACGATTAACGGTCCAAAAAt                     >  1:388068/1‑150 (MQ=255)
                                                                                                                                  gccACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGAACACCAGACGCTTCACAATCAGGAGCATACAGATGACGATTAACGGTCCAAAAAt                     <  1:1095972/150‑1 (MQ=255)
                                                                                                                                    cacaGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGAACACCAGACGCTTCACAATCAGGAGCATACAGATGACGATTAACGGTCCAAAAATCa                   <  1:892416/150‑1 (MQ=255)
                                                                                                                                      caGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGAACACCAGACGCTTCACAATCAGGAGCATACAGATGACGATTAACGGTCCAAAAATCATg                 <  1:1757505/150‑1 (MQ=255)
                                                                                                                                      caGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGAACACCAGACGCTTCACAATCAGGAGCATACAGATGACGATTAACGGTCCAAAAATCATg                 <  1:1057428/150‑1 (MQ=255)
                                                                                                                                        ggCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGGCCGACCACCAGCCGCTTCACAATCAGGAGCATACAGATGACGATTACCGGTCCAAACATCATGAc               >  2:1640300/1‑150 (MQ=255)
                                                                                                                                          cAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGCAGCCGCCGTTCAGGTCAGAGACAGCGCCCAGACCGAACCCCAGACGCTTAACATTCAGGAGCATACAGATGACGATTAACGGTCCAAAACTGATGCCCg             >  1:1387802/1‑150 (MQ=255)
                                                                                                                                           aaaGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGAACACCAGACGCTTCACAATCAGGAGCATACAGATGACGATTAACGGTCCAAAAATCATGACCgg            >  2:1267277/1‑150 (MQ=255)
                                                                                                                                             agtcGGTGCAAATCAGCAGGTAATACGCGATGAACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGAACTCCAGACGCTTCACAATCAGGAGCATACAGAAGACGATTAACGGTCCAAAAATCATGACCGGtt          <  1:669217/147‑1 (MQ=255)
                                                                                                                                             aGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGAACACCAGACGCTTCACAATCAGGAGCATACAGATGACGATTAACGGTCCAAAAATCATGACCGGtt          >  1:1665033/1‑150 (MQ=255)
                                                                                                                                                     ccAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGAACACCAGACGCTTCACAATCAGGAGCATACAGATGACGATTAACGGTCCAAAAATCATGACCGGTTTACtgatg  <  1:1799448/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GACCCACGTCGATAGTGATCGACAAGCCACCCAGTGAGTAACCAAACAGACTCGCCAACAGCGCCGGACGCACCAGCGGATGGTATTGCCCACGGTTAAAGACGTATACCGCCCACGCCAGCGCCCAGCCGCCACAGGCAAAGCCGGTGCCAATCAGCAGGTCAAACGCGATCCACACGCCCCACGGGAAGCCGCCGTTCAGGTCAGAGACAGAGCCCAGACCGAACACCAGACGCTTCACAATCAGGAGCATACAGATGACGATTAACGGTCCAAAAATCATGACCGGTTTACTGATG  >  NC_000913/3143831‑3144129

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